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SUMMARY 


Detailed measurements of the mean and turbulence characteristics of 
the surface wind measured near the Atlantic coast at Wallops Island, 
Virginia, are presented In this comprehensive final report. This final 
report consists of five articles and a detailed report of the 
experimental results for the U.S. Department of Energy under contract 
no. AC 06-79 ET 23007. Earlier reports dealing with the 
Instrumentation, calibration, data handling, data analysis and early 
results are not Included In this final report. Detailed information on 
these items are available as reports from the National Technical ! 

Information Service under NASA-CR-62093, 1973, NASA-CR-1 37456, 1974, 
NASA-CR-1 40586, 1974, and NASA-CR-1 41 422, 1977. 

* 

The experimental data were acquired from a 76-meter tall [ 

/ 

Instrumented micrometeorological tower. Mean wind and turbulence 
measurements were made with two types of instrumentation consisting of 
cup-vane and temperature probes, primarily used for mean profile 
measurements of velocity and temperature respectively. The second 
system, a hot film and thermocouple system, was used for measurement of 
turbulence variances and covariances and spectra. The cup-vane system 
was used to acquire data from all wind directions, while the hot-film 
system was primarily used for turbulence measurements from the two 
prevailing wind directions, south and north-west. 

The micrometeorological tower is a self-standing non-guyed tower 
with five working platforms at 15.2m (50 ft.) intervals, with cup-vane 
and aspirated temperature probes mounted at each platform. This 
instrumentation system is used by N.A.S.A. in conjunction with its 
sounding rocket launching operations. With this system data could be 


acquired without interruption for about eight hours at a rate of one 
sample each two seconds. A total of 195 digital tapes were generated 
with this system during the period of July 1974 and December 1978. 
Approximately 300 data records were analyzed, each varying between 43 
and 85 minutes in duration. 

Six three-dimensional split-film anemometers (T.S. I. -1080) were 
used for the turbulence measurements. These Instruments were mounted on 
1.8m booms at the same levels as the cup-vane instruments and also at 
the 9m (30 ft.) level. The electronics for this system as well as the 
data-acquisition and data-handling system were located in an 
instrumentation trailer parked at the base of the tower. This sytem was 
designed to handle outputs from six split-film anemometer systems, 
sampled at a rate of 200 samples per second for a period of 
approximately one hour. A total of 32 one-hour data runs were acquired 
with the hot-film system between July 1976 and December 1978 and 
analyzed for variances, covariances and spectra. 

All data were acquired under predominantly moderately strong wind 
conditions with an hourly mean wind speed between 10 m/s and 20 m/s at 
the 76.2 (250 ft.) level. Under these conditions, the atmosphere near 
the surface cannot be considered neutrally stable, and convective 
effects must be included in the analysis of the data. For on-shore 
winds, the surface flow is more complicated as the result of the 
development of an internal boundary layer. The latter is initiated as 
the air crossing the beach from the water generally experiences a change 
in surface roughness and in surface temperature. The internal boundary 
layer has a height between 15m and 30m at the tower location depending 
on the development distance from the beach which varies with wind 


direction. Also for southerly winds during the day time, the warmer air 
flows over the cooler water, allowing the formation of a surface»based 
inversion of variable depth. Under these conditions a low-altitude 
maximum velocity (surface jet), occasionally below the highest 
observation level of 76m, has been observed. Under extreme stable 
conditions at hourly mean velocities in excess of 10 m/s at the highest 
observation level, the turbulence has been observed to vanish 
completely. Moreover, low-frequency internal gravity waves have been 
observed to co-exist with the turbulence. 

For continental winds under convective conditions, usually 
encountered during the day time, the planetary boundary layer (P.B.L.) 
is usually characterized by a strong upward heat flux and by strong 
vertical mixing due to the positive buoyancy forces. As the result of 
this intense mixing process a relative sharp inversion layer separates 
the mixed layer from the free atmosphere. This inversion delineates the 
depth of the P.B.L. Under uniform surface conditions the surface flow 
can be modeled adequately by the Morin-Obukhov (M.O.) similarity theory. 

Westerly winds approaching the Wallops site pass in succession over 
the mainland, the Chesapeake Bay, the "Delmarva" peninsula and the tidal 
marsh. For the last 100m-300m the air travels over the island before it 
arrives at the tower location. As the result of the nonuniformity of 
the surface roughness and the surface temperature encountered by these 
westerly winds, the appl icability of the M.O. theory is questionable. 

As the air flows over water the upward heat transfer is interrupted and 
the downward heat flux from the inversion at the edge of the P.B.L. « 
penetrates down to the surface. This downward heat flux is often 
observed at the tower site under convective conditions, while at the 


same time the temperature gradient would suggest an upward heat flux. 

More complications of the surface flow may be encountered when 
clouds are present and condensation of water vapor occurs within the 
P.B.L. Under these conditions latent heat is released and an inversion 
layer develops just below the cloud base which is responsible for the 
suppression of the turbulence below this inversion. Also for westerly 
winds between 300° and 330°, the instantaneous wind direction 
experiences frequently large fluctuations toward the north. These 
direction fluctuations are not normally distributed, providing an 
explanation for the unusual high lateral turbulence intensities and the 
non-zero lateral horizontal shear component 

The experimental results of this research program clearly indicate 
that there is no single and no simple flow model available to describe 
the mean and turbulent flow near the surface under moderately strong 
wind conditions at the Wallops Island site. The presence or absence of 
appreciable low-frequency velocity fluctuations causes the statistical 
parameters describing the turbulence to vary a great deal. The non- 
uniform surface conditions and the presence of large-amplitude low- 
frequency fluctuations, either In the convective boundary layer or in 
the on-shore winds in the form of interval gravity waves, cause the 
observed surface flow to be quite different from the M.O. boundary layer 
flow model . 
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AN INSTRUMENTATION SYSTEM FOR THE MEASUREMENT OF AT- 
MOSPHERIC TURBULENCE 
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and Mechanics, Blacksburg, Virginia 24061 (USA) 

(Received February 10, 1976) 


Summary 

A method is described for the detailed multipoint measurement of the turbulent veloc* 
ity components and temperature in the atmospheric surface layer. 

The TSI three-(^{Mensional split-film anemometer system is used as sensor and a special 
designed data<ucquisition and data-handling system Is used for the sampling, digitizing and 
storing of the data on digital magnetic tapes. Further processing of the data is performed 
in a digital computer and programs are developed to allow efficient calculation of the sta- 
tistical quantities such as, mean values, covariances, power- and cross-spectral density 
functions. At the present, six anemometer systems have been mounted on the 30-, 50-, 100% 
150-, 200- and 250-foot levels of the 250-foot Meteorological Tower at NASA, Wallops 
Flight Center. The tower is located at Wallops Island which is one of the barrier islands 
protecting the so-called *'Delmarva'* petdnsula from the Atlantic Ocean. The island, which 
consists of a narrow strip of dunes approximately three meters above sea level, is situated 
in a northeast— southwest direction and is separated from the mainland by a three-mile 
wide tidal swamp. 


Introduction 

The aerodynamic properties of the lower part of the atmospheric boundary 
layer up to a height of approximately 500 feet, are usually measured with in- 
struments placed on masts, meteorological towers or other tall structures. For 
measurements at higher altitudes, instruments placed on meteorological bal- 
loons or low-flying aircraft can be used. During the last two decades a number 
of research programs dealing with these measurements and interpretation of 
the aerodynamic properties of the atmospheric boundary layer have been 
carried out. A discussion of some of the tower measurements can be found in 
the literature and a comprehensive review of these studies has been given by 
Panofsky [1] , Busch [2] , Haugen [3] and Panofsky [4] . Except for a few 
cases, most of these research programs have been limited in scope because of 
the lack of simultaneous and multi-point velocity and temperature measure- 
ments over a wide frequency range. For example, in many cases either the 

^ At present graduate student at Johns Hopkins University, Materials Science and Mechanics 
Department, Baltimore, Maryland. 
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vertical velocity componfiftl or the fluctuating temperature are missing or 
measurements are made with instruments whose frequency response is limited. 

The objective of this NASA-supported research program is to develop a 
system which will be able to obtain a complete and an extensive set of multi- 
point measurements of the three velocity components and the temperature 
in the lower atmosphere. These measurements will allow the accurate deter- 
mination of mean values and statistical information of the low level atmo- 
spheric turbulence. 

Instrumentation 

The total vector anemometer 1080 system manufactured by Thermo-Sys- 
tems, Inc. was chosen for this research program since it seems to have the ad- 
vantage of small physical size, fast response and high sensitivity over a wide 
range of velocities. In addition, the probe consists of three split-film sensors 
as shown in Fig.l, which allows the determination of magnitude and direc- 
tion of the instantaneous velocity vector. 



Fig.l. Schematic diagram of the sensor oriented coordinate system, ABC, and the probe 
oriented coordinate system, XYZ. 


The three sensor rods are mounted mutually perpendicular to form a Car- 
tesian coordinate system from now on referred to as the sensor-oriented co- 
ordinate system. The active sensing element consists of a thin platinum film 
of approximately 1000 Ingstrdro. thickness and 0.080 inch long, which is 
placed on a 0.2 inch long and 0.006 inch diameter cylindrical quartz rod. 

A thin coat of quartz over the film is used to provide protection from the 
environment. The platinum film on each rod consists of two segments. 
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separated from each other by two longitudinal splits 180 degrees aparf**. The 
orientation of the splits on each sensor rod is such that it allows the detenni* 
nation of the octant in which the instantaneous velocity vector is located. For 
minimum support interference, the three senser rods are mounted at only one 
end to a thin supporting structure in which the electrical leads are embedded. 
The active elements on each rod are electrically heated to the same constant 
temperature by separate anemometer circuits, A copper— constantan thermo- 
couple is mounted close to the sensors and measures the ambient temperature 
with a maximum frequency response of 25 Hz. The manufacturer’s perfor- 
mance specifications list a response up to 1 kHz depending on the cable length 
used. For tower mounting of the probes, cable len^hs of approximately 350 
feet and 200 feet were used. The three sensors and the tl^ermocouple are pro- 
tected, while not in use, by a shield which can be placed over the sensors and 
retracted pneumatically. An air stream is continuously supplied through the 
shield, in order to keep the sensors free of contamination and moisture while 
not in operatiorji, 

The data acquisition and data handling system was designed to handle out- 
put data from one to thirteen anemometer systems, sampled at a rate of 200 
samples per second, for a duration of approximately 45 minutes. All necessary 
equipment and controls are fixed in an instrumentation trailer placed under 
the meteorological tower. The system consists of two main parts: (a) the mul- 
tiplexing and analog recording system, schematically shown in Fig,2 and (b) 
the demultiplexing, digitizing and digital recording system, shown In Fig.3. 


7 V.C,0.‘» 



Fig.2. Multiplexing and analog recording syitem. 


^Patented by Themio-Systenu, Inc. 
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OtSCRIMINATORS 



Fig.3. Demultiplexing, digitizing and digital recording system. 

Each anemometer provides seven analog voltages, one from the thermo- 
couple circuit and six from the anemometer circuits, all in the range of 0—5 
volts, All voltages are frequency modulated by voltage-controlled oscillators, 
each with a different center frequency. For the output of each anemometer 
there is one set of voltage-controlled oscillators with center frequencies of 
8, 12, 16, 20, 24, 28 and 32 kHz, respectively. All voltage-controlled oscU; 
lators operate at a deviation from their center frequency of ±1 kHz for a 
m axim um voltage input of ±5 volts. The seven frequency-modulated signals, 
together with a reference signal of 100 kHz are summed by a summing am- 
plifier to produce one single multiplexed signal of approximately 3 volts 
peak-to-peak. The multiplexed signals from the six anemometers are simul- 
taneously recorded on six separate channels of a 14-channel analog tape re- 
corder. The time-of-day is recorded on a separate channel of the analog tape- 
recorder and serves as a reference for the recorded data. 

At a later time, the analog tape is replayed and the multiplexed signal from 
the desired channel is demultiplexed to its seven components after passage 
through seven discriminators. A reference discriminator senses any deviation 
in the reference frequency of 100 kHz caused by fluctuations in the tape 
speed of the analog tape recorder and corrects the output of the seven dis- 
criminators accordingly. Consequently, the six outputs corresponding to the 
six anemometer bridge voltages are each passed through a 100 Hz low-pass 
filter. This is necessary to avoid aliasing (distortion) of the velocity spectra 
at firequencies below 100 Hz from contributions of frequencies above 100 Hz. 
Next, the seven analog voltages are passed through the analog-to-digital con- 
verter where the analog signals are sampled at a rate of 200 samples per second 
and each sample is converted to an 16-bit (plus sign) digital word. The total 
scan, settling and conversion time for the sampling of the seven voltages is 50 
microseconds, which is approximately 1% of the time interval between two 
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scans. The recorded time-of-day is converted to a 32-bit binary coded decimal 
signal by a serial-to-parallel converter which is interfaced with the controller 
coinputer by a custom-designed interface. The computer (DEC Model PDP 
11/20) controls the multiplexing, analog-to-digital conversion as well as the 
digital recording. The computer is programmable in DEC PAL-11 assembly 
language and has an 8 K memory of 16-bit-words. Access to the computer is 
obtained throu^ a teletypewriter. 

The data conversion starts at a time-of-day prescribed by the operator, and 
the: analog-to-digital converter performs successive scans and conversion of 
the seven analog voltages into two-byte (16 bit) words at a rate of one scan 
every 5 milliseconds. These words are stored in one of the four buffers of the 
computer which in turn dumps the data on a 9-track digital magnetic tape. Each 
buffer has a capacity of 209 scans which represents 1.05 seconds of data. 

A tot?i of 2500 records make up a data sample, extending over a time period 
of approximately 45 minutes. 

Operation 

The total heat flux for a split-film sensor, where both films are operating 
at the same temperature is given by 

+ llRfi )liTf-T„) = C + DipUef' (1) 

where Ii and I 2 are the heating currents of films 1 and 2, respectively. Rfi 
and Rf'i are the hot film resistances, Tf is the film operating temperature, Tg 
is tlie ambient temperature, p is the air density, Ue is the effective cooling 
velocity which is a function of the velocity magnitude and the sensor yaw angle, 
(^, and n is the exponent in the heat transfer equation. Introducing the relation 
for: the effective cooling velocity as suggested by Hinze [5] and relating the 
velocity to standard pressure, Pg - 29.92 inches of mercury and standard tem- 
perature Tg - 530 degrees Rankin, the heat flux equation becomes 

{l\Rfi + l\Rfi )l(Tf-Ta) » C + DU^ (cos* 0 + fe* sin* 0)«/* , (2) 

Where Ug is the velocity related to standard conditions and fe is an experimental- 
ly determined coefficient of the order of 0.1. 

Introducing the bridge voltages £1 and ^- ^ox films 1 and 2, respectively, equa- 
tion (2) can be modified to 

Q/A T = (^i CFiE\ + K 2 CF-,E\ )HTf-Ta) = Dl72(cos* 0 + fe* sin* 0)"/* , (3) 

where Q is a measure of the total heat transfer, CF\ and CF^ are temperature 
correction factors for films 1 and 2, respectively, to account for cable and 
ambient temperature changes as discussed and defined by Tieleman and 
Tavoularis [6] . The constants and X 2 should ideally be the ratio of the 
film resistance and the square of the corresponding bridge resistance. How- 
ever, these constants also correct for area and resistance differences of the 
two films and must be determined experimentally. In addition, for all veloc- 
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ities of interest in this study, the constant C may be eliminated by adjustment 
ofthe constants ICi andKi. 

Assuming that the three sensors A, B and C are mounted mutually perpen- 
dicular, the following trigonometric relation holds for the sensor-yaw an^es. 

sin* + sin* + sin* * 1* (4) 

The standard effective-cooling velocities are related to the magnitude of 
the standard velocity Ug and the sensor-yaw angles, i^y the ^ixpression 

(cos* <t>a + sin* 0a), (5) 

where the subscript a refers to either one of the sensors A, B or C. Assuming 
that identical values of kg, exist for all three sensors and using expressions (4) 
and (5), the magnitude of the standard velocity can be expressed in terms of 
the effective-cooling velocities of the three sensors as follows, 

VI “ )/*av. (6) 

Here the constant kgu is of the order of 2.1 and needs to be determined throu^ 
calibration. 

The sensor-yaw angles can be obtained from expressions (4) and (5) and 
be of the following form 

|0al = arcsin[{l-(CaUes„/t/,)*}/(l-fe„)]'/*, (7) 

where the coefficients are of the order of unity but due to variations in 
heat transfer with sensor-pitch angle, (pa, these coefficients may deviate a 
small amount from unity. In reference to Fig.l, it can be seen that the sign 
of the sensor yaw angles, 0 a, can be determined by comparing the output 
voltages of sensor C. In case 0 a is zero degrees, the velocity vector is located 
in a plane formed by the sensors B and C which is similar to the plane formed 
by the two splits of sensor C. In this case, the ratio of the output voltages of 
sensor C is defined as 

Eic/Eic ~ lie- ( 8 ) 

In general, when the ratio of the film voltages, Eic/Eic is larger than Rq , 
the sensor-yaw angle, 0a, is considered positive and if the above ratio is less 
than iZc, 0A is considered to be a negative angle. Similarly, the magnitude 
of the ratios Ei/^/Eza will determine the sign of the sensor-yaw 

angles 0B and' 0c, respectively. 

The actual velocity components in the sensor-oriented coordinate system 
can now be calculated since the magnitude of the velocity and its direction 
with respect to the sensor-oriented coordinate system is known. The actual 
velocity can be obtained ftom the standard vdocity using the continuity 
equation and the equation of state for a perfect gas as 

U^(P,/P)(T/T,)U„ (9) 


where P and T are the ambient pressure and ambient absolute temperature of 
the air at the time the measurements were made. 

The correction factors CFi and CFj in expression (3) need to be intro- 
duced to account for cable temperature changes as well as ambient tempera- 
ture changes. According to Tieleman et al. [7] , the values for diese correc- 
tion factors for the six different bridge voltages can be calculated as follows 

CF^i = [1 + (^RcJRU^ -Ta) . (10) 

Again, the subscript at refers to either sensor A, B or C and the subscript i > 

1,2 indicates the appropriate film. The zero-superscribed quantities refer to 
the calibration conditions in the TSI laboratory. ARc represents the change 
in cable resistance, Rf represent the hot film resistances. The second factor 
in (10) represents the correction for ambient temperature and film tempera- 
ture. 

Tieleman et al. [7] tested the accuracy of the three-dimensional split-film 
anemometer in the low-speed wind tunnel at VPI and SU for a wide range of 
velocities and angles of attack. In order to improve the accuracy, changes 
were made in the so-called “improved accuracy” as suggested by Olin and 
Kiland [8] . hi addition, a temperature correcfion was introduced to correct 
the output voltages for variations in heat transfer and cable resistance as a 
result of variations in ambient temperature. These changes were not sufficient 
since new tests indicated rarely a percentage error of less than 5 percent for 
the calculation of the velocity components. This high percentage of error was 
due to a number of problems which were prominent over different ranges of 
probe-yaw angles. The heat transfer is not constant when, the sensor rod is 
rotated about its axis which results in errors in the evaluation of the cooling 
velocities, total velocity and sensor-yaw angles. Secondly, large errors are made 
in estimation of the sensor-yaw angle, (p, when the absolute value of this 
angle is less than 10 degrees which occurs when the probe-yaw angle, is 
approximately ±40 degrees. Errors due to convection of heat from an up- 
stream sensor to a sensor located directly downstream results in a lower es- 
timate for the cooling velocity of the downstream sensor. This occurs when 
the probe-yaw angle, p, is ±90 degrees. For certain probe-pitch angles, the 
thermocouple senses an ambient temperature which is too high as a result of 
convection of heat from an upstream sensor. In addition, a different method 
for the evaluation of the velocity components as suggested by Olin and 
Kiland [8] was tried with slightly better results. This method was later re- 
jected because the required calibration procedures were too lengthy. From 
all these tests it was concluded that the TSI three-dimensional split-film 
probe cannot measure velocity and direction over the full 360’’ solid angle 
in three-dimensional flow fields within the accuracy as was claimed by the 
manufacturer. It was finally concluded that this instrumentation can only be 
used with the mean velocity closely, parallel to the probe axis if an accuracy 
of less than 5 percent is required. In order to satisfy this requirement, it was 
decided to use the three-dimensional split-film probes in conjunction with a 
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lotatjjig mechanism which allows the probe to be rotated about a vertical 
axis so that the average probe-yaw angle will be close to zero. 

Calibration 

Based on the conclusions reached as the result of extensive testing of the 
probes, new calibration procedures were developed by Tieleman and 
Tavoularis [6] . The calibration information is obtained from wind tunnel 
testing of the three-dimensional split-film probe as follows. First, data is 
taken over a range of velocities between 3.5 and 50 feet per second with each 
of tt»e three sensors perpendicular to the flow. With the probe-yaw angle, |3, 
equal to zero and probe pitch angle, a = —35,26 degrees, sensor A is perpen- 
dicular to the flow. Sensors B and C become perpendicular to the flow when 
the probe pitch angle a, equals zero and (3 is +39.25 degrees and ^ is —39.25 
degrees, respectively. Additional data are taken for a = 0 degrees and probe- 
yaw angles in the range of /3 » ±50 degrees in steps of 10 degrees. In addition, 
measurements have to be made of the cold film resistance and cable resistances 
in order to determine the correction factors. The various constants that are 
necessary for the calculation of the velocity components are determined from 
the above measurements. First, the heat flux from each sensor is calculated 
when the sensor is perpendicular to the flow according to expression (3), 
using values for the constants Jf , and Ki those suggested by the manufac- 
turer as a first approximation. These constants for two of the three sensors 
are then adjusted so that the heat transfer versus cooling velocity curves col- 
lapse to one single curve, see Fig.4. The next step is to determine the proper 
values for coefficients D and exponent n in equation (3) with the sensor-yaw 
angle, (j>, equal to zero. It turns out that no single set of values for D and n de- 
scribes the calibration curve adequately. The calibration curve has to be divided 
into three ranges with different values for D and n in each range. Once the 
heat flux versus effective cooling velocity relation has been established, the 
standard cooling velocities, Ucs, can now be calculated for the set of data where 
the probe-yaw angle is varied between ±50 degrees. Assuming that an expres- 
sion like (6) exists, an average value of the coefficient kgv can be calculated 
since the magnitude of the standard velocity Ug is measured independently by 
a Pitot static tube during the wind tunnel testing. The final step is to deter- 
mine the calibration constants Ca and ka for expression (7). These constants 
for sensors B and C are obtained by fitting a curve of the form of expression 
(7) to the calibration data. The latter are obtained by plotting the velocity 
ratio Ues/Ug versus the predetermined sensor-yaw angles, as shown in Fig.5. 

The yaw angle for sensor A can be obtained by using the trigonometric rela- 
tion (4) after the values of angles and <pc have been determined. 

As an example of the above calibration procedure and the accuracy with 
which the velocity components can be obtained, the calibration constants for 
the probe mounted on the IQO-foot level of the meteorological tower are given 
in Table 1. 

Figures 6 and 7 show the accuracy with which the velocity components 
can be obtained as a function of the probe-yaw angle, 0. For small probe-yaw 
angles, the velocity components can be determined t^'ithin 5% with respect to 
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TABLE 1 


Cillbrmtioii coniUnU for the TSl No. 11^ 3 An«nr>meter 


S^inior 

Film 

n 

Tf 

^^old 

oCoid 


T? 


K 

k 

c 

R 



"F 

"F 

“F 

Ohm 


"F 

Ohm 






1 

76 

633 

73.5 

43.8 

0.001087 

77 

Y.33 

2.0804 



0 9022 

A 

2 

76 

633 

73.5 

44.5 

0.001074 

77 

7.31 

2.1540 





1 

76 

521 

73.6 

42.7 

0.001254 

77 

7.39 

2.0793 

A 1 9 

1 moi 

1.0574 

B 

2 

76 

521 

73.5 

42.9 

0.001273 

77 

7.26 

2.2088 





1 

76 

516 

73.5 

49.6 

0.001258 

77 

7.37 

2.0280 

A 1A 

1.0309 

0.9821 

C 

2 

76 

516 

73.6 

42.9 

0.001295 

77 

7.28 

2.1748 

v.Xv 


Q/AT > 0.52, D - 0.1299 and n - 0.44 
0.38 < Q/aT < 0.52, D - 0.1723 and n - 0.36 
Q/aT < 0.38, D - 0.1761 and n - 0.34 
kav ■ 2.086B 


is the temperature at which the cold>film resistance, , is measured and af is the 

temperature coefficient of resistivity of the film. 7^ is the temperature at which the resis- 
tance of the <£abie, R^, is measured. K represents the film constants as in equation (3), and 
k the constant in the expression for the effective cooling velocity (2). C is a constant in 
the expression for the sensor-yaw angle (7). R represents the ratio of the two film voltages 
of each sensor at zero, sensor-pitch angle. D and n are a coefficient and exponent in the heat 
flux expression (3) and fcau is a constant in the expression for the determination of the 
standard velocity. 

the magnitude of the velocity vector over a velocity range data are ex- 
pected to be gathered. For probe-yaw angles of about plus and minus 
40 degrees, large errors may occur as a result of the low sensitivity of the 
ratio Ues/Ug with respect to the sensor-yaw angles <^c ^utid respectively. 

The general calibration procedure which was suggested by TSI proved to 
be complicated and inaccurate and as a result this new calibration procedure 
had to be developed. This calibration was based on the recognized fact that 
the three-dimensional split-film anemometer gives the most accurate results 
when the probe operates with its axis parallel to the direction of the mean- 
wind vector and the sensors directed into the wind. 

All measurements that are necessary for the calibration are performed at 
the Quality Verification and Calibration Facility at NASA WaUops Fli^t 
Center. The entire procedure is streamlined into a standard form, allowing 
a reasonably fast and efficient calibration of the large number of anemometers 
used in this research program. Based on these calibration procedures, a com- 
puter program was developed, which converts the output voltages of the 
anemometers into velocity components and allows the determination of all 
requited calibration constants. 
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Pig.6. Percent error of the velocity components versus probe*y«w angle for V ■ 6 ft./sec 
(probe No. 1193). 



Fig.7. Percent error of the velocity components versus probe-yaw angle for V > 34.5 ft./sec 
(probe 1193). 


Comparison test 

The three-dimensional split-film anemometer was tested in the atmosphere 
and the results were compared with the results firom two Gill propeller anemo- 
meters. These propeller anemometers were motmted in such a way that one 
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was parallel to the TSI probe axis and the other perpendicular to the first one 
and both in a horizontal plane and located adjacent to the TSI probe. Not 
only the mean velocity components were compared, but also variances, co- 
variances and power spectra measured with the two different types of instru- 
mentation were compared. The results, of these tests are listed in a report by 
Tieleman and Chen [9] , and the results of one set of data are given in Table 2. 

TABLES 


Mean wind componants in tha proba orientad coordinata tyitam 



TSI No. 1192 

Gill anemometers 

U' (ft./sec) 

25.93 

27.91 

V* (ft. /sec) 

6.70 

6.25 

0 (degrees) 

-14.5 

-12.6 

Meant, variances and covariances of the velocity components in the mean-wind coordinate 
system 


TSI No. 1192 

Gill anemometers 

U (ft./sec) 

26.78 

28.60 

V (ft./stc) 

0.0 

0.0 

W (ft./sec) 
(ft./sec)’ 

5.92 

— 

35.91 

26.37 

ulT (ft./sec)* 

2.46 

0.32 

y* (ft./sec)* 

32,93 

22.68 

(ft./sec)* 

25.32 

— 

uw (ft./sec)* 

-4,34 

— 


The velocity components in this table are given in the mean-wind coordinate 
system, which is determined by the direction of the vertical plane containing 
the sample mean velocity vector. The 40-minute data run is broken up in 40- 
second blocks which is necessary for the fast Fourier analysis. The block 
means of the longitudinal velocity component measured by the Gill anemo- 
meter are consistently higher than those measured by the TSI probe; the same 
is true for the sample mean velocity. This is due to the fact that the propeller 
anemometer accelerates faster as a result of a sudden velocity increase and 
decelerates slower as a result of a similar velocity decrease. However, due to 
the limited response characteristics of the Gill anemometers, the turbulence 
quantities obtained from the Gill anemometer are consistently lower than 
those obtained from the TSI anemometer. The same phenomenon is well il- 
fr^strated in Figs. 8 and 9 where the power spectra of the longitudinal and the 
lateral velocity components are compared for the two types of anemometers. 
These results indicate that the Gill anemometer has the proper response capa- 
bilities up to 1 Hz. The spectra measured by t^vc Gill anemometers fall off at 
a faster rate at frequencies higher than 1 Hz while the split-film anemometer 
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spectra show an extensive region where the spectrum function varies as the 
frequency raised to the *~3/3 power. Details of the calculation of the statisti* 
cal quantities and the spectral analysis technique are described in references 
[6] and [9] . 

Conclusions 

A complete data-acquisition and data handling system was designed to take 
high frequency (100 Hz maximum) wind and temperature data measured by 
TSI 1080 three-dimensional split-fflm anemometers. At the present, six probes 
are mounted at 50 foot intervals on the 250-foot meteorological tower located 
at Wallops Island on the Atlantic coast. All probes can be rotated so that their 
axis can be placed closely parallel to the mean wind direction. Hie time dura- 
tion of each data sample is approximately 45 minutes. The output voltages 
from all six probes are simultaneously recorded on analog magnetic tapes. At 
a later point in time, the analog tape is played back and the analog signals 
from each anemometer are digitized, at a rate of 200 samples per second. 

The three-dimensional split-film anemometer system provides seven output 
voltages from which the velocity components in the sensor-oriented coordi- 
nate system and the temperature can be calculated. The general calibration 
procedure as suggested by the manufacturer proved to be complicated and in- 
accurate, as a result, a new complete calibration procedure has been developed. 
This calibration procedure is based on the fact that the anemometer gives the 
most accurate results when the probe operates with its axis parallel to the 
direction of the mean-wind vector and the sensors directed into the wind. The 
data measured with these anemometers are compared with the data measured 
with a set of propellor-type Gill anemometers. Mean values, variances, covari- 
ances and spectra of the two horizontal wind components, measured simultane- 
ously with these two types of anemometer systems, are compared. The discre- 
pancies of the measured quantities can be attributed to the varying and 
limiting response characteristics of the Gill propeller anemometers. 

As a result of carefully carried out calibration procedures, accurate statistical 
estimates of long time series describing the fluctuating wind components can 
be obtained with the TSI 1080 three-dimensional split-film anemometer 
system. 
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A STUDY OF THE MEAN AND TURBULENT FLOW IN THE LOWEST 250-FEET OF AN ON-SHORE 
WIND 


H. W. Tieleman and D. B. Derrington 

Engineering Science and Mechanics Department, Virginia Polytechnic Insti- 
tute and State University 

ABSTRACT 

Three-dimensional velocity measurements and temperature measurements were 
obtained on the 250-foot meteorological tower located near the Atlantic Ocean 
at Wallops Island, Virginia. The wind approached the tower from the ocean, 

As the wind passes the shore line at an oblique angle it experiences a change 
In roughness and possibly a change in surface temperature, resulting in a 
complex three-dimensional type flow. 

Measurements of the mean and turbulent flow quantities were made with a 
special computer-controlled data-acquisition system. Data analysis includes 
the following statistical parameters: mean velocity and temperature, vari- 
ances and covariances including heat fluxes and Reynolds stresses. Nine, one- 
hour runs were analyzed and the results agreed closely with the suggested 
model . 

INTRODUCTION 

The Monin-Obukhov similarity theory is based on the assumptions that the 
flow in the surface layer can be treated as being steady in the mean and takes 
place over a. horizontal and homogeneous surface. The assumptions that 
Coriolis and pressure-gradient forces are negligible and the flow is parallel 
to the surface, lead to a one -dimensional treatment of the problem. Real at- 
mospheric flows, where inhomogeneity of the surface is a common occurrence, do 
not correspond to this ideal situation. Abrupt changes in surface roughness 
and surface temperature lead to horizontal fluxes of momentum and heat. 

An example of this type of flow is the oblique passage of an on-shore 
wind from the ocean to the land surface. Over the ocean a large part of the 
radiation is spent on evaporation of the water and much less for heating the 
water. During a warm surnner day, the warmer air blows over the cold water 
surface and the lower layer of the atmosphere is cooled and an inversion 
develops. For an on-shore wind this stable layer is heated from below after 
it crosses the coastline, and in addition the surface roughness changes from 


smooth to rough, This leads to development of an unstable internal boundary 
layer (IBL) with stable air above it. The depth of the IBL increases with 
increased distance from the shore (Fig. 1). In addition, the overland flow 
becomes three dimensional if th** on-shore wind passes the beach at an ob- 
lique angle. After a certain distance inland from the shore, the IBL may • 
reach some degree of equilibrium when the advection of momentum and heat in 
both the direction of the mean flow as well as in the lateral direction has 
decreased sufficiently in comparison to the vertical heat flux. At that 
point the one-dimensional model could describe the flow adequately, and, 
theoretically, the M-0 similarity theory can then be applied to the lower 
part of the IBL. 

An experimental study is currently underway to measure the mean and turbu- 
lent flow quantities for such flow. Velocity and temperature measurements are 
made with fast- and slow-response type instruments mounted at six levels on a 
250-foot meteorological tower at NASA-Wallops Flight Center. The tower is 
located at Wallops Island, Virginia which is bordered by the Atlantic Ocean on 
the east and a tidal flat on the west. Wind and temperature sensors are 
located at the 30 and the 50-foot levels and from there at 50-foot intervals 
along the height of the tower. The island is one of the barrier islands 
along the coast of the Oelmarva peninsula and consists of a narrow strip of 
dunes 10 feet above sea level and directed appro^^imately along a line from 
North-East to South-West. The distance from the tower to the shore along the 
prevailing summer-wind direction (North-South) is approximately 1000 feet 
(Fig. 2). 

A special computer-controlled data-acquisiti on system is capable of 
sampling the six fast-response anemometer systems at a rate of 200 samples per 
second for a period of about one hour. Data analysis includes the determina- 
tion of the following statistical parameters: mean values, variances, co- 
variances (vertical heat flux and turbulent stress), spectra and cospectra. 
Nine, one-hour runs with reasonable stationary characteristics have been 
analyzed, and the results closely agree with the suggested flow model, 

MODEL FOR THE ON-SHORE FLOW 

The basic governing equations for a thermally-stratified atmospheric 
boundary layer, where temperature does affect the dynamics of the flow, are 
the b&ussinesq equations. In order to deal with the mean flow and the turbu- 
lence separately, Reynolds averaging is introduCied. This process leads to 
the well known Reynolds equations which include the turbulent stresses and the 
turbulent heat fluxes. It is assumed that humidity variations do not affect 
the dynamics of the flow and therefore can be considered as a passive part of 
the air. 

In order to facilitate the study of the flow, let us consider the follow- 
ing simplifications: 

a, Mean and statistical flow parameters do not vary with time. 

b, Mean flow parallel to the surface. 

c, Molecular fluxes negligible compared to corresponding turbulent ones. 
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d. Vertical component of the Coriolis force negligible compared to the 
gravitational force. 

e. Gradients of the pressure field are neglected. 

With these imposed conditions, the equations for the mean flow and mean 
temperature take the following form. 

Momentum equation: 


x-direction (horizontal and in direction of the surface stress) 
y-dr recti on (horizontal) 

dV 

z-direction {vertical) 

0 = - 

Continuity: 

Thermodynamic energy equation: 
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This is the set of equations that governs the flow in the IBL and the over- 
lying inversion. The Coriolis terms can be neglected in the IBL but not al- 
ways in the overlying inversion layer. In the absence of large-scale vertical 
eddies in this layer, no turbulent mixing can occur with the adjacent layers 
and the effect of Coriolis forces can no longer be ignored. 

In the case the flow direction makes an oblique angle with the coast line, 
the flow becomes three-dimensional in the IBL and horizontal momentum and 
heat fluxes exist, although these may be small compared to the vertical 
fluxes. As the flow moves inland the lateral fluxes as well as the longi- 
tudinal fluxes will decrease, and the flow in the IBL gradually recovers its 
one-dimensional character if no other inhomogeneities are encountered. 

The question is whether or not experimental data obtained in this three- 
dimensional flow can be analyzed in the Monin-Obukhov framework, which 
realistically applies only in the constant flux layer. After the IBL has 
developed for a certain minimum distance, it can be assumed that the hori- 
zontal fluxes become relatively unimportant in comparis on t o the vertical 
fluxes. In this case, the local shear velocity, U* = /-uw , the local 
characteristic temperature, T* = - w6/kU*, the buoyancy parameter, g/To, and 
the height z should be the relevant independent variables influencing the 


statistical parameters which describe the turbulence. It is assumed that the 
Coriolis parameter, f, only affects the mean flow and does not affect the 
turbulence appreciably. With this approach one can define a "local” stability 


length, L = 


^ _ 

K g/To we 


, which would convert into the surface stability length. 


Lot as the surface layer is approached. Consequently, the dimensionless 
statistical parameters describing the turbulence should be functions of z/L 
only. In order to conclude whether or not the M-0 similarity hypothesis can 
be applied to the IBL and overlying inversion, the experimental results 
should be compared carefully with the results in the surface layer over uni- 
form terrain. For this purpose the Kansas data will be used as reported by 
Kaimal et al. [2], 


INSTRUMENTATION, DATA ACQUISITION, DATA HANDLING AND DATA ANALYSIS 


The data are obtained with high frequency three-dimensional split-film 
anemometers (TSI-1080 D) and a low-frequency system consisting of cup-vane 
instruments and aspirated platinum-resistance temperature probes. The TSI 
probes are mounted on the south face of the 250-foot tower at heights of 30, 
50, 100, 150, 200 and 250-feet. The cup-vane instruments are mounted on both 
the South-East and the North-West corner of the tower at 50 foot intervals. 
The temperature sensors are mounted on the North side of the tower at 50 foot 
intervals. The data obtained with the low-frequency instruments are recorded 
on digital tape at a rate of 1 sample every two seconds. 


The operation of the three-dimensional anemometers capable of making 
high-frequency velocity and temperature measurements is discussed in detail 
by Tieleman and Tavoularis [6]. Each probe is mounted on a rotor which is 
capable of rotating the probe about a vertical axis so as to align the probes 
approximately into the mean wind direction. A mini -computer controlled data- 
acquisition and data-handling syste" is placed in a trailer located directly 
beneath the 250-foot meteorological tower. 


The output of each anemometer consists of seven voltages, one from the 
thermocouple and six from the hot films. These signals are low-pass filtered 
at 100 Hz, and then sampled and digitized at a rate of 200 samples per 
second. The digital results are analyzed on the IBM-370 computer at VPI and 
SU. A fast Fourier transform algorithm, developed by Tieleman and Chen [5], 
is used to calculate the spectra and cospectra. In view of the computer 
capabilities, the number of data points per block is chosen to be N = 2^^ 
{40.96 seconds in real time). A total of 75 blocks (51.2 minutes in real 
time) were analyzed per run. Because of this treatment of the data, the tur- ' 
bulence can only be analyzed in a frequency range from 0.0244 Hz to 100 Hz. 

The high-frequency limit of the thermocouple is approximately 16 Hz although 
its output is sampled at 200 Hz and analyzed up to 100 Hz. Consequently, 
variances and covariances as well as spectra and cospectra involving the 
temperature have to be analyzed differently in order to account for the noise 
contribution in the frequency range between 16 Hz and 100 Hz. 


Four computer programs for the analysis of the data have been developed, 
rhe first prog^’am converts the voltages to velocity components and tempera- 
ture, and in addition records those quantities on digital tape. These 
velocity components are obtained in the TSI-sensor oriented coordinate sys- 
tem l 5»6]. The second program calculates block n«ans and block variances of 
the velocity components and temperature. In addition, the probe-yaw angle, 

6, (horizontal angle between mean wind and probe axis) aiid the total number 
of reverse arrangements [1] of the means and variances, are calculated at 
this point. Inspection of the block means and block variances and their 
total number of reverse arrangements is the test for stationarity of the run. 
Based on this information a decision is made whether or not to continue the 
analysis. The third program converts the velocity components in the mean- 
wind coordinate system (Fig. 2), calculates the block means and removes the 
block means from the data in each block. Next, all the possible variances 
and covariances are calculated for each block and averaged to obtain the mean 
variances and covariances for the entire run. Also the velocity and tempera- 
ture data with the mean removed are recorded digitally once more. With the 
fluctuating velocity components available in the mean-wind coordinate system 
and with the fluctuating temperature, the spectra and cospectra of the fol- 
lowing pairs of variables, uw, wo and wv are obtained. The variance of the 
temperature and we covariance were obtained by integration of the e spectra 
and we cospectra in the frequency range from 0.0244 Hz to 16 Hz. 

DISCUSSION OF RESULTS 

Until recently, data were only collected during the summer months. Since 
the prevailing winds are from the south during the summer, the split-film 
probes were mounted on the south side of the tower. Originally, thirteen one- 
hour data runs were taken; however, four of these runs were rejected as being 
nonstationary. Table 1 shows the results of the mean flow and the turbulence 
of the nine runs which were analyzed. The values of the direction angle, 
indicate the direction of the mean-wind vector measured clockwise with respect 
to north. 

The results show clearly the presence of an unstable inner layer and a 
stable outer layer. For each run, the vertical heat-flux changes direction 
between 50 and 100 feet. The results (not all listed in Table 1) indicate 
that the flow is three-dimensional, fluxes of heat and momentum in the hori- 
zontal direction are appreciable, but are usually less than the vertical 
f 1 uxes . 

For most of the runs, the velocity profiles when plotted in the usual 
semi -logarithmic fashion show one or more "kinks." Extrapolation of the 
velocity from the 50 and 30 foot elevations lead to extremely high roughness 
lengths, Zq, especially for runs 2, 7, 8 and 9 (zq > 1 foot). These values 
of Zq are typical for flow over moderately rough terrain, such as wooded 
areas and suburban areas, but are too high for the terrain at Wallops Island 
in this case. Consequently, these extremely high roughness lengths indicate 
the presence of an abnormal boundary layer. The velocities at the higher 
elevations are much higher than under normal conditions, when large-scale 


turbulence provides, as a result of turbulent mixing, a much more uniform 
velocity profile. In the absence of large-scale eddies, Coriolis forces may 
become sufficiently important to affect the flow. The results do indicate 
this, runs 1, 7 and 8 show a maximum clockwise rotation of the mean flow 
direction of 4 to 10 degrees between the 50 and 250- foot levels. 

The mean-temperature profile for the nine runs show some unusual "kinks" 
mostly between the 100 and 150 foot levels. The most interesting sequence of 
mean- temperature distributions is from runs 7, 8 and 9 (Fig. 3), which were 
taken on the afternoon of August 13, 1976 starting at 2:12 P.M., 3:37 P.M. 
and 5:05 P.M., respectively. The early afternoon runs show a decrease in 
temperature up to the 100-foot level and an increase in temperature for 
higher elevations. In the early afternoon the layer above the 100-foot level 
shows a rather strong positive temperature gradient. As the afternoon pro- 
gresses this temperature gradient decreases dramatically. The run starting 
at 5:05 P.M. shows a decrease in surface temperature with a subsequent cool- 
ing trend in the lower 50 feet. 

Many of the afternoon runs show a decrease in temperature with height up 
to the lOO-foot level, then a slight increase between the 100 and 150 foot 
levels and decreasing temperatures above the 150-foot level. The temperature 
profiles obtained outside the IBL are governed by the depth of the inversion 
layer over the ocean. Apparently, this depth was larger than 250 feet dur- 
ing the early afternoon (Run 7); the second run of the afternoon (Run 8) 
still shows a positive temperature gradient but of a much smaller magnitude. 
For the third run (Run 9) the temperature gradient between 100 and 150 feet 
is slightly positive but above the 150-foot level it is negative. This indi- 
cates that during the afternoon the inversion layer is decreasing in depth 
but also that its degree of stability is decreasing. 


No simple relationship exists between the standard deviations of the 
three velocity components normalized with the local shear velocity and the 
stability parameter z/L for extreme stable conditions, z/L > + 2.0 (Figs. 

4, 5 and 6). For near neutral conditions the average values of oj^* are 
1.75, 1.5 and 1.4 for a = u, v and w, respectively. The values of the vari- 
ances in stable and near-neutral air depend on the frequency response of the 
instruments with which the velocities are measured and the manner in which 
the data are analyzed. Low-frequency contributions affect the horizontal 
variances significantly, and consequently the way in which the means are re- 
moved from the original time series can influence the magnitude of these 
variances a great deal. The results also indicate that the ratios cr^/U* are 
very sensitive with stability, and increase drastically for deviations from 
neutral conditions in both the stable and unstable regime. The standard 


deviations of the velocity fluctuations, 0 „, show maximum values in the IBL 
and decrease rapidly across the interface between this layer and the overlying 


inversion. In the stable air the vertical variances are usually smaller than 
the horizontal variances. Inspection of Table 1 reveals that the vertical 
fluctuations almost disappear for those runs with clockwise rotation of the 
mean wind with height. This indicates that turbulent mixing in this layer 
with adjacent air does not exist and Coriolis forces become important. 


The ratio, ae/lT#!, shows a great deal of scatter close to neutral thermal 
Stratification as T* should approach zero (Fig. 7). Values of 2.0 or larger 
are common in the IBL. In the stable air above this layer, values of this 
ratio decrease sharply with increased stability to an average of 0.5. 

Panofsky [4] reports that aa/lT*! shows an average of 1.8 in stable air and 
decreases with increasing stability similar to the results obtained at Wallops 
Island. In this study the temperature ratio, oh/jT*!, was obtained with the 
local scaling temperature T* rather than the value of this quantity in the 
surface layer. 

The values of the uw covariances are all negative as expected, except for 
run 8, 150-foot level. The local shear velocities obtained from these co- 
variances are successfully used as reference velocities for the calculation 
of the scaling temperature, T*, the stability parameter, z/L, and the flux 
Richardson number. The uw covariances show maximum values at the 30-foot 
level except for run 8 and 9. Extremely small values for uw were obtained in 
the stable air above the IBL. The lowest value was obtained at the 250-foot 
level for run 7, for which the largest rotation of the mean wind direction 
was measured. The results in Table 1 indicate that extremely small values 
of uw not automatically mean absence of large scale eddies. Runs 1 and 8 
which show clockwise rotation of the mean-wind direction with height, also 
show larger uw covariances than those from other runs which do not experience 
any appreciable rotation of the mean-wind vector. These results seem to in- 
dicate that the contribution to the shear for runs 1 and 8 must come from 
the velocity fluctuations with smaller wave lengths. 


The we covariances show positive heat flux in the IBL and negative values 
in the overlying inversion layer. The values of these covariances are of the 
same magnitude, indicating that the upward heat flux in the IBL is of the 
same order as the downward heat flux in the stable layer. Because of the 
absence of large scale eddies in runs 1, 7 and 8 it can be concluded that 
heat is more effectively transported by the smaller eddies. Considerable 
horizontal heat transfer occurs at the interface between the IBL and the 
overlying stable layer. 


The characteristic length-scale of the energy-containing eddies can be 
obtained from the normalized autocorrelation function using Taylor's 
hypothesis. 



R(x)dT 


( 6 ) 


In Table 1 the values of Lj^ were obtained from the measurements made with the 
NASA cup-vane instruments. These instruments do not respond to high-frequency 
velocity fluctuations and consequently, the calculated values for the integral 
length-scale, Lj^ may be somewhat high. 

An alternate method for obtaining a measure of the scale of the energy 
containing eddies is the wave length corresponding to the reduced peak fre- 
quency of the logarithmic longitudinal velocity spectra: 


( 7 ) 


Kaimal [3] gives the following relation for the two length scales 

xj; = 2. (8) 

According to the results listed in Table 1, there does not seem to exist a 
simple relation between and . For most runs, X™ > for the unstable 
flow in the IBL. This trend is reversed for the stable flow above the IBL; 
the integral length scale, Ljj , is in most cases larger than Ajt* by a factor 
of 1.5 - 2.0. The length scale, , is smaller in the stable air for those 
runs which show rotation of the mean wind, as is to be expected. The block 
length used for the analysis of the cup-vane data is about 512 seconds, and 
removal of the block means results in high-pass filtering at a lower fre- 
quency. Consequently, large scale fluctuations are included in the determi- 
nation of the autocorrelation, resulting in large values of the turbv?ience 
integral scale. These low frequency fluctuations which can be identified as 
gravity waves, are present in run 4 at all levels above the IBL and at the 
250-foot level of run 7. 

CONCLUSIONS 

The results of this experimental study indicate that the stable on-shore 
flow, modified by changes in surface roughness and surface temperature, does 
not lead to an ordinary boundary layer with the usual verticaT turbulent mix- 
ing. The measurements support the physical model which include the develop- 
ment of an unstable internal boundary layer (IBL) with a stable overlying in- 
version. The magnitudes of the variances of the velocity components and the 
uw covariances decrease rapidly across the interface between the two layers. 
The vertical heat flux is positive in the IBL and negative in the overlying 
inversion. The roughness lengths, Zq, are unusually high, indicating that 
the velocity profile is not the normal profile one would encounter in an at- 
mospheric boundary layer with large-scale turbulent mixing. Under certain 
conditions these large-scale eddies are absent in the inversion layer and the 
Coriolis forces gain importance resulting in a noticeable clockwise rotation 
of the wind direction above the IBL. 

The air in the inversion layer and in the IBL are essentially independent, 
and ordinary boundary-layer models based on turbulent mixing cannot be ap- 
plied. Consequently, in absence of turbulent momentum exchange with adjacent 
layers, the wind velocity of the stable layer is larger than the velocity in 
the adjacent layers. This explains the existence of a low level jet with a 
maximum velocity at an elevation of approximately 300 feet. It is believed 
that the presence of the low-level jet depends on the depth and the degree of 
stability of the inversion layer above the ocean just before it is modified 
by travelling over land. 

The presence of the low-level jet with the corresponding clockwise rota- 
tion of the mean wind vector, depends on the existence of large-scale eddies 
which promote the turbulent mixing between the layers. Large positive values 
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of the stability parameter, z/L and flux-Richardson number are necessary but 
not sufficient for the absence of these large-scale eddies. At the present 
time it is not clear under what conditions this type of flow will occur. In 
order to shed some light on these questions it will be necessary to study the 
stable boundary layer over the ocean before it is being modified by the land 
surface. 
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Figure (1) - Model for the temperature profiles for on-shore flow and the inter 
nal boundary layer. 
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Figure (5) - Ratio of standard deviation of lateral velocity and local 
friction velocity as a function of z/L, 



Figure (6) - Ratio of standard deviation of vertical velocity and local 
friction velocity as a function of z/L. 



Figure (7) - Ratio of standard deviation of the temperature and local scaling 
temperature as a function of z/L. 
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Atmospheric Turbulence Below 75m in the Convective 
Boundary Layer (Strong Wind Conditions) 

by 

H. U. Tieieman and S. E. Mullins 
Virginia Polytechnic Institute and State University 
Department of Engineering Science and Mechanics 
Blacksburg, Virginia 24061 

Summary 

Mean as well as turbulent velocity and temperature measurements obtained from 
the 76m micrometeorological tower at Wallops Island, Virginia are discussed. The 
measurements were obtained with cup*vane instruments and aspirated temperature 
probes (5 levels) for mean profiles as well as with hot-film probes (6 levels) for 
turbulence parameters. The wind conditions during data acquisition can be classi- 
fied as moderately strong (at least 10 m/sec at the 15m level) with winds coming 
from directions varying between north and west. These conditions are usually en- 
countered following the passage of a cold front. Observations have been compared 
with strong wind results as listed in several review papers. A new scheme is 
developed for the normalization of the velocity spectra using a reduced frequency 
related to the characteristic length scale of the velocity components associated 
with the inertial subrange. 

Introduct ion 

In the field of wind-engineering dealing with wind loads on structures, it is 
generally assumed that under strong wind conditions the flow in the atmospheric 
boundary layer is controlled by the mechanically produced turbulence [1,2,3]. In 
general, the turbulence in the planetary boundary layer (PBL) is either mechani- 
cally produced or buoyancy produced. Recent observations throughout the PBL have 
indicated that significant entrainment of heat and momentum from the free atmos- 
phere is still another source of production of turbulent energy, specifically in 
the low-frequency range [4,5,6]. Depending on the atmospheric conditions and 
surface conditions any combination of these three turbulence-producing mechanisms 
can play an important role in the development of the mean flow and the turbulent 

structure in the PBL even under moderately strong wind conditions. 
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Based on observations [4,5] and theoretical considerations [7]> the convective 
boundary layer over un i form terrain has been divided into three different layers: 

1. The surface layer (lowest 10% of the PBL or less) where the wind shear 
plays a dominant role and the momentum and the heat fluxes are reasonably constant 
with height. Monin-Obukhov similarity applies in this layer to the statistical 
parameters involving the vertical velocity component, but not to those of the 
horizontal velocity components. Panofsky [8] has suggested that in the surface 
layer, proper scaling of the standard deviation, a, of the horizontal turbulence 
components is of the following form 

o/U* - f(z./L) (I) 

where U* is the surface friction- veloci ty, z. is the height of the PBL del ineated 
by the location of the lowest inversion layer and L the Mon i n^Obukhov length. 

2. Above the surface layer a free convection layer exists where the con- 
vective turbulence dominates the mechanical turbulence (-z/L>l). 

3. The remaining part of the PBL is called the mixed layer, where, as a 
result of strong mixing, the mean flow as well as the turbulence are no longer 
functions of height and the surface stress. 

Overlying this three-layered structure is an inversion layer through which consid* 
erable amounts of heat and momentum enters the PBL from the free atmosphere. 

Observations made over extremely flat and uniform terrain in northwestern 
Minnesota [4] show a net downward heat flux in the upper parts of the mixed layer. 

It was also observed [5] that the entrained warm air descends in the form of dis- 
crete plumes and upward heat flux from the surface takes place in the form of dis- 
tinct thermals In an alternating fashion. Moreover, the spectra of the horizontal 
velocity components in the surface layer follow two different similarity schemes [6]. 
The inertial subrange of these spectra follows Monin-Obukhov similarity, whereas 
the low-frequency range follows mixed-layer similarity with z. as the controlling 
parameter. 

How the flow in the convective PBL will be modified as a result of nonunifor- 
mities in upstream terrain is difficult to predict. The terrain variations have 

a more pronounced effect on the flow as it passes several times from land to water 
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and back to land. In such a case the air will be In a transient state, being 
affected by variations in the surface fluxes of moisture, heat and momentum. The 
northwestern flows observed at Wallops Island will have crossed in succession the 
U.S. mainland, the Chesapeake Bay (20-50 km), the "Oelmarva” Peninsula (20-50 km) 
and for the last 3~5 km the tidal swamp separating Wallops Island from the penin- 
sula. The last 60-100 m, the air travels over the island before it arrives at 
the tower. During the daytime, the upward heat flux from dryland will be much 
higher than from water or moist swamp terrain. Consequently, the upstream flow 
will experience a discontinuity in surface heat flux as it crosses the tidal swamp, 
and the downward heat flux, which is the result of entrainment of warm air from 
the free atmosphere, can be expected to penetrate further downward into the PBL. 

This paper deals with the discussion of mean flow and turbulence measurements 
obtained from the 76-meter micrometeorological tower at Wallops Island, Virginia. 

The data were obtained under strong wind conditions of at least 10 m/sec at the 
I5"m level and from directions varying between north and west. These conditions 
are usually encountered following the passage of a cold front. In addition to 
daytime data, some nighttime observations are reported for comparison. One set 
of observations was obtained from a direction where a 6.5"m high rocket-fuel 
storage bunker is approximately 85 m upstream of the tower The strong wind ob- 
servations have been cpmpartsd with the results from the review papers ESDU, 197^ [2] 
and Couni han [3]. 

A new scheme is developed for the normalization of the velocity spectra in 
the high-frequency range using a reduced frequency related to the characteristic 
length scale of the turbulence associated with the high-frequency range. More 
stable estimates of the characteristic length scale of the energy-containing eddies 
can be obtained in this manner as compared to the traditional integral scales ob- 
tained from the autocorrelation function. 

Experimental Details 

The Wallops Island tower is instrumented at six levels with hot-film anemometer 
probes capable of high-frequency velocity and temperature measurements. A descrip- 
tion of these probes, data acquisition and data hand 1 ing system as well as a 


description of the data analysis can be found in detail in a report by Tleleman and 
Tavoularis [9] and in a shorter version in a paper by the same authors [10]. In 
addition, cup-vane velocity instruments and aspirated temperature sensing probes 
used for the measurement of mean qu-antitites are mounted at five levels. The data 
obtained with the hot-film probes were sampled 200 times per second; the output of 
the wind- and temperature-profile sensors once every two seconds. A total of six 
55-minute periods of observation with the hot-film instruments as well as approxi- 
mately 25 runs of varying duration (between 1*8 and 80 minutes) with the cup-vane 
and aspirated temperature sensors form the data base for this paper. The data from 
the hot-film probes are analyzed In two frequency ranges, the high-frequency analy- 
sis in the range of 0.02l*l*-l00 Hz and the low-frequency analysis in the range 
0.00153-6.25 Hz. The data from the cup-vane and aspirated temperature sensors are 
analyzed in the range of 0.00195-0.25 Hz. Computation of spectra followed the 
procedures as described in a report by Tleleman and Chen [II]. 

Experimental Results 
a. Mean Wind Profile 

Mean-wind profiles under strong wind conditions over relatively smooth and 
uniform terrain follow closely either the logarithmic law 

U/U* » l/k In z/z^, (I) 

o o 

or the power law 

U/U^ * (z/z^)® 

where U* Is the profile friction velocity, k the von Karman constant (k = 0.4), 
z^ the roughness length, Uj is the mean velocity at some reference height z^ 
and a the power index. 

The height of the layer where these profile laws apply depend on the 
wind speed, the terrain roughness as well as the surface heat flux and 
the atmospheric conditions at the upper edge of the PBL. The Wallops velocity 
profiles show that this height, even under relatively strong wind condi t ions, may 
be limited to 50 m or less (Tables I and tl). Mean temperature and turbulent 

heat flux observations show a countergradient downward heat flux to within 


a height of 9.1 m above the ground for a.M runs listed in Table M, including the 
evening run. These results indicate that warm air is transported downward to the 
surface, possibly as a result of entrainment of warm air Into the PBL from the 
free atmosphere. Low-frequency components in the measured velocity spectra, for 
all runs listed in Table II with the exception of run 16, similar to those observed 
by Kaimal [6], seem to point to the existence of a mixed layer. Such flow will 
require a transition from the velocity profile in the surface layer to a uniform 
velocity profile In the mixed layer showing up as a kink in the profile. Changes 
of this type in the mean velocity profile have been observed for several runs at 
a height of approximately 50 m. These profile changes are definitely not associ- 
ated with a change in upstream roughness as is observed for run 186 A with the 
bunker about 85 m upstream from the tower (Figs. 1 and 2). The variation of the 
power index with the corresponding roughness length for all Wallops data tend to 
fal'l between the relations suggested by ESDU 1972 (12] and Counihan [3] (Fig. 3). 
This Indicates that the Wallops wind profiles below the kink fall in the strong 
wind category, 
b. Turbulence ParameterE 

1. Turbulent Stress 

The measured local shear velocity defined as 

U* « [(uw)^ + (3) 

obtained from the hot-film probes in the lowest 30 m and analyzed in the high-fre- 
quency range (0.02A4j<n£l0O Hz) compared quite well with the profile friction 
velocity U*. This indicates that the velocity profile up to 30 m is matched with 
the turbulent stress obtained in the high-frequency range, implying that the 
foreign low-frequency components do not affect the velocity profile appreciably 
in this layer and that the mechanical turbulence dominates. Moreover, the results 
indicate that the profile friction velocity is an accurate measure of the 
local shear velocity in this layer. 

2. Standard Deviations and Turbulence Intensities 

Standard deviations of the velocity components were obtained from the cup- 
vane instruments (Oy, o^) in the frequency range of 0.00195"0*25 Hz, and from the 


hot-film probes (o , o and a ) in tho frequency range of 0.00153*6.25 Hz. Be- 
cause of the limited frequency response, from the cup-vane system had to be 
adjusted by a factor of 1.17 after being compared with from the hot-film instru- 
ments. The results clearly indicate that the ratio at the 15 m level is 

consistently lower for those runs with the convective kink In the velocity profile 
(table I and It). With a kink the average - 2.5, and without a kink the 

average - 2.9, both measured at the 15 m level, and decreasing to 2.2 and 2.7 

respectively at the 76-m level. The variations in a /\J* and o /U* do not seem to 

V O wo 

follow any systematical pattern. Average values of vary from 2.5 at the 

15*m level down to 2.3 at the 76-m level, while the average value of “ 1«2 

is uniform over the depth of the layer investigated. 

Mixed layer similarity as suggested by Panofsky [8] does not work for the 
Wallops data, probably as a result of slightly nori- homogeneous upstream terrain 
but more importantly as a result of the variation in upstream surface heat flux. 
Consequently, for those involved in practical situations when perfectly uniform 
conditions are rarely encountered, rt seems more advantageous to use turbulence 
intensities. Figures 4, 5 and 6 show the averaged turbulence intensities of all 
three components compared with the intensities suggested by fSOU, 1974. The com- 
parison is excellent except for the lateral turbulence intensities which are con- 
sistently higher than those predicted by ESDU , 1974 [2]. 

3. Turbulence Scales and Spectra 

The characteristic length scales of the energy-containing eddies of the tur- 
bulence In the atmospheric surface flow are used quite extensively, especially for 
geometrical ly scaled model studies of the flow In wind tunnels. This turbulence 
length scale is also important in relation to the characteristic size of a structure 
exposed to a corresponding turbulent flow. Oetermlnation of these scales can be 
achieved by several methods: 

a. integration of the autocorrelation coefficient with respect to the time lag up 
to the first zero crossing. The autocorrelation coefficient can be obtained either 
directly from the velocity data or via the Fourier transform of the spectrum 


function. 


b. Determination of the wavelength corresponding to the peak of the logarithmic 
power spectrum. 

c. From the best fit of the measured spectrum function to the corresponding 
empirical von Karman spectrum functions. 

d. By estimation of the velocity spectrum function at zero frequency 


e. 


2 

L » U/4 I im (S (n)/a ) for a ■ u, v and w 
“ . n -^0 “ “ 

From theoretical considerations of kinetic energy and dissipation rate. 

L ■ <yjz a « u, V and w 

QL a 


ik) 
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(5) 


Ideally, these methods should give nearly identical results for a given data 
set. Most of these methods however have serious shortcomings especially when 
low-frequency components are present in the data sample. The first method is 
reasonably accurate for shear generated wind-tunnel turbulence but gives very un- 
stable estimates of the integral length scales of atmospheric turbulence es- 
pecially in the presence of low-frequency components. The integral length scaies 
do vary considerably with the method the autocorrelation coefficient is obtained. 
Results show that the scales obtained with the direct method are considerably 
smaller than those obtained via the spectrum function (Figs. 8, lO). Moreover, 
length scales obtained from the cup-vane system are on the average smaller than 
those obtained from the hot-film instruments (Fig. 8). In the case Iqw-frequency 
components are present, the logarithmic power spectrum does not have a well defined 
spectral peak and the second method gives unreliable scales. With low-frequency 
components present, the von Karman spectrum function can only be fitted adequately 
to the high frequency part of the measured spectra (Figs. 11-13). This requires 
that the velocity spectra need to be accurately measured in this range up to a 
frequency of at least 60 Hz for moderately strong wind conditions. Consequently, 
under these conditions this method can not be used with data obtained from cup 
or propeller anemometers. Again when low-frequency components are present, esti- 
mation of the spectral function at zero frequency gives rise to a great deal of 
uncertainty making the fourth method useless. 


The fifth method requires knowledge of the dissipation rate e, which can be 
obtained from the Kolmogorov's spectral law for the inertial subrange 

2/T ‘5/3 

S (k,) ■ K k, a » u, V and w (6) 

ala I 

where k^ is the wavenumber in the direction of the mean flow and is usually ex- 
pressed in terms of frequency by Taylor's hypothesis as kj ■ 2ffn/U. are the 
Kolmogorov's constants, K » K ■ ^»/3 K assuming isotropy. The universal con- 
slant, K^, has been established to be 0.5* In terms of the Euler ian frequency- 
domain the Kolmogorov's law becomes 

nS„(n) - k (7) 

a a 

where A.. ■ 0.147 and A * A * 0.196. Whether or not this relation can be applied 
U V w 

to the turbulence measured at Wallops Island remains to be seen. The dissipation 
rate is an impossible quantity to measure independently, without making use of the 
local- isotropy assumption. 

Assuming that the frequency version of the Kolmogorov's law applies, values 

of the dissipation rate e for all three velocity components have been calculated. 

The results indicate that the dissipation rates, obtained from each individual 

velocity spectrum for the same data sample, are not identical. The isotropy 

prediction of 4/3 for the ratios of either one of the lateral spectral values and 

the corresponding streamwise spectral value is clearly not satisfied in the range 

-2/3 

where the measured logarithmic spectra functions closely vary as n . These 
results indicate that the existence of isotropy and the validity of Taylor's 
hypothesis may not apply to the turbulence encountered at Wallops Island. Until 
these questions are answered adequately, evaluation of the dissipation rate in the 
inertial subrange Is unreliable, and consequently cannot be used to obtain in- 
tegral length scales of the turbulence accurately. 

The evaluations of these length scales from different instruments and by 
different methods clearly indicate that in the presence of low-frequency components, 
the predictions of the integral scale from ESOU [2] are much too low for the moder- 
ately Strong wind conditions encountered at Wallops Island (Figs. 8-10). 


Consequently I an alternate method has to be found to obtain better estimates 
for the Integral length scales of the atmospheric turbulence. Therefore It Is pro- 
posed that a scheme similar to the one used by Kalmal [14] be adopted In the field 
of wind engineering. When logarithmic velocity spectra normalized with the total 
variance are plotted against a modified frequency, spectra for like ve- 

locity components collapse In the high-frequency range (Figs. 1 1-13) « Here, f Is 
the reduced frequency nz/U and f^ Is the reduced frequency at the Intersection of 
the line of best fit In the range where the logarithmic spectrum Is proportional 
to with the line n S (n)/o^ « I. For the v and w spectra, obtained from 

data without low-frequency components, a clearly discernible logarithmic spectral 
peak Is located at [15]. If f^ denotes the peak reduced frequency then 

f** = 3-4 f** for a » V, w (8) 

m o 

corresponding to a value of 3«53 obtained from the v, w - von Karman spectrum 

function. In general no clearly discernible spectral peak can be obtained from 

the measured u spectra and therefore it is proposed to use a value for the ratio 

f /f obtained from the von Karman spectrum function for the u component 
fn o 

f“ = 3.85 f‘' (9) 

m o 

In the case low-frequency components are present, and no clear spectral peak 

exists a value of f can be obtained from the measured spectra and a corresponding 

o 

value for f can be calculated with either (8) or (9). In this case f corresponds 
m m 

to a ficticious peak reduced- frequency which would have existed in the case the 
turbulence was purely shear- gene rated and low-frequency components would have 
been absent. A related peak wave length can then be obtained from 

X * U/n_ * z/f OO) 

m ni m 

This peak wavelength provides a length scale characteristic of the turbulence in 
the frequency range corresponding to the range in which the natural frequencies of 
most structures such as buildings and bridges fall. It is the turbulence in this 
range, with wavelengthscorresponding to the characteristic length of most structures 
which Is responsible for the dynamic response of structures and Is also responsible 
for the positive and negative local extreme pressures present on the exterior 


surface of structures. From the point of view of the engineer dealing with the 
modelling of the atmospheric surface flow in the wind tunnel in order to conduct 
wind-load model studies, it is essential that the atmospheric turbulence in this 
frequency range is modelled accurately. The lotJing effect of the low-frequency 
velocity components, corresponding to wavelengths larger than can be considered 
as being the result of a non-steady mean flow. Consequently, it seems that this 
new method provides meaningful integral length scales corresponding to that part 
of the atmospheric turbulence which has a pronounced effect on the dynamic be- 
havior of structures. 

Figures il-l3 show the logarithmic spectra of the three turbulence components 
for run I9» plotted against the modified reduced frequency f/i^^ with the corres- 
ponding modified von Karman spectrum functions superimposed. Collapse of the 
spectra in the high-frequency range is of course satisfied, and the presence of 
low-frequency components in all three spectra is evident. The v spectra show a 
somewhat different shape between the high-and low-frequency range quite similar to 
the V spectra obtained at Minnesota [6]. Kaimal {6] explains this odd spectral 
shape to be a transit! ion between the high-frequency range (inertial subrange) 
corresponding to the w-spectrum and the low-frequency range corresponding to the 
u spectrum. 

The distributions of averaged values of with height for the three velocity 
components are shown in figure 14. Because of a great deal of variation in at- 
mospheric conditions and of non-uniformity in upstream surface conditions, re- 
sulting in appreciable scatter of the data, it is clear that under these circum- 
stances no simple normalization scheme is available. The values of X obtained 

(n 

over extremely uniform terrain [5] and under uniform convective conditions show a 
systematic behavior when plotted in dimensionless form using the height of the 
boundary layer, z., as characteristic length scale. 

The advantage of the use of peak wavelength over the traditional integral 
length scale is that uncertainties due to sample rate variations and method of 
analysis are being eliminated. If the spectral functions and the standard de- 
viations (including the low“frequency components) of the three velocity components 


can be determined accurately only then a reliable estimate of the peak wavelength 
can be obtained. That these peak wave lengths vary a great deal with atmospheric 
conditions, time of day and variations in upstream surface conditions is a matter 
that needs further investigation. 

Conclusions 

The results from this research clearly indicate that under moderately strong 
wind conditions from directions varying between north and west, the turbulence in 
the surface flow below 76 m at Wallops Island is not all of mechanical origin, in 
addition to the shear-generated turbulence, there exist low-frequency velocity 
fluctuations apparently due to the entrainment of turbulent momentum into the PBL 
across the lowest inversion. 

The flow is certainly not neutrally stratified although bul k- Richard son numbers 
are relatively small. However, as a result of a stretch of 3"5 km wet tidal 
swamp just upstream of the tower location, the vertical surface heat flux is being 
interrupted and warm air which has entered the PBL from the free atmosphere is 
allowed to penetrate downward against the temperature gradient to within 9 m above 
the surface. 

Consequently, the mean flow and the turbulent flow in this layer cannot be 
normalized with the usual surface-layer parameters. The flow encountered below 
76 m at Wallops Island compares in many aspects quite well with the flow measured 
at the uniform Minnesota site t^» 5, 6J. However, as the result of the north west 
wind passing twice from land to water and back to land before reaching the Wallops 
island tower, the vertical heat flux is being Interrupted giving rise to modifi- 
cations of the flow. No simple similarity scheme is available at the present time 
to describe this type of flow, even if the wind conditions can be classified as 
moderately strong. 

The following general conclusions can be made for the daytime flow below 76 m 
over relatively open terrain with extended areas of water upstream, under moderately 
strong wind conditions occurring after the passage of a cold front: 

1. The height of application of the log-law is limited, transition from this 
velocity profile to the uniform profile in the mixed layer starts occasionally 


below 76 m. 


2. The air In this layer is not neutrally stratified. As a result of reduced 
vertical convective heat flux, warm air from the free atmosphere is transported 
downward against the temperature gradient. 

3. The turbulence in this layer is not all of mechanical origin. Low frequency 
velocity fluctuations in all three components are present as a result of entrain- 
ment of turbulent momentum from the free atmosphere across the lowest inversion 
layer. 

k. Under these conditions the turbulence intensities except for the lateral 
component can be estimated quite well from ESOU, 1974 [2]. Measured values of 
Oy/U are consistently higher than those predicted by ESDU, 1974 [2]. 

5* In the presence of low-frequency fluctuations, traditional Integral length 
scales of the turbulence are much larger than those predicted by ESOU, 1974 [2]. 

6. Measured integral scales vary a great deal with instrumentation and type 
of analysis. 

7. Instead of the traditional integral length scale, it Is proposed that the 

peak wavelength, A , be used as a measure of the characteristic length scale of the 
m 

energy-containing eddies. The peak wavelength can be obtained from the reduced 

-2/3 

frequency, f^, at the intersection of a line of best fit through the n -region 

2 

of the logarithmic spectrum and the line S (n)/a * i. 

a OL 

8. With low-frequency fluctuations present the von Karman spectrum function 
does not fit the spectral data in the low-frequency range. 

S. In the case an independent method is used for the determination of the 
integral length scale, the von Karman spectrum function does not fit the spectral 
data In the high-frequency range either. 

10. Velocity spectra normalized with the reduced frequency, f , collapse to a 

-2/3 

single curve in the n region. The von Karman spectrum functions, modified In the 
same manner, fit the data well in the high-frequency range. 

11. The measured spectra of the v-component show a transition between the high 
and the low frequency ranges similar to the spectra obtained at Minnesota [6]. This 


transition is not predicted by the von Karman spectrum function for this component. 

12. In this layer the peak wavelengths, obtained from the reduced frequency 
f^, vary linearly with height. 
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Fig. 1 Typical strong wind velocity profiles, logarithmic representation. 
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Fig. 2 Typical strong wind velocity profiles, power law representation. 
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Fig. 3 Variation of power law exponent with roughness length. 
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Fig. 4 Averaged turbulence Intensity, steamwise component. 
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Fig. 6 Averaged turbulence intensity, vertical component 








PAGE JS 
Of POOR QUAt»Tf 





Jijfk kk 'Ifk«, ■,,.■• 
S' > 1 \, I 1'J» ‘ 



-41- 


4|- 


'( -'ll! I % 

RECORD OF STREAMWISE VELOCITY COMPONENT, 
RUN 19, 250 ftrLEVEL, U»I3.9 m/sec , (Ty =1.60 m/sec, 
3-25-'77, 14:21:00, lJ *306 m 

RECORD OF LATERAL VELOCITY COMPONENT, 
RUN 19, 250 ft.-LEVEL, Oy =l.65m/sec, Lj»479m 







-4I- RECORD OF STREAMWISE VELOCITY COMPONENT, 
RUN l6,250ft.-LEVEL,U=ll.4m/sec, 0 -^= 0.81 m/sec , 
3“23“'77, 20:30:00, L5=77.9m 



RECORD OF LATERAL VELOCITY COMPONENT, 
RUN 16, 250 ft - LEVEL, (Ty =0.62 m/sec, Lj=23.3m 


Fig. 7 Records of velocity components. 
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ABSTRACT 


VcXocity and temperature moasuromente wore obtained from an instrumented 76-meter 
meteorological tower located near the Atlantic Ocoatv at Wallops Island, Virginia, 
The instrumentation consists of a slow-rcsponse cup-vane and resistance temper- 
ature (profile) system and a hot-film, thermocouple >systcm for turbulence measure- 
ments. Results are presented for moderately strong winds from westerly directions 
(category I) and for on-shore winds from southerly directions (category IX), 

Uesults from category I winds indicate the presence of low-frequency velocity 
fluctuations affecting all turbulence parameters similar to the observations made 
at Kansas and Minnesota. 

Winds of category II eKpericivce a change in surface roughness and surface temper- 
ature as they cross the coast line, resulting in the development of an internal 
boundary layer (IBL). The stable ocean air above the IBL shows extremely low 
turbulence levels, often loss than 3%, Because of th:i lack of turbulent mixing 
with adjacent layers, Coriolis effects are important and velocity profiles with 
a maximum velocity (low-level jet) at heights between 60-200 m exist. 

introduction 

Structures and structural elements arc usually deslguod to support the maximum 
loads tl\ac are likely to bo imposed upon them within the lifetime of the structure. 
Under the strongest wind conditions the assximption of a neutrally stable plane- 
tary boundary layer (PBL) Is probably satisfied. However, problems of fatigue 
and serviceability of a structure arc often overlooked in its design. These 
problems can become accute under much less severe or xmder moderate wind condi- 
tions, Consequently, in the field of wind loading on structures and their re- 
sponse, one needs Co study the flow in the HBb not only under extreme wind con- 
ditions but also under more conunon moderate conditions. Flow characteristics of 
the PBb under chesQ conditions, mostly corresponding to flow over relatively 
homogeneous terrain have been svmmmrized in several review papers (1,2,3, 4). 

These results have formed the basis for empirical models of the wind structure 
which are also prosentod in these papers and others (5,61, Most of the proposed 
models, with exception of the model proposed by Panofsky 15], are limited to the 
following conditions: X, Strong, stationary winds, 2. Neutral thermal stability 


146 


H. Wt Tielcman and S. E. Mullins 


3, Level terrain with unform roughness. 

Upon reading of the cited references one does not always get a clear picture 
when neutral conditions exist. Deaves and Harris [6] do not give any conditions 
fcr neutrally stable flow other then the mean velocity has to be strong enough^ 
Counlhan [4] considers only those data which specifically Indicate adiabatic or 
near adiabatic conditions, or which have wind speeds in excess of 5 m/s at a 
heigh of 2 -lOm. ESDU [2,3] considers neutral stability to exist when the hourly 
mean wind at a height of 10 m is greater than 10 m/s. Using the gradient 
Richardson number Ri * (g/0) (d0/dz)/(dU/dz)2 as a measure for the degree of 
thermal stability, near neutral conditions are suggested to exist by Teunlssen [1] 
for 1 ri| <0.03 and by Panofsky [5] for |ri| < 0.01. 

Daytime observations throughout the PBL under moderately strong ylnd conditions 
at Minnesota [7,8] indicate that the PBL consists of a surface layer, a mixed 
layer and an overlying inversion layer delineating the height of the PBL, z., above 
which the flow is relatively undisturbed by the presence of the earth's surface 
(free atmosphere) . The height of the PBL increases during the day as the earth 
surface heats up. The large-scale convective turbulence structure in the PBL 
interacts with the air in the free atmosphere causing the formation of an inver- 
sion through which warmer air of higher velocity and with a clockwise directional 
shear is entrained into the mixed layer [7,8,9] . Full-scale measurements [8], in, 
as well as laboratory simulation [10] of the convective PBL seem to suggest the 
existence of an inverted thermal plume structure throughout the PBL. These large- 
scale structures contribute significantly to the variances of the turblence com- 
ponents down in the surface layer [8,11]. 

Even over the roughest surfaces and under the strongest wind conditions, the wind 
shear, dU/dz, decreases with height and the convective turbulence gains Importance 
with height. The height above which the convective turbulence cannot be neglected 
varies with wind speed, surface roughness, time of day, and a variety of atmosphe- 
ric conditions such as cloud cover, precipitation etc. 

Also deviating from the neutrally stable model is an on-shore wind [12] . As the 
usually stable air passes a coast line, the flow is modified as a result of both 
a change in surface roughness and surface temperature, giving rise to the develop- 
ment of an internal boundary layer (IBL) . Because of the lack of turbulent mixing 
with adjacent layers, the flow in the inversion layer above the IBL can manifest 
itself as a surface Jet with a maximum mean velocity at a height varying be- 
tween 60-200 m [12,13]. 

In this article mean velocity and turbulence measurements which were obtained 
under moderately strong wind conditions varying in direction from southwest to 
north (category I) and from southerly directions (category II) are described and 
compared with empirical relations for the PBL models as listed in the previously 
mentioned references. 

EXPERIMENTAL DETAILS AND DATA HANDLING 

Mean wind and turbulence measurements have been obtained with two types of 
instrumentation from the 76-meter (250 foot) meteorological tower at Wallops 
Island, Virginia. Wallops Island is one of the barrier islands protecting the 
so-called "Delmarva” peninsula from the Atlantic Ocean. The island, which consists 
of a narrow strip of dunes, approximately 3 meters above sea level is situated 
in a northeast-southwest dire**tion and is separated from the peninsula by a tidal 
marsh about three kilometers wide (Fig. 1). The island is used by the National 
Aeronautics and Space Administration as a sounding rocket launch facility. 
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Profile Instrumentation 


Measurements were obtained with cup-vane velocity instruments and aspirated tem- 
perature sensing probes at 5 levels mounted at 15 m (50 feet) Intervals. This 
system is primarily used for profile measurements » the outputs are sampled once 
every two seconds. The data are analyzed in blocks of 2® - 256 samples repre- 
senting a data record of 512 seconds in duration. Inspection of block means 
leads to the selection of stationary data samples usually six to ten blocks in 
length (51 to 85 min). Sample averages of velocity » direction and temperature 
are obtained as well as variances and covariances of the horizontal components. 

In addition, a 256-pQint autocorrelation function of the u component is calculated 
from which the turbulence integral scale, is obtained. Data from this in- 
strumentation system can only be analyzed in the frequency range of 0.00195 - 
0.25 Hz. 

Turbulence Instrumentation 

Six three-dimensional split-film anemometers (TSI-1080 D) are used for turbulence 
measurements. These instruments are located at the same levels as the cup-vane 
instruments and also at the 9 m (30 foot) level. The electronics for the split- 
film anemometers as well as the data- acquisition system are located in an instru- 
ment trailer which is parked at the base of the tower. For a detailed description 
of the calibration procedures and the data-acquisition system consult references 14 
and 15. The digitized data (sample rate 200 Hz) are recorded and consequently 
analyzed on an IBM-370 computer. Four main computer programs have been developed 
to calculate the following statistical parameters describing the turbulent flow: 
mean values, variances and covariances, spectra, cross spectra and integral 
scales [16]. 

The frequency range in which the data are analyzed is covered in three stages. 

The highest range (0,0244 - 100 Hz) is obtained by dividing each data record into 
80 consecutive blocks of 8192 data points each. For the middle range (0.00153 - 
6.25 Hz) and the lowest range (0.00031 - 1.25 Hz), the original data are subjected 
to a 16-polnt and an 80-point non-overlapping average respectively. The spectra 
computed from the averaged data inherently show some scatter in the low-frequency 
range but the agreement in the overlap regions is excellent. 

RESULTS AND DISCUSSION 


Terrain Description 

Winds corresponding to category I (west-north) are usually encountered following 
the passage of a cold front and will have crossed in succession the mainland, the 
Chesapeake Bay (20-50 km), the ”Delmarva*' peninsula (20-50 km) and the tidal 
marsh (3-5 km). Depending on the direction, for the last 200-300 m the air travels 
over dry land before it arrives at the tower location (Fig. 1). The tidal marsh 
consists of shallow areas of water interchanged with swamp vegetation, mostly 
grass of a maximum height of 1 m. Some taller vegetation consisting of bushes and 
brush of a maximum height of 5 m exists in several upstream directions. For winds 
in the sector between 20 degrees south of west and west, a 6.5 m high rocket fuel 
storage bunker is approximately 85 m upstream. An elevated roadway 2 m above the 
surrounding terrain passes the tower on the west side with 200 m. 

f . 

Southerly winds (category II) approach the tower from the Atlantic Ocean and cross 
the beach approximately 300 m upstream of the tower. In the southerly direciton a 
6 m high rocket assembly building is about 100 m upstream possibly affecting the 
measurements at the two lower levels. Comparison of the measurements from another 
meteorological tower (91-m) at Wallops Island shows that the effects on the mean 


ORIGINAL PASS fg 
OF POOR QUALiTY 


148 


H. W. Tielemaiv and S. E. Mullins 


flow from a 30-m tall rocket launching tower about 300 m upstream In a direction 
15 degrees west of south can be considered to be negligible. 



Fig. 1. Tower location on V/allops Island and coordinate system orientation 


Category I Winds 

Mean Wind and Mean Temperature Profiles 


For flow over homogeneous terrain under near neutral conditions, the velocity pro- 
file should obey the logarithmic law: U = (U*/ic) IhCz/zq). In practical situations 

these conditions are seldom met and profile variations usually exist as a result 
of non-neutral conditions and/or non-uniform upstream terrain. 


The logarithmic velocity profiles of category I winds over near-uniform terrain 
are either linear within the height of observation (76 m) or show two different 
linear segments (Fig. 2). All upper segments show a tendency towards a uniform 
profile. The heights of the kinks (intersection of the two linear segments) in 
the logarithmic velocity profiles seem to vary with the degree of thermal stabi- 
lity in the surface layer. However, values of the profile shear-velocity, Ug, 
obtained from the lower linear segments compare quite well with the shear 
velocity obtained from turbulence measurements in the surface layer. These ob- 
servations were made in the stability range Ri > -0.06 with Ri evaluated at 
15 m and with the shear velocity based on the horizontal shear stress according 
to U* = [(uw)2 + (vw)2]l/^ and analyzed in the high-frequency range (0.0244 - 
100 Hz). The implications of these observations are that in this stability range, 
values of z can be obtained which correspond to those obtained under neutral 
conditions and that the low-frequency velocity fluctuation components do not 
affect the velocity profile in the surface layer. 


Although the scatter of the data is appreciable, profile results Indicate that 
for the stability range Ri > -0.15, the average value Zq = 0.034 m corresponds 
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Fig. 2. Typical category I logarithmic 
wind profiles 



Fig. 3. The roughness parameter, 

versus Ri at z « 15.2 m for 
category I winds 


quite well with predicted value from ESDU [2] for similar terrain (Fig. 3). For 
Ri < -0*15, Zq decreases rapidly as the surface layer becomes more unstable. 
Figure 4 shows the variation of z with wind direction, <j>, for category I winds, 
clearly indicating the location or the bunker at » 90 ** and increased roughness 
for ^ > 170*. 



Fig. 4. The roughness parameter, z , versus wind direction for 
category I winds 


For non-neutral conditions the velocity profile is given by U*=(U*/k)[ln(z/z 
where is a universal function depending on the stability parameter z/L which 
in turn depends on Ri [5]. Theoretically, if the velocity, U, is plotted as a 
function of Inz - values of Zq corresponding to neutral conditions could be 
obtained. This technique applied to the category 1 wind profiles at Wallops 
Island leads to large corrections resulting in unrealistic estimates 
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Under moderately strong wind conditions » typical daytime temperature profiles 
have a maximum negative gradient near the surface and typical nighttime 
tei!n>erature profiles show a maximum positive gradient near the surface 
both decreasing in magnitude with height. During the evening and night the 
temperature profiles often show kinks at heights varying between 15 and 45 m. 

No significant variation in z can be detected for daytime or nighttime data 
provided they fall in the stability range -0.15 < Ri <0.1 with Ri measured 
at z * 15 m. 

Turbulent Stress and Turbulent Heat Flux 

The horizontal turbulent shear » [(uw)^ *f (vw)2j^/^, is nearly constant throughout 
the height of the lower linear segment of the logarithmic velocity profile. The 
magnitude of the individual components vary considerably with height, although uw 
is consistently negative, vw has been observed to be either positive or negative. 
When analyzed in the middle- frequency range, the magnitude of either shear com- 
ponent is 25% larger In comparison with the high-frequency analysis. The local 
shear velocity, U*, defined as the square root of the horizontal turbulent stress 
obtained at the lowest three observation levels and analyzed in the high-frequency 
range compares well with the profile friction velocity, Ug. 

Day time observations show negative values for the vertical turbulent heat flux, 
w0, while the horizontal components (u6 and v6) are either positive or negative. 
The simultaneous observations of a downward turbulent heat flux and a negative 
temperature profile lead to the problem that the gradient Richardson numbers (Ri) 
and flux Richardson numbers (R^) have opposite signs. A possible explanation 
for this phenomena is that category I winds experience a discontinuity in sur- 
face heat flux while crossing the Chesapeake Bay and the tidal marsh. Conse- 
quently, the downward heat flux as the result of entrainment of warmer air from 
the free atmosphere can be expected to penetrate deeper into the PBL giving rise 
to the counter-gradient heat flux. 


Variances 

Standard deviations of the horizontal velocity components (a , a ) were obtained 
with the cup-vane instruments in the frequency range of 0.005f95-8.25 Hz, and of 
all components from the hot-film probes in the middle frequency range of 0.00153- 
6.25 Hz, The contributions to the variances outside the middle frequency range are 
negligible. Because of the limited frequency response of the cup anemometers, 

0 from the cup-vane system had to be adjusted by a factor of 1.17, which was based 
on comparison with o obtained from the hot-film instruments. For logarithmic 
velocity profiles wi¥h a kink, the average value A = without a 

kink in the profile, the average value of A ■* a /u* *= 2.9, measured at the 15-m 
level and decreasing to 2.2 and 2.7 respectively at the 76-m level. Average 
values of A = a /U* vary from 2.5 at the 15-m level to 2.3 at the 76-m level, 
while the average value of A^ = ^w^^o “ uniform over the height of obser- 

vation. These results obtained near°the surface are in fair agreement with the 
average values of these coefficients obtained from several sources [1,4 and 5] of 
2.5, 2.0 and 1.2 for A , A^ and A respectively. The variances of the horizontal 
velocity components show a^great deal of variability, and on the average the 
ratio A of the lateral velocity (wind direction) at Wallops Island is 
larger ^han at other sites. Inspection of several years of strip-chart recordings 
of wind speed and ditection obtained from a propeller-vane anemometer at Wallops 
Island Indicates that for NW winds, the wind direction shows frequently large 
fluctuations towards the north, explaining the large values of A^ (Fig. 5). 
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Theta rtcordlngi indicate that the fluctuations in wind diraction are larger than 
normally eneountarad and moreover, are not normally distributed which exii/lains 
the non~zero turbulence stress » w. 

The average turbulence intensities are in fair agreement with those suggested by 
ESDU [3], with the expect ion of the lateral turbulence intensities, a /U, whose 
measured values are consistently higher (Fig, 6). ^ 



Fig. 5i Typical records of wind 
direction and wind speed 
for category I winds 



Fig. 6. Variat;ipn of turbulence In- 
censlty (lateral component) 
with height for category X 
winds 


Velocity Spectra 

For frequencies higher than n *= 3U/z, measured logarithmic velocity spectra vary 
as n*“2/3 according to the Kolmogorov's law for the inertial subrange: 
kS (k) = K Here k is the wave number in the streamwlse direction 

(k=2TTn/U by Taylor's hypothesis) and e is the dissipation rate. With the 
universal constants * 0.5 and introducing the dimensionless 

dissipation, < z e/(U*)3, the^inertiai subrange logarithmic spectra nor- 

malized with the variances have the following form; 

nS (n)/o^ - (0.27 /aJ) and nS (n)/oJ = (0.36/A^)4 for a - v,w 

where f is the dimensionless frequency f = nz/U. With ^ * 1 (assuming near 
neutral conditions) and with values of A , based on either the local shear 
velociiiy, U*, or the profile shear velocity, U*, the Kolmogorov spectrum functions 
represent the inertial subrange u, v and w spectra quite well for all elevations 
up to 76m. 


For frequencies lower. than n=3U/z, the observed, u, v and w spectra vary a great 
deal so that no simple spectrum expression seem to exist for this frequency 
range. Kaimal et al [17] obtained empirical spectrum formulae based on measurements 
in the near-neutral surface layer over smooth and homogeneous terrain. These 
spectrum equations represent the spectra at near-neutral conditions approached 
from the stable regime (O'*’). The Kansas results [17J show that the near neutral 
spectra of the horizontal components (u, v) approached from the unstable regime 
(0 ) have much higher spectral values at low frequencies than predicted by the 
Kaimal spectrum expressions. The discontinuity observed in the u, v spectra at 
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near neutral conditions does not exist for the w componentt and consequently the 
Kaimal spectrum expression should fit all w spectra at near neutral conditions. 

Consequently^ for near neutral conditions, u, v spectra as well as the variances 
show a great variability depending on the presence of low-frequency components. 

The latter seem to be the result of a large-scale turbulence structure observed 
in the mixed layer [8] during the daytime and genetrating down into the surface 
layer. The near neutral u, v spectra in the 0 category do not follow a single 
similarity law. In the high-frequency range the Kolmogorov law vith the height z 
as length scale is appropriate. On the other hand in the low-frequency range the 
height of the mixed layer, z^, seems to be the appropriate length scale [11]. 

All category I daytime u, v spectra belong to the 0* category while only one 
data run (run 16) taken during the early evening belongs to the 0*^ category. 

The difference in the near-neutral u and v spectra obtained in the surface layer 
and belonging to either the 0* or the 0**’ categories is shown in figure 7, which 
also shows the excellent fit of the Kaimal expressions for the u, v spectra in 
the O'*' category and for w spectra from either category. All the measured spectra 
in the surface layer and up to a height of 76m, follow the Kolmogorov law or the 
Kaimal expressions in the high-frequency range, although the constants often need 
to be adjusted. Based on these observations it cannot be automatically assumed 
that under tooderately and possibly strong wind conditions the horizontal turbu- 
lence components give rise to spectra belonging to the 0'*' category. 

Integral Scales 

The longitudinal scales, for a » u, v, w, in the direction of the mean flow are 
often calculated from the normalized autocorrelation function, assuming that 
Taylor's hypothesis is valid. A second method for obtaining these scales is using 
the best fit of measured velocity spectra to the corresponding empirical von Karman 
spectrum functions. 

Figure 8 shows the measured scales, by these two methods for two different 
data runs. One data run (run 18) wltii, and a second run (run 16) without low- 
frequency velocity velocity fluctuations. In the case, low-frequency velocity 
fluctuations are present, the von Karman spectrum functions represent the data 
only in the high-frequency range. Consequently, Integral scales obtained with 
the second method are based on the turbulence in the inertial subrange, excluding 
the contributions of the low-frequency fluctuations which are Included when 
integral scales are obtained from the auto-correlatlon function. In the case 
the low-frequency fluctuations are absent the two methods predict almost equal 
scales (Fig. 8). The ESDU [3) predicted scales for corresponding terrain fall 
in between the measured scales for the two different data runs. 

Probability Density Function 

The probability density functions of each of the velocity components are usually 
considered to be Gaussian in the range of + 3a^ [3,4]. Careful Inspection of 
many hours of category I wind data show large asymmetrical fluctuations in 
direction towards the north (Fig. 5). Measured probability density functions 
indicate that the fluctuations of the velocity magnitude in contrast to the 
direction fluctuations are usually normally distributed. This was especially 
observed for daytime data, while for the evening run without low-frequency 
fluctuations (run 16) both magnitude and direction are almost normally distri- 
buted. Several probability density functions of direction fluctuations show a 
secondary peak or the nucleus of a secondary peak at -1.5o (Fig. 9). These 
peaks could be caused by a downward flow of air with a more northerly direction 
probably as the result of entrainment. 
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Fig. 7a. 


U- spectrum» run 16, z»9.1 m, 
U»6.23 m/sec a •0.91 in/sec 
Kaimal [17]. 


Fig. 7b. V- spectrum, run 16, 
2*9.1 m, U“6.23 m/sec 
0^ *0.59 m/sec 

Kaimal [17]. 


I. ■ 
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Fig. 7c. W- spectrum, run 16, z*9.1m, 
U»6. 23 m/sec o *0.53 m/sec 

w 

Kairaal [17]. 


Fi|5. 7d. spectrum, run 18, 
2*9.1, U*10.l4 m/sec 
0^*188 m/sec 

Kaimal [17]. 
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Fig. 7e. V- spectrumi run 18, z“9.1 m, Fig. 7f. W- spectrum, run 18, z»9.1 m 
U“10.14 m/scc c 7„»1.42 m/sec U*10.14 m/sec a «0.75 m/sec 

V w 


Kalmal [17]. 


Kaimal [17]. 


-^ESDU [3] 

O RUN 18 METHOD I 
• RUN 18 METHOD 2 
RUN 16 METHOD I 
A RUN 16 METHOOt 
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Fig. 8. Integral scales obtained from 
the autocorrelation function 
(method 1) and from the 
von Karman spectrum f unction 
(method 2) for two runs, one 
with low- frequency components 
(run 18) and one without low- 
frequ’ency components (run 16) . 


Fig. 9 . Probability density function 
of the wind direction fluctu 
ations of run 18, 76 m level 
o^“7.97®. Solid line repre- 
sents the Gaussian distribu- 
tion function. 
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Category II Winds 

Mean Wind and Mean Temperature Profiles 

The category II velocity profiles when plotted in the usual logarithmic 
fashic£ 3 < tre either straight or curved opposite to the category X logarithmic 
profile or show a maximum below the 76-m level (Fig. 10) • The roughness lengths 
obtained from either the straight logarithmic profiles or from the lowest two 
observations of the curved profiles vary a great deal and a::e too high for the 
terrain at Wallops Island. Corrections of these pro files » in order to obtain 
values of corresponding to the terrain » as suggested by Panofsky [5] are 
inadequate. 

The mean* temperature profiles show some unusual kinks mostly between the 30 
and 45 meter levels. The most interesting sequence of mean*teiQperature distri* 
butions is from three runs taken on one suimncr afternoon (Fig. 11). The early 
afternoon data show a decrease in temperature up to 30 m with a rather strong 
positive tempersture gradient at higher elevations. As the afternoon progresses 
this temperature gradient decreases drastically. The majority of observations 
show a positive temperature gradient: above 30 m of varying intensity and below 
30 m the temperature gradient is elcher positive or negative depending on the 
time of day and the time of year. 

Turbulent Stress and Turbulent Heat Flux 

The measured turbulent stress uw has a maximum value at the 9 m level t and de- 
creases rapidly with height. The vertical turbulent heat flux w 6 Is positive 
(upward) below 30 m and is negative (downward) at this height and above, corres- 
ponding to the observed summer tf?mperature profiles. These results Indicate a dual 
type bourdary layer within the height of observation (76 m) , with unstable 
conditions near the surface and an overlying inversion. The unstable layer near 
the surface corresponds to the internal boundary layer which deve 7 tops as the 
stable on-shs^re wind is heated from below after it crosses the bf ch. The depth 
of the IBL Increases with distance from the shoreline. 

Variances 

The standard deviations of the velocity fluctuations, are largest In the IBL 
and decrease rapidly in the overlying inversion. The ratios, a^/U*, are 
quite sensitive to stability and increase for deviations from neutral conditions 
in both the stable and unstable regime [12]. For near-neutral conditions In the 
IBL, the average values of o are 1.75, 1.5 and 1.4 for d ** u, v and w res- 
pectively. Turbulence intensltlc,*; O /U in the stable air above the IBL (25 m or 
higher) are about half of tVte ESDU [ 3 ] predicted values with * 0.001 m (Fig. 12) 

In the IBL with z ■ 0^01 m, the ESDU predictions [3] are comparable to the 
measurements. The hot- film data shown in Figure 12 were acquired during hot 
and sunny summer afternoons with mean velocities at the 76 m level varying 
from 6.7 m/sec to 12.7 m/sec. The cup-vane data were obtained during raornihg 
and early afternoon in the winter and early spring, with mean velocities at the 
76 m level varying from 10.7 m/sec to 20.6 m/sec. For 15 of a total of 38 runs 
investigated* the streamwise turbulence intensity at the 76 m level is less than 
3%. The lateral and vertical turbulence intensities vary in a similar fashion. 

Included in Figure 12 are the turbulence intensities for run 7 obtained under 
extreme stable conditions (R^ +25,5, U - 10.7 m/sec at z=61 m) during a sunny 
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Fig. 10. Typical category II logarithmic wind profiles. 




Fig. 11. Temperature profiles for 
runs 7.8 and 9. 


Fig. 12. Variation of turbulence in- 
tensity (streamwise compo- 
nent) with height for 
category II winds. 
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warm afternoon in August. Under these conditions the turbulence almost disap- 
pears. During a trial run with the hot-film system under similar conditions, 
the turbulence disappeared completely at the two highest observation levels for 
at least 20 cnlnutes. For many of the runs, a distinct clocHwiae rotation of a 
maximum of 20 degrees of the mean velocity vector with height is observed. With 
stable sir coming from the ocean, maximum velocities at elevations below 200 m 
have been observed visually from exhaust trails of launched rockets and occasion- 
ally measured with the tower instruments below 76 m. (Fig. 10). 

Based on these observations it can be concluded that the stable air above the 
IBL with low turbulence Intensities and small turbulence integral scales, loses 
dynamic contact with the air above and below it. Little mixing and little momentum 
transfer takes place between the stable air and the adjacent air so that Coriolis 
forces become important and the stable air encoimters little resistance resulting 
in the low-level Jet. 

Velocity Spectra 

-2/3 

Measured logarltlimic velocity spectra of category II winds vary as n in the 
inertial subrange and those spectra originating in the IBL belong to the 
category, varying as in the low frequency range. However, for elevations 
above the IBL, low-frequency components are present in all three components. 

Under stable conditions, wave phenomena have been observed near the surface 
[13, 18] and can be recognized from the turbulence for those flows where the 
turbulence levels are low. Waves are not random and usually occur at a different 
frequency band so that a distinct spectral gap exists between the spectral peak 
due to the turbulence and the spectral peak at lower frequencies due to the waves. 

The part of the spectra due to the turbulence in both the IBL and overlying inver- 
sion are well represented by the Kaimal expressions, although the constants in 
these expressions vary with stability [12], The presence of the low-frequency 
waves varies a great deal with atmospheric conditions and their effect on the 
spectra seems to grow with height, as the fluctuations due to turbulence and wave 
motion fall in similar frequency bands. 

Integral Scales 

Integral scales obtained from the integration of the normalized autocorrelation 
functions vary a great deal with the presence of the low-frequency components^ 

Values of the integral scales obtained with this method for the consecutive 
afternoon runs 7, 8 and 9, measured at the 76-m level, vary from 406, 162 to 
56 m respectively. The wavelengths corresponding to the peak of the logarithmic 
spectra vary substantially with stability and^height [12]. For near-neutral 
conditions in the IBL, the peak-wave length equals z/0.15, z/0.35 and z/0.15 
for a » u,v and w respectively. Normalized spe^^tral values nS (n)/U*^ evaluated 
at f“10, vary a great deal with stability [12]. Consequently, ^integral scales 
obtained from fitting of the von Karman functions vary with stability as well. 

CONCLUDING REMARKS 

The wind measurements made at this eastern coastal region under sustained moderately 
strong wind conditions show marked deviations from observations obtained in purely 
shear-generated turbulent flow. For near-neutral conditions of cateogry 1 (west- 
north) winds, low-frequency velocity and direction fluctuations Influence all 
parameters describing the turbulence. The large direction fluctuations can have 
a pronounced effect on the response of tall buildings especially in the case the 
mean wind blows diagonally onto the building [19,20], In wind tunnel model tests 
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the non-mechanlcally produced turbulence cannot be Rlmulated adequately re* 
suiting in a discrepancy between model and tests [20], 

For moderately on-shore winds from southerly directions^ a surface based in- 
version is present even with the hourly mean wind speed over 20 m/sec and a 
maximum velocity (low-level jet) may exist at heights between 60-200 m. The 
turbulence intensities in the stable layer are often less than 3% and under cer- 
tain atmospheric conditions the turbulence vanishes completely. Such flow could 
affect the response of suspension bridges located near coastal areas especially 
when the bridge decks are at the sane height as the almost turbulence-free low- 
level jet. 

Based on the observations made at Wallops Island under moderately strong wind 
conditions, the turbulence is not of purely mechanical origin. Moreover, as a 
result of non-homogeneous upstream surface conditions at this site, similarity 
theories are not applicable. 
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ON THE WIND STRUCTURE NEAR THE SURFACE AT A 
MID-ATLANTIC COASTAL SITE 

by 

Henry W. Tieleman* 


Introduction 

This paper describes the mean flow and turbulence measurements 
made in the atmospheric surface layer up to a height of 76 m at 
Wallops Island, Virginia. The observations were made under moderately 
strong wind conditions with an hourly mean-wind velocity of at least 10 
m/s at a height of 76 m. The measurements were made with cup-vane and 
resistance temperature sensors primarily for profile information and 
with hot-film and thermocouple instruments primarily for turbulence 
measurements. [1,2] 

Site Description 

Wallops Island consists of a narrow strip of dunes, approximately 
3 m above sea level, and is situated in a northeast-southwest 
direction. The island is separated from the **Delmarva” peninsula by a 
tidal marsh on the west side, and with the Atlantic Ocean on the east. 
The terrain roughness at this site (Fig. 1) is not nearly as uniform ac 
for other sites usually selected for ground wind experiments. 

At the Wallops site, even under the strongest observed wind 
conditions, neutral thermal stratification has never been observed for 
any length of time. For most directions the atmospheric winds 
experience a change in surface roughness and surface temperature before 
the tower location is reached. 

*Professor, Department of Engineering Science and Mechanics , Virginia 

Polytechnic Institute and State University , Blacksburg, VA 24061. 
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Discussion of Results 

The diabetic wind profile 

U » (Uq*/K) [ln(z/zQ - >l;(z/L)l (1) 

is used to determine the roughness length, zq and the profile friction 
velocity Uq*, as outlined by Panofsky [3). Figure 2 shows the 
variation of roughness length, Zq, averaged for each 15* wind sector. 

An insufficient number of records in the wind sector 75<4‘''i95*, and 
upstream buildings made it impossible to get consistent results in this 
sector. In Figure 3 the profile friction velocity, Ug*, is compared 
with the friction velocity, U*, obtained from the Reynolds stresses 
either as U*= /3w or U* = [uw^ + In general the values of Ug* 

are larger than those for U*. For northwesterly winds Ug* compares 
reasonably well with the friction velocity, U*, based on the total 
horizontal Reynolds stress. 

Observed values of o^/U* for a = u,v and w are consistently lower 
than the usual near-neutral predictions of o^/U*-2.5, o^/U*=l.9 and 

0 /U**1.3. Figure 4 shows the average ratios of standard deviations 

w 

a /a and a Jo for each 15* wind sector. The observations of a Jo 

V u w u V u 

for wind directions parallel to the island are less than the expected 

value of 0,8 and for south winds and northwest winds larger than 0.8, 

Under convective conditions, occurring frequently during daytime for 

northwest winds, a^and are approximately equal except at the lowest 

observation level, but the ratios ojo are close to the expected value 

w U ' 

of 0.53 (Table 1). For southwest vnnds both ratios a Jo and a Jo 

V u w u 

conform closer to their respective expected values of 0.8 and 0.53. 
Variations of average turbulence intensities of the horizontal velocity 
components obtained with the cup-vane system are shown together with 
their standard deviations and their extreme values for the highest 
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observation level (2*76,2m) in Figures 5 and 6. In the wind sectors 
90*<^<165* (on-shore winds) and 0®<<j)<15* an insufficient number of 
records were available and only average turbulence intensities are 
shown. Note the low turbulence intensities for southerly winds usually 
observed under stable conditions above the internal boundary layer. A 
near constant turbulence intensity o^/U of 14% is observed in the wind 
sector 225®<(|)<360'^ , in contrast with o^/U which reaches a maximum value 
of 147 j for a wind direction of 300*. This observation is consistent 
with the maximum observed value for the ratio in Figure 4 for the 

same wind direction. 

For near-neutral stability conditions a universal form for all u,v 
and w spectra has been proposed by Kaimal [4], 


nS^ (n) 0.164 f/f^ 

* — — for a = u, V and w, (2) 

0 2 l+0.164(f/f„)5/3 

a u 

where f = nz/U and fg is the scaling parameter which is the intercept 
on the reduced frequency axis of the extrapolated inertial subrange and 
the line nS^(n)/a^^ = 1. The von Karraan spectrum functions, which are 
frequently used by wind engineers may also be modified into the Kaimal 
format as: 


nS^(n) 0.156 f/fp 

aj- [l+0.108(f/fQ)2]5/6 

( 3 ) 

and 

nS„(n) l+0.679(f/f.)2 

= 0.12(f/f.) [- 9 . ] 

U+0.255(f/fQ)2}ll/S 


a=v , w 
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The two sets of spectral functions are nearly identical except for the 
values of the constants, and the Kaimal spectrum is slightly flatter 
near the spectral peak. In the absence of low-frequency velocity 
fluctuations, the shape of the observed spectra of the three wind 
components for southerly (on-shore) winds are remarkably similar and 
insensitive to stability. Either the Kaimal (2) or the von Karman (3) 
spectrum functions fit the observed spectra extremely well (Figs. 7,8 
and 9). However in the absence of low-frequency components, some of 
the V and w spectra observed at the lower elevations show that the 
Kaimal and the von Karman spectrum functions overestimate the spectral 
densities in the low-frequency range and that the spectral slope is 
larger than the predLcted value of 1.0 (Fig. 10). A secondary spectral 
peak may develop in the spectra of all three components as the result 
of low-frequency internal gravity waves (Fig. 11), 

For land winds from westerly directions similar observations have 
been made under generally stable conditions with secondary low- 
frequency peaks. Also, for some of the v and w spectra, the Kaimal and 
the von Karman spectral functions overestimate the observed spectral 
densities in the low frequency range, especially at the lower 
elevations. Under convective conditions all spectra show varying 
degrees of increased contribution to the spectral densities in the low- 
frequency range, specifically for the u and v spectra. In many cases, 
the w spectra show also an increase in spectral densities at lower 
frequencies, often increasing with height. In Figures 12 and 13 the u, 
V and w spectra are shown for z - 9. 1 m and z ” 76 m respectively in 
comparison with the Kaimal spectrum function (2) . The results also 
show larger increased spectral densities in the low frequency range for 
the observations at z=76m than for the observation at z=9.1m. 


In the high frequency range all observed logarithmic spectra vary 
with consistent with the Kolmogorov's law for the inertial 

subrange. In this frequency range, the v and w spectra densities are 
to be higher than the spectral densities of the u component by a factor 
of 4/3 as predicted by isotropy. Nevertheless, for most observed 
spectra at the Wallops site under cither stable or convective 
conditions, the u, v and w spectral densities are the same in the **2/3 
range, especially for the observations at the lower levels. For some 
of the spectra observed at the higher elevations either the v or the w 
spectral values are higher, but the expected spectral ratios 

®w^"^ 4 

■ 7 T~t — A — “ t: — 7 — ^ have not been observed for any of the sets of 

S (n) S (n) 3 ^ 

u u 

measured spectra including those shown in Figures 12 and 13. 

Conclusions 

Based on the presented observations, it can be concluded that 
there is no simple PBL flow model available to describe the mean and 
turbulent flow at the Wallops site under moderately strong wind 
conditions. Statistical parameters describing the turbulence vary a 
great deal with the presence or absence of low-frequency velocity 
fluctuations. The atmospheric surface flow at this site does not 
follow the Monin-Obukhov similarity or mixed layer similarity 
consistently, mainly due to the nonuniformity of the surface roughness 
and surface temperature. Based on the measurements made at this site 
under the strongest observed wind conditions, there is no evidence that 
the PBL flow will approach the neutral boundary-layer model under 
extreme wind conditions . 
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TABLE Ir Averages of their respective standard 

deviations for north-west winds and south-west winds obtained 
from hot-film measurements. 
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For neutral conditions without low-frequency velocity fluctuation 
predictions are <^yt/o = 0.8 and o^/o = 0.53. 
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Fig. 1. Wallops Island, Immediate Surroundings of the Meteorological Tower 
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EXECUTIVE SUMMARY 


Detailed measurements of the mean flow and turbulence have been 
made with the use of a mlcrometeorologlcal facility consisting of an 
Instrumented 76-m -tall tower located within a 100-m distance from 
the Atlantic Ocean at Wallops Island » Virginia. An Interpretation of 
the experimental results demonstrates that under moderately strong 
wind conditions (hourly mean wind speed between 10 m/s and 20 m/s at 
a height of 10 m) , the popular neutral boundary-layer flow model fails 
to provide an adequate description of the actual flow. 

For daytime westerly winds the convective boundary layer, which 
has been previously observed at sites on the continent, provides an 

adequate model for the surface flow at the Wallops Island site. 

• 

However variations ftom this model have been observed for certain 
wind directions and under certain atmospheric conditions such as 
low altitude cloud cover combined with precipitation. The observed 
low-frequency velocity fluctuations give rise to Increased turbulent 
Intensities and larger turbulence Integral scales. These low- 
frequency fluctuations also occur In the surface layer where the 
observed mean velocity profiles generally fit the logarithmic law quite 
well. 

For on-shore winds the surface flow is complicated as the re- 
sult of the development of an internal boundary layer (IBL) as the air 
crossing the beach generally experiences a change In surface roughness 
and surface temperature. The Internal boundary layer has a height 
between 15 m and 30 m at the tower location depending on wind direction 
and change In surface conditions. For southerly winds the warmer air 

v 
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flows over the cooler water allowing the existence of a surface-based 
inversion of \'ariable depth. Under these conditions a low-altitude 
maximum '.'eloclty (surface jet), occasionally below the highest 
observation level of 76 m, has been observed* Under extreme stable 
conditions at hourly mean velocities in excess of 10 m/s the tur- 
bulence has been observed to vanish completely. In addition, low- 
frequency Internal gravity waves have been observed to co-exist 
with the turbulence. 

In addition to detailed flow lnforma;,ion for all wind directions, 
averages of the Important .flow parameters used for design such as 
vertical distribution of mean velocity, turbulence intensities and 
turbulence Integral scales have been presented for wind-direction 
sectors with near-uniform, upstream terrain. Power spectra of the 
three velocity components for the prevailing northwesterly and southerly 
winds are presented and discussed in detail. 

The experimental results indicate clearly that the non-uniformity 
of the upstream surface conditions, the non-neutral thermal stratifi- 
cation and the presence of appreciable low-frequency velocity fluctu- 
ations have a pronounced effect on the surface flow. Consequently It 
is impossible to find a simple and single PBL model to describe the 
flow at this site even under moderately strong wind conditions. More- 
over, there is no evidence that under still stronger wind conditions 
(hourly mean wind speed at z-10 m over 20 m/s) the surface flow will 
alter sufficiently as to conform to the neutral boundary-layer model 
whose turbulence is of purely mechanical origin. 
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1. INTRODUCTION 

The purpose of this report Is to provide information on the local 
wind climate at a mid-Atlantic coastal Site. The acquired information 
can be used for the design of wind-turbine generators at similar sites. 
Since wind is a very Important design parameter for these generators, 
information is provided in this report on wind speed, wind direction, 
wind shear and wind turbulence. 

The data presented in this report were collected from an instru- 
mented meteorological tower, 76.2 m (250 feet) tall and located at 
Wallops Island. This Island is one of the barrier Islands at the 
Atlantic coast along the Eastern Shore of Virginia and is used by 
the National Aeronautics and Space Administration as a sounding rocket 
launch facility. The cesults acquired from this facility should be 
typical for any Atlantic coastal site, although local effects such as 
upstream buildings and obstacles and changes in surface roughness and 
surface temperature modify the flow near the surface>. At Wallops 
Island, the surrounding terrain beyond a distance of 100-300 m from 
the cower can be considered as homogeneous and uniform, so that the 
flow above a height of approximately 10-30 m should not be affected 
by local terrain nonuniformities. 

The wind and temperature data from this site were acquired under 
moderately strong wind conditions with an hourly mean velocity of at 
least 10 m/s at Che 76. 2m (250 ft) level. For wind directions be- 
tween northeast and southwest this requirement had to be reduced to 
approximately 8 m/s, since strong winds from this sector occur very seldom. 


Mean wind and turbulence measurements were made with two types of 
Instrumentation consisting of cup-vanes and resistance temperature 
probes primarily used for mean profile measurements of velocity and 
temperature respectively. In addition^ the cup-vane instruments were 
used for turbulence intensities of the two horizontal velocity compo- 
nents and horizontal and vertical turbulence integral scales. The hot- 
film, and thermocouple system was used for measurement of turbulence 
intensities, turbulence fluxes and velocity spectra in all three 
directions. The cup-vane system was used to acquire wind data from 
all directions, while the hot-film system was only used for turbulence 
measurements from the two prevailing wind directions, south and north- 
west. 

The results of this experimental research are presented in a form 
suitable for design purposes. Where ever possible the results are also 
compared with previously published results and with existing empirical 
models for near-neutrally stratified low-level winds. 

2. SITE DESCRIPTION 

Wallops Island consists of a narrow strip of dunes, approximately 
3 meters above sea level, and is situated in a northeast-southwest 
direction. The island is separated from the "Delmarva" peninsula by 
a tidal marsh on the west side, and With the Atlantic Ocean on the east. 

Winds with directj.ous varying between west and north are usually 
encountered following the passage of a cold front. Winds from this 
sector will have crossed in succession (Fig. 1) the mainland, the 
Chesapeake BSy .(20-50 km), the "Delmarva" peninsula (20-50 km) and the 
tidal marsh (3-5 km). Depending on the wind direction, for the last 


200-300 m the air travels over land before it arrives at the tower 
location (Fig. 2) . The tidal marsh between the Island and the 
peninsula consists of shallow areas of water Interchanged with swamp 
vegetation, mostly grass of a maximum height of 1 m. Some taller 
vegetation consisting of bushes and brush of a maximum height of 5 m 
exists In several upstream directions. For wind directions between 
255* and 270*, a (>.5 m high rocket fuel storage bunker Is approximately 
90 m upstream (Fig. 2). An elevated roadway (levee) 2 m above the 
surrounding terrain passes the tower on the west side within 200 m. 

Winds with directions varying between north-east and south approach 
the Island from the Atlan|:lc Ocean. 

Sectors with approximately the same Immediate upstream roughness 
have been established as shown In Fig. 3. Between 0* (north) and 
30* the upstream terrain features two bunkers within a distance of 
100 m from the tower. In addition, a few small buildings and Inter- 
mittent patches of brush are upstream as far as 750 m. Between 30* 
and 45* (wind direction parallel to the Island) many buildings are 
upstream and winds in this sector should encounter the roughest terrain 
at this site over a distance of approximately 4 km. For wind directions 
between 45“ and 210* the winds approach the island over water, and cross 
the beach at varying distances from the tower depending on the direction. 
For wind directions between 140“ and 170* a one-story rocket assembly 
building is about 100 m upstream of the tower. The prevailing southerly 
winds vary in direction between 170* and 210", however in the sector 
between 195* and 210* the 45-ra tall Aerobee tower and associated 
buildings are about 300 m upstream. For directions between 210* and 
230*, the wind direction Is approximately parallel to the island 
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with part of the Aerobee tower complex approximately 300 m upstream. 

In this same sector a few other buildings, levees and sand dunes are 
upstream at greater distances. Between 230" and 330" the upstream 
terrala Is very uniform with no big obstacles other than the afore- 
mentioned bunker and roadway. Por the sector between 330" and 360", 
several patches of brush, 2 levees and one radar building are upstream 
of the tower within a distance of 500 m with marsh at further upstream 
distances. 

3. INSTRUMENTATION 
3.1 Cup-Vane Instruments 

The 76 m (250 ft) mlcrometeorologlcal tower is a self-standing non- 
guyed tower with working platforms at 15.2 m (50 ft) intervals (Fig. 4). 
The cup-vane velocity-direction instruments and aspirated temperature 
probes both prli^arlly used for profile measurements are mounted at 5 
levels near each platform. Two sets of cup-vane instruments are mounted 
at each level on 2m booms on opposite sites of the tower (Fig. 4) . 

An automatic electronic switching circuit ensures that data are taken 
only with the Instruments on the upwind side of the tower. The 
electronics associated with this instrumentation system, together with 
a dig;ital readout panel of all instruments from one side of the tower, 
are located in a small instrumentation building at the base of the 
tower (Fig. 4). From this location the digitized data* are transmitted 
to the NASA control center at the main base on the peninsula about 
13 km to the northwest. Here the data from each level sampled at a 
rate of 1 sample each 2 seconds are recorded on digital tape. At 
this sample rate, data can be acquired without interruption for about 
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8 hours. This Instrumentation system Is used by NASA in conjunction 
with its rocket launching operations. Regularly scheduled mainten- 
ance and calibration of this system are performed by personnel under 
NASA's supervision. 

3.2 Hot-film Anemometers 

Six three-dimensional split-film anemometers (TSI-1080D) are 
used for turbulence measurements, which include turbulence intensities, 
turbulence fluxes, spectra and cross spectra of all three turbulence 
components and temperature. These anemometer systems were chosen for this 
research program since they have the advantage of small physical size, 
fast response and high sensitivity over a wide range of velocities. The 
Instruments are mounted on 1.8 m booms at the same levels as the cup- 
vane instruments and also at the 9 m (30 ft) level. Each hot-film probe 
is mounted on a rotor, which is capable of rotating the probe about a 
vertical axis so as to align the probe axis approximately into the mean 
wind direction. The probe-rotor combination is mounted on a 1.8 m-boom, 
which in turn is mounted on the railing at each platform. The probes 
were mounted on the south side of the tower for measurement of the 
prevailing south winds during the summer and on the north side of the 
tower for measurement of the prevailing northwest winds during the 
winter. and spring. The electronics as well as the data-acqulsitlon 
and data-handllng system for this instrumentation system are. located 
in an instrumentation trailer parked at the base of the tower (Fig. 4). 

Each hot-film probe consists of three split-film sensors used for 
measurement of wind speed and direction and a copper-cons tantan 
thermocouple used for temperature measurements . Each sensor consists 
of a 0.15 mm diameter quartz rod coated with a platinum film of about 
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1000 angstrom in thickness. The platinum film on each rod consists 
of two segments, separated from each other by two longitudinal split# 
180** apart. The active elements on each rod are elecCrlcally heated 
to the same constant temperature by separate anemometei circuits. 

The total sensor length is about 5 nun, and the three sensors are mounted 
mutually perpendicular to form a Cartesian coordinate system. When the 
instruments are not used for data acquisition, the three sensors and 
thermocouple are protected by an aluminum shield which can be moved 
pneumatically to cover the sensors. As an added precaution, dry fll- 
te*-®! air is allowed to blow across the sensors when the shield covers 
the sensors. This is done> to protect the sensors from contamination in 
the salt-air environment and moisture while not in operation. For a 
more detailed review of the hot-film anemometer system the reader is 
advised to consult Reference 1. 

Calibration of the hot-film anemometers is carried out in a low- 
speed wind tunnel located at the main base. In order to obtain data of 
a desired accuracy from the hot-film Instrumentation system, a new 
calibration and operating procedure was developed. Instead of using 
the calibration constants supplied by the manufacturer, all constants 
were obtained from calibration procedures carried out in the low-speed 
wind tunnel and thermal chamber. This procedure proved to be both time 
consuming and complicated l>ut necessary. Calibration of each Instrument 
in the wind tunnel was carried out for 11 wind approach angles between 
plus and minus 50® and for 13 velocities in a range varying between 
0.3 and 15 m/s. The best accuracy of the data was obtained for wind 
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directions parallel to the axes of the instrument, and consequently 
the tower mounted instruments were rotated in the direction of the 
mean wind before the data acquisition was started. For details of 
the calibration procedure and the relations for the conversion from 
output voltage to velocity components it is suggested that the reader 
consult References 1 and 2. 

4. DATA HANDLING AND DATA ANALYSIS 
4.1 Cup-vane Instruments 

The output signals from the cup-vane Instruments and temperature 

sensing probes are sampled and digitized at a rate of 1 sample per 

second. This Information is transmitted to the control center of 

the main base, where every other sample is recorded on digital tape. 

The data from these tapes, each capable of storing up to 8 hours of 

data, are then analyzed on the HW-625 computer at NASA, Wallops Flight 
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Center. The data are analyzed in blocks of 2-256 samples, representing 
a data record of 512 seconds. For each sample the east-west and north- 
south velocity components are calculated and averaged over 256 samples 
from which the mean velo>„ity and the mean direction for each block are 
obtained. Also a block mean for the temperature is calculated. 

Reasonably stationary sample records of 5 ts 10 blocks in length are 
sele.cted for further analysis. This selection is based on the inspection 
of the printout of the block means of velocity, direction and temperature 
for all five levels. 

Next the east-west and north-south velocity components and tempera- 
ture are averaged for the selected sample, from which the sample mean 
velocity and sample mean direction are calculated. This direction 
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defines the mean-wind coordinate system with the x-axls parallel to the 
direction of the sample mean wind, the y-axls In the horizontal plane 
perpendicular to the x-axls and the z-axis vertically upward. For all 
the data points In each block the velocity components In the mean-wind 
coordinate system are calculated and averaged to obtain the block means. 
After the block means were removed from each set of components, variances 
and covariances are calculated for each block. Sample varlaiices and 
covariances are obtained by averaging of the block variances and covariances 
over the total number of blocks. The covariances calculated In this manner 
Include all the combinations of like velocity components at the different 
levels, allowing for the calculation of the vertical turbulence Integral 
scales of both the u and v components. In addition, the autocorrelation 
function, of the streamwlse velocity Is calculated from which the 

turbulence Integral scale, L^, Is obtained as follows : 

(o ) 0 

u 

where is the time delay for which the first zero-crossing of the 
calculated autocorrelation function occurs. The turbulence data 
acquired with .the cup-vane system are analyzed in a limited frequency 
range of 0.00195-0.25 Hz. 

A total of 195 digital data tapes were generated during the period 
of July 1974 and December 1978. Approximately 300 data samples were 
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analyzed each varying between 43 and 85 minutes. Initially data were 
acquired with the cup-vane system only, It was not until February 1976 
that the temperature system came on line. However this system Is not 
too reliable and often temperature at one or two levels is missing 
as the result of the equipment being down or out of calibration. Before 
a lightning-arrestor system was installed on the tower, excessive 
amount of damage was inflicted on all systems during thunderstorms as 
a result of line power surges and voltage induction in the cables that 
connect the instruments on the tower to the electronics at the base of 
the tower. During the summer of 1976 a thermograph for recording the 
air temperature at grouted level was added to the system. 

Occasionally when the equipment on the 76m (250 ft) tower was down, 
data-acqulsitlon was switched to the 91 m (300 ft) tower located at the 
north end of the island. This tower is instrumented with cup-vane 
systems at six levels but has no temperature instruments. Its location 
from the beach is 280 m as compared to the 76 m (250 ft) tower which is 
approximately 150 m from the beach. No major buildings or other obstacles 
exist between the 91 m (300 ft) tower and the beach. However, for ocean 
winds the overland distance is longer and more modification of the undis- 
turbed ocean winds can be expected at the 91 m (300 ft) tower. 

4.2 Hot-film Anemometers 

The data-acqulsition and data-handling System is designed to 
handle output from six split-film anemometer systems, sampled at a rate 
of 290 samples per second for a period of approximately one hour. 

This system consists of two main parts: (a) the multiplexing and 

analog recording system and (b) the demultiplexing, digitizing and 
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► 

i 

1 I digital recording system* The seven output voltages from each anemo-*^ 

I y. . ^ , 

1 

j meter are frequency modulated by voltage-controlled oscillators each 

i'j 

jj with a different center frequency. There Is one set of voltage-con- 

5 : ■ : ■ 

,i ttolled oscillators for each probe. • The seven frequency -modulated 

n signals together with a 100 kHz reference signal are fed Into a summing 

! amplifier to produce one single multiplexed signal. The multiplexed 

I signals from each Instrument are recorded on separate channels of an 

! 

I analog tape recorder together with tlme-of-day» which serves as a 

I reference for the recorded data. 

At a later time, each of the multiplexed signals Is demultiplexed 
into Its seven analog components after passage through seven dlscrl- 

: 

j mlnators. In order to avoid aliasing of the velocity spectra the 

I ^ ' ' 

I six output voltages corresponding to the six split films are passed 

through a 100 Hz low-pass filter. Next the analog voltages are sampled 
at a rate of 200 Hz, digitized and recorded on digital tape. 

A mini-computer (DEC Model PDF 11/20) controls the multiplexing analog- 
to-digital conversion and the digital recording. Access to the mini- 
computer is obtained with a teletypewriter. The data conversion starts 
at a time-of-day prescribed by the operator, and the analog-to-digital 
converter performs successive scans and conversions of seven analog vol- 
tages into 16 bit words at a rate of one scan each 5 milliseconds. 

These words are stored in one of the buffers of the mini-computer which 
in turn transfers the data to a 9- track digital magnetic tape. Each 
buffer has a capacity of 209 scans representing 1.05 seconds of data. 

A total of 3300 records make up a single sample record over a time 
period of slightly less than one hour. The tapes with the digitized 
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data are taken to VPI and SU where the data are analyzed on an IBM-3/0 
computer. Four separate computer programs have been developed to 
cnJ.culate the following major statistical parameters: mean values, 

variances, covariance spectra and cross spectra. 

The first step In the data-analysls procedure Is to convert the 
seven output voltages from each film to three velocity components In 
the sensor-oriented coordinate system and temperature, using the con- 
stants obtained from the calibration data. The converted data are 
transferred on another magnetic tape to await the next step of the data 
reduction. 

In the second program velocity and temperature data are analyzed 
13 • 

in blocks of N°°2 =8192 data points, representing nearly A1 seconds 

of data. For each of these blocks of data mean velocity components, 
mean velocity and direction, mean temperature and the four standard 
deviations are calculated. A total of 80 data blocks (almost 55 

minutes) are analyzed In this manner. A statlonarlty trend test Is 

performed on each of the calculated parameters to check for un- 
acceptable nonstationarities. Also inspection af the printout of the 

block parameters helps in the decision whether or not to continue with 

the statistical anlaysls. At this point blocks with unrealistic data 
can be recognized and omitted from the data sample in future analysis. 

The sample mean velocity components are obtained by averaging the block 
means , allowing the calculation of the horizontal angle between the 
sample mean-wind direction and the probe axis. In the following step 
this angle is needed to tranfer the original velocity Components in 
the sensor-oriented coordinate system into u, v and w velocity components 
of the mean-wind coordinate system as defined previously. 
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Block means are calculated for the temperature and for the velocity 
components In the mean-wind coordinates system and removed from the data 
in each block. The resulting fluctuating components are recorded on 
magnetic tape for further analysis. Also variances and all covariances 
for each block are calculated and averaged over 80 blocks to obtain the 
sample variances and covariances. The statistical parameters, including 
the spectra to be calculated in the next step, contain only contributions 
from the fluctuations in the frequency range of 0.0244-100 Hz. In order 
to include contributions from frequencies below 0.0244 Hz, sixteen con- 
secutive data points are averaged into one data point to form a new 
data record, which is also recorded on magnetic tape. This averaging 
is performed after the data are transformed into the mean-wind coordi- 
nates and before the block means are removed. In this way only 5 blocks, 
each 10.92 minutes long, are analyzed allowing for data analysis in the 
frequency range of 0.00153-6.25 Hz. For these new data records, block 
variances and covariances and sample variances and covariances are also 
calculated. For the lowest frequency range the data, after transformation 
into the mean-wind cooordinates, are subjected to an 80-point non-over- 
lapping averaging for analysis in the frequency range between 0.00031- 
1.25 Hz. 

The last step of the data analysis is the spectral analysis of the 
high, middle and low frequency data in the frequency range of 0,0244- 
100 Hz, 0.00153-6.25 Hz and 0.00031-1.25 Hz respectively. Spectral 
estimates are calculated for each block using a specially developed 
Fast Fourier Transform algorithm [4] . The combined averaging 
technique is employed, averaging first all the block estimates at a 
given frequency (ensemble averaging) and then averaging these results 
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over appropriate frequency Intervals (frequency averaging) . 

In total, 24 one-hour data records were generated with the hot- 
fllm system. Nine runs were generated during warm summer, afternoons 
of the year 1976. This set of data was acquired for southerly winds 
only, and the detailed results are presented In Reference 3. The re- 
maining 15 runs were acquired during the spring of 1977 for winds 
of northwesterly direction. For some of these data records, data 
were acquired simultaneously with the cup-vane system. 

The hot-film system Is extremely delicate. Lightning and powerline 
fluctuations have often caused difficulties with the operation of the 
system. It is very seldcfci that the entire system is fully operational 
at one particular time. * The hot-fllms also have a tendency to undergo 
resistance shifts. If an appreciable shift Is detected, the probes 
are recalibrated. Corrections have to be made for changes in cable 
resistance due to changes in ambient temperature. Similarly, heat 
transfer corrections have to be made for changes In temperature. 

Because of the uncertainties in these corrections and other variations, 
the mean velocity and mean temperature measured with this system 
are not reliable. Resistance shifts in the active part of the system 
(films and cables) result in a parallel shift of the heat transfer 
(voltage) /velocity calibration curve. This shift of course will 
affect the mean quantities, a great deal but should not affect the 
calculated turbulence quantitltes as much. As pointed out in 
Reference 5, the results from these instruments become less accurate 
for wind directions Of +40° and +90° with respect to the probe 
axis. Consequently, the probes are rotated in the direction of the mean 
wind prior to data acquisition. For southerly winds this is no problem 
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since these winds especially at the higher elevations are very steady 
and have low turbulence levels. For northwest wlndt? the alignment 
with the mean wind Is more of a problem because of the presence of long- 
period fluctuations In direction. Precautions were made as much as 
possible to ensure that data of the highest quality were acquired, 
and It Is believed that the measured turbulence quantities fall within 
an accuracy level of less than 10%. 

5. PLANETARY BOUNDARY LAYER (PBL) 

The planetary boundary layer (PBL) may be defined as that part 
of the atmosphere where the effect of the earth's surface Is directly 
felt. The flow structure of the boundary layer Is extremely complex 
due to the variability of surface roughness changes In terrain, changes 
In surface temperature, variability of water vapor, presence of clouds 
and the fact that the flow is turbulent. Consequently, a simple model 
describing all the variables in the PBL such as velocity, wind direction, 
temperature and humidity, and covering all possible conditions is still 
not available. 

The unstable or convective PBL is characterized by a strong upward 
heat flux from the surface and by strong vertical mixing due to positive 
buoyancy forces/. Under these conditions above the surface layer, a 
well-mixed layer exists with an almost uniform potential temperature 
and an almost constant wind speed and direction. Due to the convective 
mixing a relatively sharp inversion is created on top of the mixed layer 
that delineates the depth of the PBL. Above this inversion the atmo- 
sphere is relatively undisturbed by the presence of the earth's surface 
and only gravity waves are generated. During the course of a sunny day 
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Che convective PBL Increases In depth as the land surface heats up and 
extends up Co the inversion layer, which is the result of the convective 
mixing of unstable air below with warmer stable air above. Under these 
conditions the convective PBL can be characterized by three different 
layers, which are, starting from the earth's surface: 

1) a surface layer, 2) a mixed layer, and 3) a capping inversion. The 
existence of this model for Che convective PBL is based on theoretical 
considerations [6), numerical flow modeling [7,8,9,10] , laboratory flow 
modeling [11] and actual Observations [12,13^14]. 

In the case a layer of stable air is present at the surface, and/or 

one or more stable layers exist at higher altitudes, no capping inversion 

exists. In stable air, turbulence is suppressed and may vanish completely 

under extremely stable conditions, and the air in the different layers 
becomes uncoupled as a result of reduced mechanical mixing. No simple 
model is available to describe this situation. A surface-based inversion 
with or without one or more stable layers at higher elevations has been 
observed over water surfaces, as warm air flows over the copier water 
during either day or night [3,18]. Under these conditions a low-level 
wind maximum (surface jet) is usually observed and the interpretation of 
the wind fluctuations near the surface is often complicated as a result 
of the co-existence of turbulence and internal gravity waves. These two 
phenomena have quite different properties although a nOh-llnear inter- 
action may exist between them. 

More complications are introduced in the flow analysis of the PBL 
when clouds are present and condensation of water vapor occurs within 
the PBL. As condensation takes place in the layer where the clouds are 
located, latent heat is released and an inversion layer develops below 


the cloud base that :1s responsible for some degree of suppression of 
the turbulence. Under these conditions the large-scale plume structure, 
normally present In the convective PBL, does not develop resulting In 
a considerable suppression and reduction In low-frequency velocity 
fluctuations. This phenomenon In turn results in reduced velocity 
variances and an Increased roll-off In the low-frequency range of the 
spectra of the horizontal velocity components. 

Just above the earth's surface the shear-production of turbulence 
dominates the buoyant production or suppression of turbulence. HOwever 
the Importance of the shear-produced turbulence diminishes with height 
as the effect of buoyancy on the turbulence gains Importance. In the 
layer just above the earth*'S surface where shear-produced turbulence 
dominates and the Reynolds stress Is nearly constant, the velocity Is 
adequately expressed by the well-known logarithmic law, provided the 
roughness of the terrain is uniform. As buoyant production or buoyant 
suppression of turbulence gains importance with height relative to shear 
production, modification of the logarithmic velocity profile will result. 
The buoyant production of turbulence Is associated with the upward heat 
flux of sensible and latent heat, which usi>'illy occurs when the atmosphere 
neat the earth's surface is unstably stratified ^ 0)* The buoyant 
suppression of turbulence is associated with the downward heat flux 
which usually occurs when the atmosphere is stably stratified (~ > 0), 
where 6 Is the local potential temperature. 

Under strong-wind conditions on a sunny day, the layer in which the 
shear-produced turbulence dominates Increases In height, and the logarith- 
mic velocity profile exists to higher elevations before it is eventually 
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modified as a result of convective activity, even under the strongest 
wind conditions. 

Following this discussion of the PBL, It is unlikely that Its 
flow characteristics are based exclusively on mechanically or shear- 
produced turbulence. Experimental results describing the mean flow 
and turbulence In the PBL, under a variety of stability conditions 
ranging from unstable via near-neutral to stable conditions, clearly 
Indicate that some basic differences exist between the PBL and the 
typical zero-pressure gradient wind-tunnel boundary layer. The most 
noticeable differences are observed In the evolving convective PBL 
on a sunny day with the mixed layer occupying about 90^ of the height 
and the surface layer consisting of the lower 10% of the PBL. In the 
mixed layer the velocity and wind direction are practically uniform 
[12,13]. Also, the downward entrainment of heat and momentum at the 
base of the capping Inversion Is a phenomenon that makes the assumption 
of vanishing turbulence at the edge of the boundary layer untenable 
[12,13,14]. 

The entrainment of heat and momentum of different magnitude and 
direction into the mixed layer is mainly responsible for the low- 
frequency fluctuations observed lu the horizontal velocity components 
down to near ground level [19,20]. These low-frequency fluctuations 
are not present under stable conditions, although as mentioned before, 
low-frequency fluctuations associated with gravity waves have been 
observed under stable conditions [15,18]. 

An abrupt discontinuity (increase) in low-frequency content has 
been observed In the spectra of the horizontal velocity components as 
the thermal stratification changes from stable to unstable [19,20]. 
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Similarly, the variances of the u and v velocity components and the 
reduced peak frequencies of these spectra increase abruptly as the 
stratification changes from -‘table to unstable. These abrupt 
changes have been observed in the horizontal velocity components 
u and V only, and do not occur in the vertical velocity component w. 
Experimental evidence of the abrupt increase in variance of the u 
velocity component and of the total velocity variance as the flow 
changes from slightly stable to slightly unstable Is clearly presented 
by Busch [21]. For slightly stable stratifications the ratio of the 
standard deviations of the vertical and longitudinal velocities 
o^/o^ « 0.36 and for slightly unstable stratifications this ratio 
is reduced to 0.25. : 

It has been observed that shortly before sunset the convective 
boundary layer disintegrates suddenly in a matter of minutes as a 
surface-based Inversion begins to develop [22,23]. During the 
night (or when there is a heavy cloud cover and water-vapor condensation 
in the PBL) the large-scale turbulence structure of the convective PBL 
is suppressed and no low-frequency components are present in the u and 
V velocity spectra near the earth's surface. 

Some of the experimental evidence, used previously In support of 
the flow description of the PBL, was obtained under relatively 
strong wind conditions. For example, during the Minnesota experiment 
[12] two data records (2A1 and 2A2) were acquired each with an hourly 
mean velocity of approximately 10 m/s at a height of 10 m. For these 
two data records the mean velocity Is practically uniform above 60 m, 
which is typical for the mixed layer. Also the standard deviations 


of the three velocity components do not change appreciably with elevation 
up to an elevation of 1220 m. Indicating that there Is little evidence 
of vanishing turbulence at the edge of the PBL. In addition, the 
horizontal turbulence stresses uw and vw exhibit minima near the 
surface and a considerable upward heat flux Is observed up to the top 
of the mixed layer where the heat flux changes direction due to entrain- 
ment of warm air at the base of the capping Inversion. The gradlent- 
Rlchardson numbers for these two records of Reference 12 are approximately 
-0.30 and -0.A6 respectively, which excludes these records from the near- 
neutral stability category. 

In another example, extremely strong winds with mean velocities up to 
28 m/s have been observed at 10 m heights over water [24]. The turbulence 
Intensities of the u and v components show an abrupt Increase for mean 
velocities In excess of 12 m/s from 3% to 15%, o^/U from 8% to 

11.5% and o^/U remains constant at 5%). There Is strong evidence that 
helical vortices, which form over the ocean under near-neutral conditions, 
are responsible for this Increase In turbulence Intensity of the 
horizontal velocity components. Also nocturnal jets near the earth's 
surface with a jet veiocity of approximately 20 m/s and a gradient 
Richardson number of 40.5 have been observed In Nebraska [16]. 

Based on these examples and other experimental evidence presented 
In the cited references, there Is no basis for the assumptions that 
under strong wind conditions the turbulence In the PBL is purely of 
mechanical origin and that the PBL flow-structure is similar to that 
of a zero-pressure gradient wind-tunnel boundary layer with vanishing 
turbulence at the free stream. 
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For the benefit of structural engineers and wind engineers who 
deal with a large number of proble'iiis requiring the knowledge of the 
mean flow and turbulence In the atmosphere, the atmospheric flow 
characteristics near the earth's surface under strong wind conditions 
have been reviewed and summarized In several review papers [References 25 
through 31] . The Input data to most of these review papers were obtained 
from many different sources and Include multi-level tower data and airplane 
data usually taken over horizontal terrain with near-uniform roughness. 

The results from these strong-wind experiments have formed the basis 
for empirical models of the wind structure near the earth's surface 
and of the atmospheric flows at higher elevations . In many of the 
review papers it Is assumed that under strong wind conditions the 
turbulence Is purely mechanically generated, and that buoyant produc- 
tion and suppression of the turbulence can be neglected. Under these 
conditions it Is then assumed that the atmosphere Is neutrally stable 
and the mean and turbulent flow behave similar to the flow in a 
zero-pressure gradient turbulent t mdary layer developed in a low-speed 
wind tunnel. 

The review papers generally Indicate tht the PBL under strong wind 
conditions is neutrally stable but may be modified as the result of 
thermal effects. However, it Is generally assumed that under strong 
wind conditions there 1® sufficient mixing that thermal effects can be 
completely Ignored. The neutrally stable PBL-flow model, which is 
generally assumed in the review papers, has never been experimentally 
verified at different sites for a large variety of strong wind con- 
ditions. Nevertheless, it Is widely assumed by many engineers and 
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scientists in the Clcld on wind engineering tluit this model provides 
an adequate description o£ the flow in the I’Bl under strong wind 
conditions. 

Upon reading of the review papers one does not always get a clear 
picture when neutral conditions exist. Deaves and Harris [30] do 
not give any conditions for the existence of neutrally stable flow 
in tite PBL other then the mean velocity has to be strong enough. 

Counlhan [26] considers only those data which specifically Indicate 
adiabatic or near adiabatic conditions, or which have wind speeds in 
excess of 5 m/s at a height of z »• 10 ra. ESBU [27,28,29] considers 
neutral stability to exist when the hourly mean wind at a height of 
10 m is greater than 10 m/s. Using the gradient Richardson number 
as a measure for the degree of thermal stability in the PBL, near- 
neutral conditions are suggested to exist by Teunissen [25] for |K1| <0^.03 
and by Panofsky [31] for |Ri|<0.01 for flow near the surface. 

Panofsky [31] also states that wind shear, and tViercfore mechanical 
production of turbulence, decreases rapidly with height and that the 
effect of buoyancy becomes progressively more important. Consequently 
thermal effects can no longer be neglected even under strong wind 
conditions for heights above approximately 50 m. 

Most of the observations at Wallops Island presented in this 
report have been made with the mean wind speed at the 76 m (250 ft) 
level between 10 m/s and 20 m/s and therefore can be considered to be- 
long to the strong-wind category. However the results of these ob- 
servations will be presented and discussed in the framework of non- 
neutral PBL- flow. 
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6. RESULTS ANO DISCUSSION 

This chapter deals with the discussion of the mean flow and 
turbulence measurements obtained with the two Instrumentation systems 
mounted on the 76 m (250 ft) meteorological tower at Wallops Island > 
Virginia. 

Data were acquired with the cup-vane and resistance temperature 
probes during a period of more than <i years (August 1974 - December 1978). 
Most of these data records were taken during daytime, although some 
night records are included. The prevailing wind directions at this 
sJi-o are southerly during the summer and fall, end between west and 
north during winter and spring. Consequently, for these two direction 
sectors many data records are available but the data base for on-shore 
winds in the sector between northwest and south is limited. For this 
set of data, mean velocity U, mean direction (j), mean temperature T, 
the turbulence variances and o^*, and the turbulence integral scales 

X z z 

^u* ^u* ^v obtained for each data record, based on the measure- 

ments from 5 levels at 15.3 m (50 ft) Intervals on the 76 m (250 ft) 
tower . 

A second set of data was taken with the hot-film instrument 
system. Nine data records, each one hour long and for southerly wind 
directions, were acquired during 3 days in July and August of 1976 [3j. 

In addition 14 data records for winds from northwesterly directions 
were acquired during several days in March, May and June of 1977 [23]. 


*The symbol 0 refers to the standard deviation, but for simplicity 
the word variance will be used in the text. 


For the.se data records the following parameters were obtained: variances 

and covariances of the three turbulence components and temperature, 
spectra and cospectra, and longitudinal Integral scales L^, and L^, 
all measured at the same levels as the cup-vane Instruments and at 9«1 m 
(30 ft). For some of these records data were acquired simultaneously 
with the cup-vane system, allowing direct comparison of the measure- 
ments with the two Instrument systems. 

6.1 Mean Wind and Mean Temperature Profiles 

For near-neutral stratification the wind profiles over homogeneous 
terrain obeys the relation 

U =$U*/K) ln(z/z^) (2) 

where U* Is the friction or shear velocity, which Is Ideally equal to 
y^/p, where Is the surface stress and p the air density, K Is 
Von Karman's constant taken at 0*4 and z^ Is the roughness length. 

The height z below which the log-law (2) Is valid depends on the 
stability of the flow or on the relative Importance of the buoyant 
production or suppression of the turbulence with respect to the mechani- 
cally produced turbulence. The height where the measured velocity 
profile starts to deviate from the logarithmic profile varies a great 
deal and the deviation Increases progressively for Increasing heights. 
Under strong convective conditions and under extreme stable conditions 
the near neutral part of the surface layer Is well below 15.2 m (50 ft), 
the height of the lowest mean-velocity measurement. Consequently under 
these conditions it Is Incorrect to fit the log-law (2) to the measured 
mean-velocity data in order to obtain values of z and U*. 
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The theoretical velocity profile over uniform terrain for the 
upper part of the surface layer where non-neutral stratification 
exists, Is based on the Monln-Obukhov similarity theory and Is given 
by [321 

U = (U*/U) [ln(z/z^-ipl (3) 


where t Ts ^ universal function of the stability parameter z/L and 
L Is the Monln-Obukhov length defined (In the absence of moisture) as 
[321 


K g Q„ 


(4) 


where T Is the absolute temperature, U* = and Is the surface 

heat flux, which can be approximated by c^p w6 measured In the surface 
layer. The stability parameter, z/L, depends on the gradient Richardson 
number [321 



Yg + dT/dz 
(dU/dz)^ 




(5) 


where y Is the adiabatic lapse rate (0.01“C/m for dry air and 
& 

0.0065®G/m for air saturated with moisture). With the use of the 

expressions relating Rl, z/L and ip as given by Panofsky [31], the 

departure from the logarithmic profile due to buoyancy effects 

In the surface layer can be obtained so that (3) can be used to 

obtain estimates of z and U*. • 

o 

This approach can be taken for velocity profiles measured In 
the surface layer, which Is loosely defined as that part of the PBL 
where the horizontal stress and vertical heat flux are nearly con- 
stant. The height of the surface layer Is frequently estimated as 


25 


the lower lOX of the convective PBL depth* z^> and is limited to a 
height z<L [13] . 

For stable conditions L varies a great deal and is relatively 

small [17]. The surface layer is not well defined and often is only a 

few meters high [33]. Under these conditions the velocity profile for 

non-neutral conditions cannot be used for estimation of z„ and U** with 

o o 

velocity measuremente taken well above the surface layer. 

In general, Monln-Obukhov (M-0) similarity theory can only be applied 
in the surface layer, when the flow conditions are stationary, when 
surface conditions are uniform (roughness and temperature), and over 
level terrain without major topographical features such as mountains 
etc. If the measured velocity data are linear with In z, then this 
observation does not automatically guarantee that the profile is truly 
logarithmic. For the majority of sites including the Wallops Island 
site the above conditions are not met, and application of the M-0 
similarity concepts such as the empirical flux-profile relationships, 
which are based on data with the above conditions satisfied, is 
questionable and should probably be avoided. 

For neutral conditions the temperature profile should follov; 
the adiabatic lapse rate of either 0.01 or 0. 0065*’C/m for dry and 
saturated air respectively. It is very seldom, even under the strongest 
wind conditions, that a truly neutral stratification is encountered 

for any length of time. The usual daytime thermal stratification near 

dT dT 

the surface is - -j— > y (unstable stratification) and - -y < y (stable 

u Z 21 u Z 3 

stratification) for nighttime. There is a short time around sunrise 
and in the late afternoon before sunset when a neutrally stable strati- 
fication is observed near the surface in the first 10 or 20 m. The 
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stratification of the air at higher elevations varies a great deal with 
time of day, cloud cover and season. For southerly flows over the ocean 
the air temperature Is usually higher than the water temperature and 
a surface-based Inversion (-dl/dz < Y ) exists day and night. 

Changes In temperature and velocity profile occur also as the flow 
experiences a sudden change In surface roughness and surface temperature, 
and an internal boundary layer (IBL) develops as the flow adjusts Itself 
to the new surface conditions. On-shore winds at Wallops Island usually 
experience an Increase In roughness and an Increase In surface temperature 
in daytime as they cross the beach. For westerly winds the surface 
temperature changes as the wind moves from the wet marsh over the warm 
land surface during daytime*. Within the observation height of 76m at the 
Wallops Island site, temperature gradients vary greatly with height and 
even change from positive to negative or vice versa. Under these condi- 
tions the local Richardson number is not a true Indicator of the overall 
thermal stratification of the observed flow and M-0 similarity does not 
apply. 

Typical temperature profiles for southerly winds on a summer afternoon 
are shown in Figure 5. These temperature profiles clearly Indicate the 
stable conditions above the IBL in the early afternoon, and the surface 
cooling in the late afternoon. Daytime and nighttime temperature: profiles 
for strong westerly winds are shown In Figure 6. During daytime the 
surface temperature is maximum and gradually decreases with height. On 
a rainy day with low cloud cover a stable stratification Is observed above 
40 m. Neutral temperature profiles throughout the tower height are seldom 
observed; the second temperature profile shown in Figure 6 is the closest 
to an observed neutral temperature' profile. Typical nighttime temperature 
profiles show a stable stratification near the surface with near neutral 
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conditions above 15 m. 

Ratios of mean velocities at the 15.2 m (50 ft), 50.5 m (100 ft), 

45.7 m (150 ft) and 61 m (200 ft) levels with respect to the mean velocity 
at the 76.2 m (250 ft) level are shown in Figure 7. The plotted data in 
this figure represent the mean and the plus and minus standard deviation 
from the mean, as well as the maximum and minimum values of the velocity 
ratios for the data records acquired in each 15-clegree sector. Because 
the mean velocity at higher elevations is less disturbed by local surface 
obstacles and is outside the developing IBL, the mean velocity at the 
76.2 m (250 ft) level was chosen as the reference velocity. For the one 
sector between 0° and 15“ and the 5 sectors between 90“ and l65“ (east- 
southeast) the number of data records available in each sector was four 
or less, an insufficent number for determination of the standard deviation, 
and only the average velocity ratios are presented. 

Between 230“ and 360“ (southwest-north) the variation in mean wind 
speed between the different elevations is relatively minor under moderately 
strong wind conditions. Conversely, for southerly winds in the sector 
160“ < 4i < 230, large variations in mean wind speeds relative to the 
mean velocity at the 76.2 m (250 ft) level are observed. For example 
in the sector between 180“ and 195“ the mean wind-speed ratio, V50/V250, 
varies from 0.4 to 0.8, clearly indicating the great variability of mean 
velocity distribution near the surface. Relatively large variations of the 
mean velocity ratios are also observed for wind directions approximately 
parallel to the island for sectors 30"-45“ and 210“-225“. For winds from 
these two sectors a small change in wind direction gives rise to a large 
variation in upstream terrain roughness. Generally the variations of 
the mean velocity distribution near the surface for on-shore winds are 
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much larger Chan £or westerly winds. For on-shore winds, Che changing 
development length of the XBL, the thermal stratification and other 
atmospheric conditions have a great effect on the velocity distribution 
near the surface as Che data indicate. For southerly winds, a surface- 
based inversion exists during the daytime as the result of warmer air flowing 
over Che colder water [3]. Under these conditions the existence of a sur- 
face jet with a maximum mean velocity (jet velocity) within a few hundred 
meters from the surface has been observed. The jet velocity has been 
observed as low as 45.7 m (150 ft) as shown by the maximum values of 
Che velocity ratio V150/V250 of larger Chan one for wind directions be- 
tween 170° and 210°, and by one of the velocity profiles of Figure 8. 

Typical strong-wind velocity profiles for wind directions between 
southwest and north are shown in Figure 9. The stronger wind profile 
is linear over the entire observation height when presented in 
semi-logarithmic coordinates. In the second profile with somewhat lower 
velocities, a distinct "kink" is observed. The majority of the measured 
velocity profiles for westerly wind directions shows similar "kinks". 

The values of obtained from the velocity profiles above the "kink" 
are much too small for the upstream terrain. Moreover the two velocity 
profiles of Figure 9 are approximately from the same direction. Therefore 
the "kinks" in the profiles cannot be the result of an upstream change 
in roughness, but must be interpreted as the beginning of the transition 
from the surface-layer flow to the uniform mixed-layer flow [13]. 

6.2 Velocity Profile Parameters z^ and a 

The roughness length, z^, can be evaluated from fitting of either 
the logarithmic law (2) or the non-neutral profile law (3) to the 
measured mean velocities in the surface layer. The height of the 
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surface layer In which (2) and (3) are valid is roughly defined as the 
layer near Che surface which Che fluxes uw and wd are approximately 
constant and not less chan 80% of the surface stress, A Vp» and surface 
heat flux, Q^/c^p, respectively. Measurement of the turbulent fluxes at 
Wallops Island indicate that Che thickness of the surface layer varies a 
great deal with wind speed, wind direction and thermal stratification. 

For on-shore winds the top of the surface layer is generally 
below the lowest observation level of 9.1 m (30 ft) and for the 
strongest westerly winds the surface layer extends above the height of 
the tower. For on-shore winds when an IBL develops as the air cresses 
Che beach, turbulent heat-flux measurements show an upward flux at the 
lower levels and a downward flux at the higher elevations [3]. Transi- 
tion usually occurs between Che 15.2 m (SO ft) and 30.5 m (100 ft) levels, 
which should correspond to the height of the IBL at the tower location. 
The flow above the IBL is still associated with the ocean surface, but is 
well above the Surface layer, which for southerly winds extends only a 
few meters above the water surface. Consequently no values of the 
roughness length, z^, for either the land or ocean surface can be ob- 
tained from the measured velocity profiles. Despite the fact that two 
velocity profiles of Figure 8 show a linear variation of velocity with 
Inz, fitting of the log-law (2) to the velocity data leads to values 
of of the order of Im, which are much too high for either the existing 
upstream terrain or the o.cean. The parameters describing the on-shore 
turbulent flow depend In part on Che roughness of the underlying surface 
(ocean) but also to a great extent on the stability of the flow. Con- 
sequently, knowledge of Che roughness length, z^, for flows over the 
ocean is not sufficient for the prediction of the turbulence 
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parameters for higher elevations as presented by ESDU [28»34]. 

2 —2 1/4 

Por westerly winds the local shear velocity U*-[uw + vw ] 

Is approximately constant over most of the observation height* In-' 
dlcatlng that the surface layer extends well above the lowest obser- 
vation level of 9.1 m (30 ft). Also the values of the local shear 
velocity, U*, compare quite well with the profile friction velocity 
U*, obtained from mean velocity measurements below the profile "kink". 
The roughness length, z^, obtained from fitting of the log-law (2) to 
the measured mean velocities below the profile "kink", vary with up- 
stream roughness and thermal stability of the flow In the surface layer. 

Figure 10 shows the variation of z^ with wind direction, clearly 
Indicating the effect of the 6.5 m high storage bunker for mean wind 
directions of 270°. In the sector between 350° and 30° the upstream 
terrain Is much rougher which results In higher values of z^. For the 
sectors between 230° and 260°, and 280° and 340°, the terrain Is 
reasonably uniform and the average value of the roughness length, 
z^=0. 0$4 m , corresponds quite well to the predicted values of 
ESDU [28,29] fcr similar terrain. In Figure 11 all values for z^, ob- 
tained for wind directions in the before-mentioned sectors, are shown 
as a function of the gradient Richardson number evaluated at z=<15.2 m 
(Rij^j) . The values of were obtained by fitting of the log- law to 
the velocity measurements below the profile "kink". The scatter of the 
data is appreciable for Rij^j>-0. 15, however the values of z^ decrease 
rapidly for Rij^^<-0.15. At this point the mean velocity 
profile In the surface layer can no longer be estimated with the log- 
law (2), as stability effects start to dominate. An attempt was made to 
estimate the roughness length by fitting the non-neutral profile (3) 
to the mean velocity data above the "kink" as outlined by Panofsky [31]. 
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However • this correction technique leads to unrealistic estimates of 
which seem to indicate that the velocity profile above the "kink" 
is not part of the surface-layer profile but Instead is the transition 
of the Surface layer profile to the mixed layer profile. 

In Table I mean profile parameters and turbulence parameters are 
presented for strong-wind data records with wind directions in the 
above-mentioned sectors with near-uniform, upstream terrain and measured 
with the cup-vane system and the hot-film anemometers. For the data 
records with a thermal stratification close to neutral the semllogarlth- 
mlc velocity profiles (S) are generally linear throughout, and for the 
records for which the flow is more unstable, kinks (K) start to appear 
in the profile. . 

Engineers generally favor the power law as an expression for the 
velocity profile through the entire PBL over uniform terrain 

U/Uj^ = (6) 

Although under certain conditions this law has fitted observed profiles 
over uniform terrain quite well, it is now generally accepted that the 
log-law (2) is preferable over the power law (6) for flow in the sur- 
face layer. Also the power law (6) cannot be expected to be a good 
approximation for wind profiles in the convective PBL. The latter seems 
to be prevalent over the North American continent even under the strongest 
wind conditions. 

The mean velocity profiles of records 2A1 and 2A2 acquired during 
the Minnesota 1973 atmospheric boundary layer experiment [12] as well as 
the velocity profile acquired at the Boulder Atmospheric Observatory 
(300 m mst [35]} on the morning of the 11th of September 1978 under 
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extremely strong wind conditions (ZO-minute average velocity of 18.67 m/s 
at z*10 m) show clearly that the power law with one constant exponent 
does not fit the velocity measurements. The measured strong-wind 
velocity profile at the Boulder site starts to deviate from the log-law 
(2) at z*22m and becomes more or less uniform above 200 m. 

On the other hand mean-velocity measurements made at 3 levels (10m, 

80 m and 200 ro) of the 213-m mast at Cabauw, the Netherlands, under 
extremely strong westerly winds (30-minute average velocity of 22.2 m/s 
at z *10 m) [36] indicate that the power law (6) fits these data quite 
well. The Cabauw tower is operated by the Royal Netherlands Meteorological 
Institute and is located about 50 km east from the North Sea coastline. 

For westerly winds the upstream terrain is flat low-lying pastureland with 
a very shallow water 'table. Simultaneous temperature measurements 
made at 8 levels between 2 m and 200 m Indicate a near- neutral thermal 
stratification with possibly a minimal upward heat flux from the mostly 
wet upstream terrain that inhibits the formation of a convective PBL, 
which is typical for the observations of the two North American sites. 

As Panofsky [31] has pointed out the power law can be expected 
to fit the velocity data in the surface layer only over a limited height 
range. In the near-neutral surface layer where the logarithmic law 
applies the power-law exponent can be approximated by [31] 

a= [ln(A^z^/z^) ]“^ (7) 

clearly Showing the variation of a with roughness:, z^, and geometric mean 
height, where z^ and Z 2 represent the elevation boundaries of 

the layer over which a near constant a may be expected. Near the sur- 
face a has a relative large value depending on the value of z^, but 
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decreases with height and should approach zero outside the surface 
layer as transition to the mixed layer takes place. 

For westerly winds at the Wallops Island site, power law exponents 
based on the mean velocities of the three lower levels (z^>15.2 m, 

m) are shown in Figure 12 as a function of the roughness length, 
z^. Similarly, power law exponents based on the mean velocity measure- 
ments under strong-wlnd conditions (no "kinks") at all five levels and 
the three highest levels are shown versus roughness length, z^, in 
Figures 13 and 14 respectively. Panofsky's relation (7) fits the 
measurements extremely well except for the results based on the mean 

velocities of the three highest levels (Figure 14). The reason for this 

#: ' 

is that on a semi-logarithmic plot of U vs. In z, the results of the 
higher levels fall close together and may seem to vary in a linear 
fashion. However in reality the data already deviate from the log-law (2) 
and all three prediction methods, which are based on known roughness 
lengths, overestimate the power law exponents obtained from the measured 
velocity profiles. No attempt was made to obtain power law exponents 
for the velocity profiles of on-shore winds, because of the variability 
of the profiles, the presence of an IBL and the absence of values for 

*o- 

6.3 Turbulence Intensities 

Average turbulence intensities of the horizontal velocity components 
u. and V measured with the cup-vane instruments at 5 levels, are shown in 
Figures IS and 16 for each IS-degree sector. In addition to the average 
turbulence intensities, the maximum and the minimum observed turbulence 
intensities for each sector and the mean intensity plus and minus the 
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standard deviation are shown. In the secto'i^ and the 5 sectors 

between 90* and 165* only average turbulence intensities are shown, 
because the number of data records available in each of these sectors 
is four or less. For on-shore winds o^/U is approximately lOZ at 
z*15.2 m decreasing to 5% or less at the highest observation level 
(z"76.2 m). For westerly winds the o^/U decreases from about 20% at 
the lowest level to about 14% at the highest level (Fig. 17). In 
general the variation of the turbulence intensity, o^/U, in each 
section is relatively small, except for the two sectors between 195* 
and 225*. 

The turbulence intensity of the v component, o /U, shows much the same 

. ^ 

pattern, although some differences can be observed. For on-shore winds 
the average value of o^/U at the lowest observation level (in the IBL) 
is about 8%, decreasing to about 4% at the highest level (above the IBL). 
For southerly winds in the three sectors between 165® and 210* turbulence 
intensities for both u and v components of less than 2% have been ob- 
served at the highest observation levels above the IBL. These observations 
have been made under stable conditions, and on several occasions the 
mechanical turbulence has been observed to vanish completely at these 
elevations at mean velocities of 10 m/s and higher. 

The average turbulence intensity o^/U is nearly constant at each 
level for wind directions between 240* and 345*. The average turbulence 
intensity o^/U at each level- increases with wind direction from 240* to 
about 300* where a maximum is reached and decreases gradually with wind 
directions from 300° to about 50*. Figure 17 shows the variation of the 
average turbulence intensities of the u and v components with height 
in the two sectors 240*-255° and 300*-315*. The averaged results are 
also compared with the estimates of Teunissen [25] and ESDI) [28] both 
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based on Che roughness length, z^, whose average is 0.037 tn for either 
one of these sectors [Fig. 10]. Although the upstream terrain for the 
two sectors Is about Identical » the values of o^/U In the sector be* 
tween 300* and 315* are considerably higher than chose for the sector 
between 240* and 255*. The average values of a^/U for the latter sector 
correspond extremely well with Che Teunlssen estimate, but for the 
sector between 300* and 315* Che average values of o^/U are about 3 
to 4% higher. On Che other hand the average values of 0^/0 for these 
CWQ sectors are about Identical and fall between the estimates of 
Teunlssen [25] and ESDU [28]. 

A possible explanation for this unusual behavior of o^/U can be 
derived from the visual Inspection of several years of strip-chart re- 
cordings of wind speed and direction obtained contlnously from a 
propeller-vane anemometer located about 1 km northeast from the Cower 
on Wallops Island. These recordings clearly show that for northwesterly 
winds Che instantaneous wind direction experiences frequently large 
fluctuations Coward Che north (Fig. 18), These direction fluctuations 
are larger than usual and are not normally distributed as shown In 
(Fig. 19), which explains the large values of both o^/U and the co- 
variance vw measured In this sector. 

The turbulence intensity, o^/U, which is constant over near-uniform 
terrain for westerly wind directions (Fig. 17), Instead varies with time 
of day. Nighttime measurements are systematically 1% higher than the 
daytime results (Fig. 20) . The vertical distribution of the vertical 
turbulence intensity, o^/U, based on five strong-wind data records ob- 
tained with Che hoC-f Urn system , is shown in comparison with the ESDI! 
[28] and Teunlssen estimates in Figure 21. 
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The vertical distributions oi: the turbulence Intensity of all 
three components for southerly winds are shown in Figures 22 » 23 and 
2A. For these wind directions an internal boundary layer (IBL) develops 
as the air crosses the beach. Based on the change in direction of the 
daytime vertical heat flux w9, which is expected to be positive in the 
IBL and negative above it [3J, the height of the IBL at the tower must be 
between the 15.2 m (50 ft) level and 30.5 m (100 ft) level. This 
observation agrees with the relationship given by Elliot [37], which 
predicts a height of the IBL of approximately 26 m, based on an upwind 
roughness length over the ocean of *=0.001 m and a downwind roughness 
length of z^"=0.01 m and a development length of 300 m. The average 
turbulence intensities p£ the u and v components from nine data records, 
obtained with the hot-film Instrumentation during summer afternoons, com- 
pare surprisingly well with the average turbulence intensities from 29 
data records, obtained during winter and early spring under strong wind 
conditions with the cup~vane system (Figs. 22,23,24). The turbulence 
Intensities in the IBL compare reasonably well with the ESDI! [28] pre- 
dictions based on a roughness length of z^=0.01 m. However above the 
IBL the ESDI! [28] predicted values of the turbulence intensities, based 
on a roughness length z^=0.00l m, overestimate the actual measured values 
by a factor of two. The turbulence intensites for run 7 [3] obtained 
under extremely stable conditions (R£=+25.5, U=10.7 m/s at z=61 m) , which 
are included in Figures 22 through 24, show lower than average values above 
the IBL. Under similar conditions it has been observed that the turbulence 
vanishes completely for some time. The variation of the ratios of 

average turbulence intensities, o /a and o fa , for all wind directions 

V u w u 

is shown in Figure 25. The measurements indicate that the values of 
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these ratios are more or less independent o£ height, and the data 
shown In Figure 25 are not only the ratios of the average turbulence 
intensities of all data records in each sector but are also averaged 
over all observation levels. For most wind directions the values of 
o^/o^ are between 0.75 and 0.85 in comparison with the predictions 
of 0.75 and 0.80 by Counlhan [26] and Teunlssen [25]. The estifflates 
of ESDU [28] are not single valued but are dependent on height and rough- 
ness length. Smaller values than 0.75 for are observed for wind 

directions approximately parallel to the Island at AsAO" and <^-230**. 

Values of 0^1 a ^ higher than 0.85 are observed for southerly winds for 
4>“170“ and for northwesterly winds in the sector 280“<4><350“ . That 
the values of are relatively high in this last sector is no 

big surprise since often large direction fluctuations have baen 
observed la this sector (Fig. 18). Values of and obtained 

from the hot-film data records are also shown in Figure 25. The values 
of o^/cr ^ obtained with this system compare quite well with the cup-vane 

results. The values of the ratio, a /a , for winds from the sector 

w u 

between south and southwest and northwesterly winds fall between 0.55 
and 0.60 as compared to values of 0.50 and 0.52 as predicted by 
Counlhan [26] and Teunissen [25] respectively. 

The results discussed so far in this report show clearly that at 
the Wallops site large variations in mean as well as turbulent flow 
occur varying with wind direction. The Important observed deviations 
from the simple neutral boundary-layer models are: 

1. The development of an internal boundary layer (IBL) for 
on-shore winds as they cross the beach. 

2. The existence of a surface Jet for southerly winds with 
extremely low turbulence intensities (2% or less and occasionally 
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vaivlshlng under stable conditions) coexisting with gravity 
waves . 

3. Large direction fluctuations towards the north for north- 
westerly winds* observed during the daytime, which are 
responsible for higher than usual lateral turbulence intensities. 

These observations were made under strong wind conditions and there is 

no reason to believe that for still stronger wind speeds (maximum observed 

speeds) these deviations from the neutrally stable boundary-layer model 

will suddenly vanish. 

6.4 Turbulence Integral Scales 

X z z 

6.4.1 Integral Scales, L^, and from cup-vane data 

In this section the distribution of the turbulence integral scales 
obtained from measurements ^with the cup-vane system will 
be discussed and compared ‘with the estimates from several review papers 
[25,26,28,29]. The streamwise turbulence Integral scale of the u 
component, L^, is calculated from the autocorrelation function (1), 

R^(t), assuming that Tayloc's hypothesis is valid. Averages of all 
scales, L^, obtained from the data records in each 15-degree sector, 
are plotted along with the maximum and minimum and plus and minus the 
standard deviation from the average value in Figure 26 for each ob- 
servation level. Averages only are plotted in the sector 0®-15“ and 
between 90° and 160** because of the limited number of available data 
records in these sectors. The Integral scales, L^, increase 
systematically with height and show a great deal of variabilJ ty for all 
wind directions. In general the magnitude of these integral scales 
is larger for westerly winds than for on-shore winds except for southerly 
wind directions. At the highest observation level two distinct extremes 
for the maxinui can be observed in the two sectors between 180° and 210“, 
and between 300° and 330°. For souchetly winds between 180° and 210° 
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the air Is frequently stably stratified which under certain conditions 
stay lead to the coexistence of turbulence and low-frequency gravity 
waves. If gravity waves are present the Integral scale obtained from 
the autocorrelation function can be expected to be high [3]t ^ maximum 
value of 800 m at the 76.2 m (250 ft) level has been observed. On the 
other hand, if the low-frequency gravity waves are absent and only tur- 
bulence of less than 5% Intensity Is present (Fig. 15), minlmiim turbu- 
lence Integral scales of less than 100 m have been observed for the 
southerly winds. In the sector between 300" and 330" large variations 
In are also present, which are the result of either the presence 
or absence of low-frequency velocity fluctuations as can be observed 
from measured u-spectra. 

Vertical Integral scales of Che horizontal velocity components 
obtained from the measurements with the cup-vane system can be cal- 
culated by Integration of the vertical correlation coefficients of either 
the u or V components. 


^1 “ f Rii(*')d2’ 

O 


where 


( 8 ) 


Rji(z*) 




i I 


u (t,z) u.(t,z+z*)dt 


(9) 


and 1 ■ o r v 

z 

Rj^^(z') Is the vertical correlation coefficient of either the u or v 
velocity fluctuations measdured at two different levels separated by 
a distance z'. For this research program with the Wallops Island 
tower the separation distance, z', can be either 0, 15.2 m (50 ft), 
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30.5 m (100 ft)» 45.7 m (150 ft) and 61 m (200 ft). The integration 


(8) should be performed to the point where the correlation coefficient 
z 

Ruu(z’) changes for the first time from positive to negative. However 
only a maximum of five values of the correlation coefficients are 
available for either upward or downward integration according to ex- 
pression 8, and often the correlation coefficient has still a large 
positive value for the maximum separation distance z'. 

In order to arrive at a reasonable estimate of the vertical scales, 
it is assumed that the vertical correlation coefficients of the hori- 
zontal velocity Components u and v decay exponentially in the same 
manner as has been observed by Dryden et al. [38] for high Reynolds- 
number grid turbulence, according to 

R^^(z') = exp[-z'/L^] • (10) 

Vertical integral scales of either the u or v velocity components can 
then be obtained from a least-squares fit of (10) to the available 
measured correlation coefficients. 

Because of the non-symmetric flow in the boundary layer, integral 
scales obtained from upward and downward integration of the correlation 
coefficients with the origin at a common point cannot be expected 
to be the same. Instead a slightly different definition for the 
vertical Integral scale is used as suggested by ESDU [29], where the 
correlation coefficient is. defined as 

R^i(z*) = - Y I v^(t,z+z*)v^(t,z-z')dt (11) 

(a^(z)}2 Jq 

where 1 - u or v. 
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With the use of this definition only one vertical Integral scale Is 
defined for each height z. However this definition can be used for 
the evaluation of the Integral scales of the u or v component at the 
45.7 m (150 ft) level only. 

The distribution of the vertical turbulence Integral scales of 
the horizontal velocity fluctuations and obtained with the cup-vane 
system Is shown In Figure 27 and 28 respectively. Calculation of 
these scales Is based on expression (10) . Integral scales obtained 
from upward integration of the correlation coefficients are shown 
as LUZf or LVZ+ on the figures and as or in the text. For 

the Integral scales obtained from downward integration the direction 
of the arrow is reversed.' The two-sided integral scales obtained 
according to the definition (11) are shown as LUz| or LVz| on the 
figures and as L^j or L^l in the text. Averages of all the Integral 
scales for each sector as well as the maximum and minimum, and the 
mean plus and minus the standard deviations for each each sector are 
shown in these figures. In the sectors 0“-15‘’ and between 90“ and 
165“ not enough data records were available to calculate a standard 
deviation and only average values are plotted. 

■ ■ Z ' 

The integral scales Ly for on-shore winds are generally smaller 

than those for westerly winds by less than a factor of two. The 

• ' . ' z 

integral scale measured at the 15.2 m (50 ft) level and the scale 

z ■ 

measured at the 76.2 (250 ft) level are of the same magnitude. 

The smallest observed value of is about 10 m for southerly winds 
without the presence of gravity waves. The larger values of are 
also observed for southerly winds when gravity waves are present. 
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Large values of are also observed for northwesterly wind directions 

between 300* and 320*. The upward integrated scales, , obtained 

from the lower three levels for wind directions between 0* and 100* 

show much less variation than the downward Integrated scales, , 

from the upper three levels in the same sector. Comparison of the 

three different scales, L^f , L^'t- and L^I of the measurements at 

u u u ▼ 

the 45.7 m (150 ft) level shows that for winds between 0* and 100* 
the upward and downward Integrated scales are of the same magnitude 
(37 m) with the scale having an average value of approximately 

43 m. For westerly winds between 250* and 360*, Is about 20 m 

larger than L^'t’ . 

The vertical Integral scale of the v-component, L^, behaves In a 

similar fashion as the scale L^, but the latter Is generally 

twice as large. The magnitude of for winds from westerly directions 

jr 

is about twice the value of for on-shore winds.. Minimum values of 

Just a few meters are observed for on-shore winds. Maximum values of 

these Integral scales are associated with winds from northwesterly 

directions (4>-310*) for which large direction fluctuations have been 

observed (Fig. 18). For southerly winds ((}i=180“) no large maximum 
z z 

values for L have been observed as for L . For on-shore winds the 
V u 

magnitude of the three different scales at the 45.7 m (150 ft) level 

are about the same but for westerly winds L 't' > L 't' with the value 

V v 


of L* I in between. 

Figures 29 and 30 show the variation of the relative magnitude 
of the average turbulence Integral scales L*/L^ with wind direction 
for each observation level. For southerly winds between 135“ and 225“, 
the ratio Increases with height from about 3.7 at the lowest 
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level Co values between 10 and 15 ac Che three highest observation 

levels. For all wind directions outside this sector the ratio Ly/by 

is 3.3 and 3.7 at the 15.3 m (50 ft) level and the 76.2 m (250 ft) 

level respectively. The magnitude of Che ratio of average integral 
z z 

scales, is 1»9» generally independent of height and wind 

direction. 

6.4.2 Comparison of the Cup-Vane Scales with Predicted Values 
from References 25, 26 and 29 

From the above discussion it is clear that Che measured turbu- 
X z z 

lence integral scales and vary significantly .with terrain 

roughness, wind direction, time of day, thermal stability and other 
atmospheric conditions. In this section the measured integral scales 
are compared with the predicted scales from either Teunisseh [25], 
Counihan [26] and ESDU [29J. 

In Figure 31 the averaged turbulence Integral scales, L^, ob- 
tained from the cup-vane data are shown for two westerly wind- 
direction sectors over near-identical terrain and are compared with 
the estimates of References 25, 26 and 29. At the higher elevations 
the values of in the sector 300‘’<iJ><315° are approximately 100 m 
larger than those in the sector 240“<iti<255® . The Counihan [26] pre- 
dictions match the measurements below 30 m but for higher elevations 
all three references under-estimate the measured integral scales, 

‘ 2 Z 

Averaged values of L and L for the same two wind direction sectors 
are shown in Figures 32 and 33 and compared with predicted values from 
Counihan [26] and ESDU [29], The results show the decrease in magni- 
tude when the change is made from upward to downward integration. 
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Values of in the sector 300‘*<<t><315* are about twice as large as those 

in the sector 240"<i^<255"» which is in contrast with the values of 

which are approximately identical in the two wind-direction sectors. 

The ESDU [29] predictions of both and $iatch the measured data 

in the wind sector 240^<4i<255" reasonably well, hut underestimate 

the scales for wind directions between 300” and 315” fur which large. 

wind-direction fluctuations have been observed [Fig. 18]. 

Figures 34, 35 and 36 show the averaged values of and 

obtained from the cup-vane data, in comparison with the predicted 

values for southerly winds. Since no measured values of roughness 

length are available for this wind direction, values of were 

selected in accordance with the nature Of the upstream terrain and 

Table 1 of Reference 28. Below 20 m, for the flow in the IBL, a 

roughness length z^=0. 01 m is appropriate and for the flow above 

the IBL a roughness length of z^=0.00I m was selected. The measured 

longitudinal scales, L*, are generally larger than the ESDU [29] 

z z 

predicted values, while the measured vertical scales and are 

generally smaller than the ESDU [29] predictions. These results support 

the likelihood that internal gravity waves exist in the surface-based 

Inversion over the ocean. The extent of the waves is much longer 

in the direction of the flow than in the vertical direction because 

of the suppression of vertical velocity fluctuations by buoyancy 

forces in stable air. The results of Figure 30 show clearly the 

X 

large extent of the streamwise integral scale, L^, relative to the 

vertical scale, for southerly winds only. For these winds the 

scale ratio L^/L^ above the IBL varies between 10 and 15, for all 
u u 
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other wind directions the magnitude of this ratio is of the order 
of 4. 

Longitudinal integral scales, obtained for westerly winds 
over near^unlform terrain arc dependent on time of day (Fig. 37). 

The scales acquired during nighttime are approximately 50% longer 
than those obtained from morning records with the scales from 
afternoon data falling in between. The increase of the scales 
during daytime can be explained with the Increase of the height 
of the mixed layer (See Fig. 5, Reference 13). However, no definite 
explanation for the large Integral scales observed during nightime 
is available. 

In Figures 38 and 39* individual values df the turbulence integral 
scales, L^, at the 15.2 m (50 ft) and the 45.7 m (150 ft) levels 
are plotted versus roughness length, z^. These results are for 
near-neutral strong wind data records for wind directions varying 
between 41=^250® and 4)"-30®. The scatter of the data is appreciable 
as is to be expected since the previously discussed results also 
Indicate variation of with wind direction (Fig. 31) and time of 
day (Fig. 37). The ESDU [29] and Teunlssen [25] predictions are 
consistently lower than the measurements, while the Counlhan predic- 
tions [26] fit the measured integral scales reasonably well. 

Similarly, the corresponding vertical scales L^| and L^| are shown 
versus roughness: length, z^, in Figures 40 and 41. The ESDU [29] 
predictions for match the measured results quite well, while the 
predictions for L^| of the same source fall generally below the 
measured data* except in the range 0.1 m<z^<1.0 m. 
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6. A. 3 The Direct and Spectral Methods for Obtaining Turbulence 
Integral Scales 

Integral scales of turbulence are defined for any correlation 
coefficient as the Integral over the entire range of the independent 
variable, which can be either time or space as previously discussed. 

In practice the Integration process is carried out between the origin 
and the first zero-crossing. Time scales are related to the length 
scales in the direction of the flow by assuming the existence of 
Taylor's frozen turbulence hypothesis. The magnitude of the time scales 
varies significantly depending on the presence of low-frequency fluctu- 
ations and or trends. In order to omit these fluctuations from the 
tlffle-correlatlon coefficients, the data should be. high-pass filtered 
at some low frequency. The filtering of the cup-vane data takes place 
as a result of the block averaging, excluding trends and low-frequency 
fluctuations below 0.00195 Hz from the sample. Similarly, the hot-film 
data are high-pass filtered at either Q.02AA Hz or at 0.00153 Hz de- 
pending on whether the data are analyzed in the high-frequency range 
(0.02AA-100 Hz) or middle-frequency range (0.00153-6.25 Hz) (See section 
A. 2). Turbulence Integral scales calculated in this way are said to 
be obtained by the direct method . 

An estimate of the size of the energy-containing eddies can also 
be obtained from the Von Karman interpolation formula for the three- 
dimensional power spectrum covering the wavenumber-range from the 
energy containing eddies to the inertial subrange [39]. 


S(k/k^) 


C(k/kg)^ 
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where corresponds to wavenumbers Ini the range of the energy-containing 
eddies. With the assumption of Isotropic turbulence* expressions for 
the physical realizable one-dlmenslonal spectrum functions S^(k^)* 

S (k.) and S (k|) can be obtained (see expressions 3-72 and 3-73* 

Hinze [40]). The wavenumber of the energy-containing eddies* k^* can 
be replaced by the Inverse of the turbulence Integral scales* L^* In the 
appropriate spectrum functions and the constants can be adjusted to fit 
the measured spectra* The one-dlmenslonal Von Karman spectrum functions 
obtained In this manner and presented by Teiinlssen [25] fit the measured 
spectra of mechanically produced wind tunnel turbulence quite well 
[41*42]. The spectral expressions for the streamwlse and lateral velocity 
components are given by. 


nS (n) 
u 


4(nLj/U) 


{l+70.7(nL*/U)^l^^^ 


(13) 


and 


nS.(n) nL 

- 4 - - 


l+188.4(2nL^/U)^ 

{l+70.7(2nL^/U)^}^^^^ 


(14) 


where 1 ■ v and w. 

For the comparison of these wind tunnel spectra with the normalized 

X 

spectrum functions* the turbulence Integral, scales* L^* were determined 

Independently using the previously discussed direct method. The 

-5/3 

Von Karmon spectral equations show the correct n -dependence In 


the inertial subrange. Integration of these expressions over the entire 
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frequency range leads to unity. At vanishing wavenumbers these 
spectral equations also predict the proper integral scale L^>S|^(0)U/4o^. 

At low frequencies the' logarithmic spectra vary as n*^^, which also fits 
the wind tunnel spectra quite well. 

However » caution must be taken with the automatic adaption of the 
Von Karman spectral functions to atmospheric turbulence under near 
neutral thermal stratification. The work of Kalmal et al [20] clearly 
points out that under near-neutral conditions the low-frequency content 
in the spectra of the horizontal velocity components can vary appreciably. 

In the convective boundary layer spectral data approach a definite shape 
when near neutral conditions are approached from the stable regime. 

Under these conditions significant low-frequency velocity fluctuations in 
the u and v components are absent and the spectral shape of these components 
is similar to the Von Karman model. On the other hand no unique spectral 
shape for the u and v component exists for near neutral conditions 
approached from the unstable regime. Under these conditions, the low- 
frequency spectral content is much higher and the Von Karman model does 
not represent the spectral data well. 

In order to check the validity of the Von Karman spectral 
functions, it will be necessary to check these relations against measured 
spectra, where the frequency is normalized with the local mean velocity 
and the turbulence Integral scale obtained independently via the direct 
method. In Figures 42 and 43 logarithmic Spectra of the u-coraponent 
measured at 5 different elevations for data record 19 (Table 1) are 
compared with the Von Karman spectral function. The difference in those 
two Illustrations is that variances and integral scales used in Figure 42 


were obtained In the mlddle-frequcncy range (0.00X53-6.25 Hz) while 
for Figure A3 the varlancca and Integral scales were obtained In the 
high-frequency range (0.02AA-100 Hz) filtering out the low-frequency 
fluctuations. In the latter case the Von Karmun expression fits the 
measured spectra quite well in the -2/3 region but In the low-frcqucncy 
range the measured spectral values arc considerably higher than the 
predicted values. Consequently no distinct spectral peak at the 
predicted value of the reduced frequency of nL^/U"0.lA6 is present, 
similarly* the measured v spectra follow the Von Karman spectrum in the 
-2/3 region only if variances and integral scales arc used from which the 
low-frequency components arc filtered out. Figures AA and A5 show 
the w-spectra for run 1119 in comparison with the Von Karman spectrum 
function for variances mtd integral scales obtained in the middle and 
high frequency range respectively. In the latter case the measured 
spectra fit the theoretical Von Karman spectrum function much better 
especially for the spectra from the higher elevations. 

Based on these results the conclusion can be drawn that the tlieoreti- 
cal Von Karman spectrum functions do not represent the spectra of atmo- 
spheric turbulence in the low-frequency range or in the high-frequency 
range when turbulence integral scales are used that are obtained via the 
direct method v hen low-frequency fluctuations are included. Better fit 
of the measured spectra in the -2/3 region is achieved when variances 
and turbulence Integt'^l scales are used that are obtained from data 
records from which the low-frequency fluctuations have been removed. 
Inversely* if the theoretical Von Karman spectrum functions are used 
to obtain stteamwise turbulence integral scales* 
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as Is suggested in references 23 and 28, Chen these scales are 
not equivalent to the integral scales obtained via the direct method. 

The Integral scales obtained via Che Von Karrnan method must be Inter- 
preted as integral scales associated with velocity records from which 
all low-frequency fluctuations with periods longer than approximately 
40 seconds have been filtered out. 

The expressions for the streamwise turbulence Integral scales, 
for l*u,v and w, listed In references 25 and 29 are based on scales 
obtained via the Von Karman method by either matching of the measured 
spectra at the peak reduced frequency or by u^lng a best overall fit 

of the spectra. Consequently the predicted Integral scales from these 

• * 
sources must be Interpreted as integral scales associated with the high- 

frequency content of the turbulence components. 

6.4.4 Comparison of liiitegral Scales Obtained Via the Direct and 

Spectrum Methods With Predicted Values from References 25, 26 
and 29. 

The turbulence Integral scales, L^, obtained from the cup-vane data 
are obtained via the direct method In the frequency range from 0.00195 Hz 
to 0.25 Hz. Integral scales obtained from the hot-fllm data and dis- 
cussed In this section, are obtained by one of the following three 
methods: 

1. The direct method In the middle-frequency (MF) range 
(0.00153-6.25 Hz). 

2. The Von Karman method, with the spectra and variances obtained 
from filtered data In the middle-frequency range (VK-MF) . 

3. The Von Karman method, with the spectra and variances obtained 
from filtered data in the high-frequency range (0.0244-100 Hz) 
(VK-HF) . 
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For method 2 and 3 values of were obtained by matching the measured 
logarithmic spectra (nS^(n)/a^ versus nz/U) at nz/U>10 to the Von Karman 
spectrum functions in the -2/3 range. Averaged and single-record stream- 
vise integral scales, L^, obtained from the cup-vane data or from the hot- 
film data, the latter derived through one of the above three methods, are 
shown for the two basic wind directions (south and northwest) irt Figures 
46 through 57 and are compared with the predicted values of references 
25, 26 and 29. 

Figures 46, 47 and 48 show the variation of the averaged Integral 

scales, with l*‘u,v and w, obtained from the hot-film data (methods 2 

and 3) versus height for southerly winds. The predicted values are based 

on a roughness leagth, a^'^0.01 m In the IBL below 20 m and on a roughness 

length z =0.001 m above the IBL. The values of L* from the cup-vane 
^ o u 

data match the ESDU [29] predictions and the scales obtained from the 
hot-film data via the Von Karman method in the middle-frequency range 
match the Teunlssen [25] predictions. Values of and obtained 
with the use of the Von Karman method in either frequency range fall 
well below the predictions. Of course, as previously discussed, the 
ratio of longitudinal scales and lateral scales is much higher for 
Southerly winds than for the other wind directions (Figs. 29, 30). 

The effect of buoyancy in the inversion layer tends to suppress the 
vertical motion and consequently the measured scales are much smaller 
than the predicted scales. It must be assumed that the predicted scales 
are based on data records taken under conditions where buoyancy had 

X X 

very little effect on the turbulence. The lateral scales, L, and L 
obtained via the Von Karman method in the middle-frequency range are 
about twice as long as the scales obtained from the same data in the 
high-frequency range. 
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Figures 49, 50 and 51 show rhe averaged Integral scales* for 

l*u* V and w* obtained from the hot-flln data and analyzed according 

to methods 1 and 2 for northwesterly wind directions. In Figure 49 

comparison of is also made with the cup-vane data. The scales 

obtained via the direct method from the cup-vane data are generally 

smaller than the scales from the hot-fllm data via the direct method 

in the middle-frequency range. The scales obtained with the Von Karman 

method In the middle-frequency range are systematically smaller than the 

scales obtained with the direct method from either the cup-vane or 

hot-fllm data In the same frequency range. The Teunlssen [25] and 

ESDI! [29] predictions only fit the measured scales obtained via the 

« 

Von Karman method below a height of 20 m* Above this height the pre- 
dicted values are systematically lower than the measured values* The 
Counlhan [26] predicted scales of fit the observed scales obtained 
via the direct method below 20 m, however above the height of 20 m the 
Counlhan prediction falls between the measured scales obtained via the 
direct method and those obtained via the von Karman method. 

The measured lateral scales L and L obtained via the direct method 

w V 

In the middle-frequency range are three to four times as large as those 
obtained via the Von Karman method In the same frequency range (Figs. 50* 
51). The latter scales match the Teunlssen [25] and ESDU [29] predictions 
very well. 

In Figures 52, 53 the Integral scales* L* , obtained from a single 
early-evening record are compared with the predicted scales and with 
the average scales from the daytime records all analyzed in the 
middle-frequency range. It has been observed [13* 23] that just before 
sunset the convective boundary layer dissolves abruptly and the low- 
frequency velocity fluctuations normally associated with the convective 
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PBL suddenly disappear as can be clearly seen from the comparison of 
the plotted records of the horizontal velocity components of record 16 
(evening run) and record 19 (afternoon run) (Fig. 54). The results 
indicate much smaller scales if the low-frequency fluctuations of the 
horizontal components are absent. For record 16 the scales obtained via 
the direct method are systematically larger than those obtained via the 
Von Karman method, although the difference is less significant for this 
early evening run than for the daytime records, which contain low- 
frequency fluctuations (Figs. 49-51) . 

Figures 55 through 58 show the averaged integral scales, and 
L*, from the daytime records and those for record 16 obtained via the 
Von Karman method in qlther the high-frequency range or the middle- 
frequency range. The Von Karman method for obtaining Integral scales 
provides near-identical results Independent of the frequency range if 
large low-frequency fluctuations are absent (e.g. record 16) or if the 
low-frequency fluctuations are removed from the data records. If the 
large low-frequency fluctuations are present (e.g. u and v components of 
daytime records Fig. 54) and are not removed from the daytime records 
the Von Karman method leads to much larger integral scales. If no 
appreciable low-frequency velocity fluctuations are present in the data 
records, the Integral scales obtained via either the direct or the Von 
Karman method are about the same in magnitude. However, if low-frequency 
components are present the magnitude of the scales depends on the method 
by which they were calculated and also depends on the frequency range in 
which the data are analyzed. 
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When the Integral scales, L^, obtained via the two- methods In 
either one of the frequency ranges, for data records with or without 
low-frequency content are compared with predicted values (Figs. 49, 52 
and 53) one observes appreciable variation. In general, the Teunlssen 

[25] predicted values correspond to scales obtained from data records 
for which the large low-frequency fluctuations are absent dr are filtered 
out,; and obtained via either method. On the other hand the Counihan 

[26] predicted values correspond more to the scales obtained from data 
records with low-frequency fluctuations and obtained via the Von Karman 
method In the middle-frequency range and to the scales obtained from 
the cup-vane data. The .magnitude of the scales obtained from the 
same data records via the direct method are generally larger than 
those predicted by Counihan [26] . The ESDI! [29] predicted 

scales fall between the two previously mentioned predictions. The scales 
obtained from data records of north-west winds with low-frequency 
fluctuations vary almost linearly with height, while the predicted scales 
show more a tendency toward independency with height at higher elevations. 

Similarly the lateral integral scales, (Figs. 50, 55 and 57) also 
show appreciable variation. The Teunlssen [25] and ESDU [29] predictions 
correspond reasonably well to the scales obtained during daytime via 
the Von Karman method in the middle-frequency range. The scales obtained 
from the same- daytime records via the direct method are significantly 
larger than the predictions. On the other hand If low-frequency compo- 
nents are absent or filtered from the data records the measured values 
for the lateral Integral scales, L^, fall well below the predicted values. 

The Integral scales, L^, (Figs. 51, 57 and 58) show a similar 
pattern, the ESDI! [29] and Teunlssen [25] predictions correspond best 
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to the scales obtained via the Von Kansan method In the middle 
frequency range. The scales obtained via the direct method In the 
same frequency range are considerably larger than the predicted values. 

If low-frequency components are absent or removed from the data records 
the measured scales are generally lower than the E3DU [29] and Teunlssen 
[25] predictions. 

In general It can be concluded that the ESDU [29] and the Teunlssen 
[25] predicted and scales should be Interpreted as scales obtained 
via the Von Kansan method for turbulence with low-frequency fluctuations. 

On the other hand the Teunlssen [25] predicted scales should be Inter- 
preted as scales obtained^ from data records from which the large low- 
frequency fluctuations are absent or are filtered out. In general the 
horizontal scales and obtained from the daytime data records via 
the Von Kansan method In the middle-frequency range vary linearly with 
height, while the predicted scales show a tendency of Independnece with 
height at the higher elevations. 

6.5 Power Spectra 

In this section the power spectra of the three velocity components 
obtained from the hot-fllm data are discussed for the two basic wind 
directions, south and northwest. Ample discussion of the spectra In the 
previous section has Indicated that the spectral shape In the low-frequency 
range depends greatly on the absence or presence of large low-frequency 
velocity fluctuations. 

The logarithmic spectra obtained from south-wind records generally 

-2/3 

show very little variation In shape and vary as f In the high-frequency 
range and approximately as f'*'^ In the low-frequency range. These two ranges 
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are separated by a distinct spectral peak. Kalmal (43) suggests that 
under stable conditions in the absence of appreciable low-frequency 
fluctuations all spectra can be brought in coincidence and .approximated 
by the empirical relation 


nS.(n) 0.164(f/f ) 

1 o 

o,^ l+0.164(f/f 

1 o 


(15) 


with i=u, V and w, 

where f»nz/U is the reduced frequency and f is the reduced frequency 

o 

at the point* of intersection of the extrapolation of the Inertial 

2 

subrange of the spectra, and the line nS^(n)/o^ =1. Kaimal’s spectra and 
variances are obtained from data in the frequency range 0.005<n£l0 Hz. 

The Von Karman spectrum functions (13« 14) can be modified into 
similar expressions 

nS (n) 0.156f/f 

U^ _ O /i ft 

9 9 S 

a ^ [l+0.108(f/f 

u o 

and 


nS^(n) 



0.12(f/y 


l+0.679(f/f )^ 

• o 

{l+0.255(f/f 

o 


( 17 ) 


with i=v and w. 

In this set of equations the reduced frequency is defined as f=nL^/U. 
However the parameter f/f^ for either the Kaimal or the Von Karman 
expressions represents the wavelength ratio where is the 

wavelength associated with the reduced frequency, f^. In both spectral 
expressions the parameter f/f^ is independent of the length scale which 
is either the elevation or the turbulence integral scale L^. 


The Kainal (15) and the Von (Carman (16,17) spectrum functions are 
nearly Identical in the inertial subrange, but the Kaimal spectrum 
predicts a slightly smaller spectral peaH> In the low-frequency range 
the u-spectra again are abotit identical, but the Von Karman v and w 
spectra fall slightly below the Kaimal spectrum^ 

In Figures 59, 60 and 61 the normalized logarithmic u, v and w 
2 

spectra (nS^(n)/o^ ) are plotted as a function of the modified reduced- 
frequency The spectra were taken from different data records for 

winds from southerly directions, which were classified according to the 
local stability parameter z/L. The spectral data and the variances 
were obtained from data analyzed in the high-frequency range 0.02A4<n<100 Hz. 
The velocity spectra obtained from stable-air records above the IBL do 

•i 

not differ from those obtained in the unstable air in the IBL and all 
fit the Von Karman and the Kaimal spectral functions remarkably well. 

The empirical spectrum functions (15,16,17) for estimation of the 
velocity spectra in the case low-frequency fluctuations are absent can 
be extremely useful if values of f^ can be predicted. Based on the 
experimental results it is obvious that f^ varies with height and with 
the presence or absence of appreciable low-frequency velocity fluctuations. 
The results did not indicate any systematic variation of f^ with stability 
as suggested by Kaimal [43]. Averaged values of f^=(nz/U)^ for each 
velocity component and obtained from normalized logarithmic spectra in 
the high-frequency range are shown as a function of height in Figure 62. 
However, values for f^ obtained from the same data records but analyzed 
In the middle- frequency range (0.00153<n^6.25 Hz) depend greatly on the 
presence of low-frequency velocity fluctuations. If no appreciable low- 
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frequency fluctuations are present, the values f^ are independent of the 
frequency range the data are analyzed. 

For all u-spectra investigated low-frequency fluctuations are 
present between the low cut-off frequencies for the middle- and high- 
frequency ranges, 0.00153 Hz and 0.02A4 Hz respectively. These spectra 
do; not exhibit a distinct spectral peak but instead values of the 
logarithmic spectra are approximately the same for f<0.1. Values of 
f^ obtained from the u-spectra analyzed in the middle-frequency range 
are generally smaller than those obtained from the same data records 
but analyzed in the high-frequency range and seem to converge to a 
general value in the range^ 0.02<f^<0.03, Independent of height. 

For those records where Internal gravity waves affect the spectra 
above a frequency of n»0. 00153 Hz, the frequency range for which 
spectral values are increased varies with elevation. At the lowest 
elevation (z*9. 1 m) only the spectral values at the lowest frequencies 
are affected and an appreciable range where the normalized logarithmic 
spectrum varies as f'*’^ is still present (Fig. 63). However for spectra 
from higher elevations the f*^^-range becomes gradually smaller as the 
effect of the waves is felt at increasin'; frequencies until no appreciable 
frequency range with a f*^^ spectral distribution is present (Fig. 64). 

For those cases values of f^ for the u-spectra seem to vary between the 
lower limit of . f *0.02-0.03 and the values of f obtained in the hlgh- 
frequency range as shown in Figure 62. Similar observations can be made 
for the V and w spectra. The values of f^ obtained from spectra analyzed 
in the middle-frequency range are generally lower than those obtained 
from the same spectra analyzed in the high-frequency range. 
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Spectra obtained from data records for westerly winds exhibit 
significant low- frequency content for all elevations and the spectral 
data fit the Kalmal or the Von Karman spectrum functions (15,16,17) in 
the inertial subrange (-2/3 region) only (Figs. 65,66,67). The v-spectra 
(Fig. 66) show a peculiar shape which is typical for spectra of the 
lateral velocity component in the surface layer of a convective boundary 
layer [19]. Kaimal’s explanation for this shape is based on the fact 
that in the inertial subrange the spectral values of the v-component 
are 4/3 times larger than the u-spectra, a requirement for isotropy in this 
range. (A similar situation exits for the w-spectra.) As the w-spectra 
reach their peak and start to roll off with lower frequencies, the v-spectra 
instead continue to increase and start to follow the u-spectra. The 
result of this is the peculiar shape of the v-spectra in the transition 
between the -2/3 range and the low-frequency range where the u and v spectra 
both are Independent with elevation but Instead vary with the height, z^, 
of the convective boundary layer. 

In the case the low-frequency fluctuations are absent as is the case 
just before sunset, as the convective boundary layer disintegrates rapidly, 
the u, v and w spectra (Figs. 68,69,70) and specifically the v-spectra 
(Fig. 69) have a completely different character. The v-spectra show a 
distinct spectral peak and a rapid roll-off at lower frequencies although 
the spectra values fall above the Von Karman prediction in this range. 
However in comparison with the spectra of the daytime run 19 (figs. 65, 
66,67], the spectra of the evening run 16 show much lower spectra values 
in the low-frequency range and the peculiar shape of the daytime v-spectra 
as discussed above has disappeared and the v-spectra resemble the Von Karman 


prediction. 

Values of the reduced frequency^ defined as the intersection 

of the extrapolation of the spectra in the Inertial subrange and the line 

nS^(n)/a^ >1, were obtained from spectra analyzed in both the high-^ 

and middle-frequency range. For the daytime u and v-spectra for which 

appreciable low-* frequency velocity fluctuations aCu px'esent, the values 

of f^ obtained in the high-frequency range are systematically three 

times as large as the corresponding values obtained from the spectra 

analyzed in the middle frequency range, while for the w-spectra the 

ratio (f )up/(f )^n, is approximately 1.6, Indicating that the low- 

frequency content is larger in the u and v spectra than in the w-spectra. 

For run 16 this ratio for the u-spectra is 1.9 and for the v and v spectra 

1.25. These results clearly indicate the effect of the low-frequency 

fluctuations on the location of the Inertial subrange when the spectra 

are presented in the logarithmic form with nS^(n)/o^ and f*nz/U as 

coordinates. The distribution of f^ with elevation for daytime spectra 

and evening spectra analyzed in the middle-frequency range are shown 

in Figure 71. The results from the daytime spectra show that the values 

of f^ for each velocity component are approximately independent with 

height and are (f ) =^0.01, (f ) ,-0.02 and (f ) «0.07. This observation 
o U O V o w 

is in agreement with some of the results obtained for southerly winds 
although the values are somewhat higher because of less low-frequency 
content. If no low-frequency fluctuations are present in the velocity 
components, the values of f^ generally increase with height (Figs. 62,71). 

The values of f^ can be used to obtain the wavelengths corresponding 
to the spectral peaks associated with either the Kaimal or the Von Karman 
spectral functions. Since for both empirical relations the logarithmic 
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spectral peak Is approximately located at f/£^~3.8, then £^-3.8f^ and 

(X_)^ ■ 0,26z/(£ ). £or t»u,v,w. (18) 

ml o 1 

Here the wavelength, peak of the Kalmal 

and Von Karman spectrum functions and to the peak of the measured spectra 
if no appreciable low-frequency fluctuations are present in the data 
records. The peak wavelengths, (^,„)i» often used in micrometeorolpgy 
as measures of the energy containing eddies, or they can be slightly 
modified to fit the original Von Karman spectrum functions (13,14) 
from which values of can be predicted as proposed by Teunlssen [25] 

= 0.146 (X ) 
u mu 

L* = jC).106 (X ) . for i=u,W (19) 

X ; m 1 

However it must be realized that these predicted scales associated with 
the empirical spectrum functions are equivalent to .the turbulence 
Integral scale obtained from correlation functions only If no large 
low-frequency velocity fluctuations are present. In the case large 
low-frequency velocity fluctuations are present and not filtered from 
the data records, the scales obtained from correlation functions are 
generally larger In magnitude (see section 6.4). 
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7. SUMMARY AND CONCLUSIONS 

A rasbiet' detailed description has been given of a micrometeoro- 
loglcal facility consisting of an instrumented 76 m (250 ft) tower 
located within a 100 m distance from the shore at Wallops Island, 
Virginia. The instrumentation system consists of cup-vane and 
temperature Instruments mainly used for profile measurements and a 
hot-film system for turbulence measurements. The data acquisition 
and handling system for the hot-tllm instruments is located in an 
instrumentation trailer located at the base of the tower. The heart 
of this system is a PDF 11/20 DEC minicomputer which controls the 
digitization of the data (200 Hz sample rate) and the data transfer 
onto digital tiipe. The digitized data have been analyzed on an IBM- 
370 computer located on the VPI and SU campus. 

Data have been acquired with the cup-vane system under moderately 
strong wind conditions for all wind directions during a 4 1/2-year 
period. Prom this data-base mean velocity and mean temperature pro- 
files and associated parameters (roughness length, z^, and powerlaw 

exponent, a) have been derived as well as turbulence intensities, 

z z 

a /V and 0 /U, and turbulence Integral scales, L , L and L . 
u V . u u v 

Averages of the calculated flow parameters from all data records in 
each 15-degree sector have been presented. In addition averaged mean 
velocity ratios intensities, o^^/U, 0^/U,and 

turbulence integral scales, L^, have been obtained for 11 
sectors each with near-uniform upstream terrain. The results provide 
information about the microclimate at this site under moderately strong 
wind conditions. This information is graphically presented in 
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Figures 72 through 75 from which average wind design data for this 
coastal site can be established. In addition , data have been ac- 
quired with the hot-film system for the southerly and northwesterly 
prevailing wind directions. With this system turbulence parameters 
such as turbulence intensities, Reynolds stresses, turbulence heat 
fluxes. Integral scales (L^, L* and L*) and power spectra of the 
three velocity components have been obtained. 

For all observations made at this site under moderately strong 
wind conditions, truly neutral thermal stratifications have never 
been encountered throughout the observation height of 76.2 m for 
any length of time. For westerly wind direction under sunny daytime 
conditions the measured velocity and temperature profiles suggest 
that the surface flow at the Wallops Island site is similar to the 
surface flow observed during the Minnesota experiment [12]. The 
observed PBL flow at Minnesota is an example of a typical convective 
boundary layer, a model of which is described in detail 1ft Reference 
13. In addition to mechanical and convective turbulence generated in 
this atmospheric boundary layer, large-scale turbulence due to the 
interaction of the mixed layer and the capping inversion (entrainment) 
affects the mean and turbulent surface flow regardless of the wind 
velocity. However, appreciable deviations from the. convective 
boundary-layer model may occur depending on atmospheric conditions, 
time of day and wind direction. 

It has been observed that just before sunset the daytime 
boundary-layer flow is modified drastically as a result of the 
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disappearance of the large-scale turbulence and appearance of a 
surface-based inversion. Similarly under conditions of low cloud 
cover combined with precipitation, an inversion below the cloud 
cover nuiy develop, impeding the regular development of the day- 
time convective boundary layer. Under these conditions the influ- 
ence of the large-scale turbulent motions on the flow below the in- 
version is reduced, resulting in an appreciable reduction in tur- 
bulence intensities and turbulence integral scales . large negative 
lateral velocity fluctuations or large wind direction fluctuations 
towards the north have been observed for northwesterly wind direc- 
tions specii-ically between *300° and 315^. For winds in this sector 
larger values of the lateral turbulence intensity and larger tur- 
bulence integral scales have been observed than for winds outside 
this sector but with the same upstream terrain. In addition, tur- 
bulence intensities and turbulence integral scales vary during the 
daytime as the convective boundary layer develops. The above 
observations have been made under moderately strong wind conditions 
with hourly mean-wind speeds between 10 m/s and 20 m/s at z=9.1 m. 
Based on all the observations made for westerly winds at Wallops 
Island, there is no evidence that similar flow variations in the 
surface layer would not exist under extreme and potentially damaging 
wind conditions with velocities in excess of 20 m/s. Consequently 
at this point in time it cannot automatically be assumed that for 
extremely strong winds from westerly directions, the PBL flow at 
the Wallops Island site is similar to the purely shear-generated. 
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neutrally-stratified boundary-layer flow model, which is so often 
advocated by wind engineers. 

The mean and turbulent flow for southerly winds also differs appre- 
ciably from that predicted by the neutral boundary-layer model. During 
the summertime the warm air blowing over the cooler ocean water gives 
rise to a surface-based inversion of variable height. Depending 
on the thermal stability, a low-level jet with a maximum velocity 
occasionally below the highest observation level has been observed. 

Under extreme stable conditions the turbulence at the two highest 
observation levels has been observed to vanish completely and 
generally internal gravity waves may co-exist with the turbulence. 

Under these conditions 1 the surface layer is very shallow, well below 
the lowest observation level. Moreover, the flow near the surface 
will also undergo a modification as soon as the ocean air crosses 
the beach and experiences an Increase in surface roughness and 
surface temperature. These modifications of the surface flow 
manifest themselves in the form of a developing Internal boundary 
layer (IBL) which at the tower location is between 15 m and 30 m in 
height, depending on the change in surface temperature and the 
• overland development distance which varies with wind direction. 

The conclusions of the boundary-layer experiment at Wallops 
Island can be summarized as follows: 

I. Westerly wind directions 

1. The observed daytime flow below 76m at Wallops Island is 
described better by the convective boundary-layer model 
[13] than by the neutral boundary-layer model. 
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2. The height at which transition occurs from the logarithmic 
velocity profile to the mixed-layer velocity profile 
varies with wind velocity surface roughness and thermal 
stability. 

3. The roughness length, z , obtained from velocity profiles 
below the transition are in agreement with predicted values 
from PBL-flow review papers. 


4. The panofsky relation, a=z /In is only useful for 

predicting values of powerSaw exponents for velocity profiles 
in the surface layer below the elevation where transition to 
the mixed-layer profile starts. 

5. Measured turbulence intensities are generally in agreement 
with predicted values, except for the lateral turbulence 
intensities, o /U, in the northwesterly wind-direction sector 
300"<4i<315'“ , Yn this sector the turbulence intensities of 
the horizontal velocity components (u and v) are of the same 
magnitude. 

♦ 

6. Turbulence intensities of the horizontal components also 
vary with time of day and atmospheric conditions, or in 
general with the absence or presence of appreciable low- 
frequency velocities fluctuating in the frequency range 
between 0.0015 Hz and 0.02 Hz. 

7. The magnitude of the turbulence Integral scales depends 
on the method (direct method or spectral method) by 
which they are calculated and also on the presence or 
absence of appreciable low-frequency velocity fluctuations. 

8. If appreciable low-frequency content is present and is not 
filtered from the data records, the turbulence Integral 
scales obtained via the direct method are larger than the 
predicted values. 

9. The turbulence Integral scales vary also with time of day, 
wind direction, surface roughness, and atmospheric 
conditions such as cloud cover combined with precipitation. 

z z 

10. The measured vertical integral scales, L and L vary with 
direction of Integration but are generality In agreement with 
predicted values. 


(>7 


11. The ratio varies Benornlly betweert 3 and A and the 

5T 2 

ratio has an approximate value of 2. 

12. The daytime turbulence spectra follow the Kaimal model 
[19] and deviate appreciably from the Von Karman model 
especially In the low-frequency range. 

13. The Von Karman spectral model does not fit the measured 
spectra if appreciable low-frequency velocity fluctuations 
are present and are included in the spectral analysis and 
in the calculation of the variance, and integral scale 
(direct method) . 


II . Southerly wind Directions 

1. For on-shore winds an IBb develops as the surface air passes 
tlie beach and experiences an Increase in surface rougitness 
and an Increase in surface temperature especially in the 
sunmier during .daytime. 

2. For southerly wind directions warmer air flows over the 
cooler water creating a surface based inversion which is 
characterised by a very shallow surface layer, low-level 
maximum velocity (surface jet) , low turbulence intensities, 
occasional vanishing of the turbulence under extreme stable 
conditions and co-existence of turbulence and low-frequency 
internal gravity waves. 

3. No simple boundary-layer flow model is available to describe 
the on-shore flow at the Wallops Island site. Variations 

in observed velocity and tempera tvire profiles, turbulence 
intensities and turbulence integral scales are extremely 
high and can occur within a very short time. 

A. Measured velocity spectra (excluding the low-frequency gravity 
waves) are independent of thermal stability and seem to fit 
the modified Von Karman spectrum model (16,17) and the Kaimal 
stable spectrum model (15) extremely well. 

As the above conclusions clearly Indicate, there is no single and no 

simple PBL-f low model available to describe the mean and turbulent 

flow near the surface under rooderately strong wind conditions at the 

Wallops Island site. The presence or absence of appreciable low- 

frequency velocity fluctuations causes the parameters describing 

this flow to vary a great deal. The non-uniform surface conditions 
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and the presence of the low-frequency fluctuations, either in the 
convective boundary layer or In the on-shore winds In the form of 
internal gravity waves, cause the observed surface flow to be quite 
different from the neutral boundary layer flow model • For the pre- 
vailing winds from either the south or westerly directions the 
experimental results do not show any evidence for the PBL flow to 
approach the neutral boundary-layer model as the wind speed Increases 
from moderately strong to extreme. 

It Is assumed that engineers, already aware of uncertainty In 
modeling the PBL flow, use safety factors In the design of wind 
turbines to allow for differences in the actual wind environment 
in comparison with the* predictions from the neutral PBL model. The 
variation of the turbulence intensities and turbulence Integral scales, 
measured under moderately strong wind conditions at the Wallops site, 
is appreciable. Consequently a great deal of difference may exist 
between actual measurements and the neutral PBL model. Experimental 
evidence does not Indicate that mean wind and turbulence parameters 
will conform closer to the neutral PBL-model under higher wind-speed 
and slightly unstable conditions. The observed differences appear at 
the lower frequencies which are pertinent to the response characteristics 
of the larger machines. Therefore, the design of large wind turbines 
may need an increased safety factor with respect to turbulence at 
frequencies below about 0.01 Hz. 
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Figure 7b. Variation of Mean-Velocity Ratio, V(30.5)/V(76.2) with Wind Direction 
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Figure 7d. Variation of Mean'-Velocity Ratio, V(61.0)/V(76.2> with Wind Direction 
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Kigure 9. Typical Strong-Wind Profiles for Westerly Wind Directions 
Gradient Richardson Number (5) Evaluated at z-15.2m. 
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Figure 11. Variation of the Roughness Length* Zq* with the Gradient 
Richardson Number (5) Evaluated at z=15.2m 



PANOFSKY [31] 
COUNIHAN [26] 
ESDU [27] 


89 



o ‘lN3NOdX3 MV1 d3A/\0d 


Figure 12. Variation of Power-law Exponent with Roughness Length, Both Based on Velocity 

Measurements Below 45.7m. Panofsky's Relation (7j Based on 2.*15. 2m and 22*45. 7m 
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Figure 15a. Variation of Turbulence Intensity, cr /U, with Direction, z=15.2m 
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Figure 15d. Variation of Ti»rbulence Intensity, a /U, with Direction, z>61.0m 
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Figure 18. Typical Strong-Wind Aerovane Record of Direction and Speed 
for Northwest Winds, (^“320*, Uj^q= 13 m/s 
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Figure 22, Vertical Distribution 
Oy/U, for Southerly W 
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Figure 23. Vertical Distribution of the Turbulence Intensity, 
0 /U, .for Southerly Winds 
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Figure 24. Vertical Distribution of the Turbulence Intensity, ff./U 
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Figure 27b. Variation of Turbulence Integral Scale, L +, with Wind Direction z*30.5m 
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Figure 28d. Variation of turbulence Integral Scale, L^|, with Wind Direction 2 « 45 . 7 m 
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Figure 34. Variation of kvexaze^ Integral Scales, L , (Cup-Vane) with Height for Wind-Direction 
Sector 180*-195* (South Winds) z<20ar+2 J^O.Ola and z>20m*z "O.OOlffl 













Variation of Average Integral Scales, Ly, (Cup-Vane) with Height for Wind* 
Direction Sector 180®-195® (South Winds) z<20nr^z -O.Olra and z>20nr>"20m=0.0i 
Arrows on Si^nbols Indicate the Direction of Integration 
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Variation of the Integral Scale, at z=l5.2m (50 ft.) with Roughness Length 
z , for all Data Records with a Westerly Wind Direction and Rij^^>-0.1 
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Figure 43. Logarithmic u-Spectra for Run 19, o and L Obtained in the High Frequency 
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Figure 48. Variation of Average Integral Scales, Lu» Obtained with the Hot-Filo 
System (Von Karman Method, High and Middle Frequency Range), with 
Height for Southerly Winds 
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Figure 49. Variation of Average Integral Scales, L^, Obtained with the Hot-Film System 
(Von Karman and Direct Methods, Middle-Frequency Range) and Cup-Vane System 
with Height for "Northwesterly Winds 
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RECORD OF STREAmWiSE VELOCITY COMPONENT, 
RUN 19, 250 ftrLEVEL, U*I3.9 m/sec , (Ty =1.60 m/sec, 
3‘25-'77, 14:21:00, lJ =306m 
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RECORD OF LATERAL VELOCITY COMPONENT, 
RUN 19, 250fl.“LEVEL, Oy =|.65m/sec, Ly*479m 





ir|n,r'll' II 


RECORD OF STREAMWISE VELOCITY COMPONENT, 
RUN 16, 250 ft.-LEVEL,U = 11.4 m/sec, av =0.81 m/sec, 
3-23-77, 20:30:00 , lJ =77.9m 




RECORD OF LATERAL VELOCITY COMPONENT, 
RUN 16, 250 ft.- LEVEL, (Ty =0.62 m/sec, Lj=23.3m 


Figure 54. Time Records of Velocity Components. Elapsed Time 55 Minute 
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Figure 55. Comparison of Integral Scales, Lv, of Record #16 (Evening Run) 

(Von Karman and Direct Methods, Middle Frequency Range) Against 
Average of Day Time Hot-Film Runs (Von Karman, Middle Frequency Range) 
for Northwesterly Winds 
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Comparison of Integral Scales, L , Obtained from the Hot-Film System and 
Analyzed in the Middle^and High-Frequency Range Using the Von Karman Method 
for Northwesterly Winds 




Figure 57. Comparison of Integral Scales, L^,, of Record ^^16 (Evening Run) 

(Von Karman and Direct Methods, Middle Frequency Range) Against 
Average of Day Time Hot-Film Runs (Von Karman, Middle Frequency Range) 
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Figure 59b. Logarithmic u-Spectra Versus Modified Reduced 


Frequency, f/f , Analyzed in the High-Frequency 
Range, South W?nds, Stable Thermal Stratification. 
0<z/L<+1.0 
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Figure 60b. 


ogarithmic v-Spectra Versus Modified Reduced 
requency, f/f_. Analyzed in 

ange, South w2nds, Stable Thermal Stratification 
/T<r4-1 .0 









ORIQINAl. PAGE (3 
OF POOR QUAUmr 







167 


ORIGINAL PAGE !S 
OF POOR QUALITY 



Figure dlb. Logarithmic w-Spectra Versus Modified Reduced 

Frequency, f/r. Analyzed in the High-Frequency 
Range, South Winds, Stable Thennal Stratification, 
0<z/L<+1.0 









Figure 62, The Vertical Distribution of the Reduced Frequent 
with i=u,v,w for Southerly Winds, Data Analyzed 
the High-Frequency Range 0.0244<n<100 Hz. 
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Figure 66. Logarithmic v-Spectra Versus Modified Reduced Frequency, Analyzed in the 

High and Middle Frequency Range, Northwest Wind, <t>*301“ 
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Logarithmic v-Spectra Versus Modified Reduced Frequency, f/f » Analyzed 
in the High and Middle Frequency Range, Evening Run, Northwest Wind (j>*292 
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Figure 71. The vertical Distribution of the Reduced Frequency (f^) 
With i*u»V»w for Westerly Winds. Data Analyzed in 
the Middle- Frequency Range 0.0015<n<6.25 Hz. 

Dark Symbols Denote the Data for the Evening Run 16. 
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Note: Type of Profile: S-single logarithmic profile, K--kink in profile. 


